D61 :: BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007792	1,817	0	C	G	70.3%	218.1 / 166.5	195	S213S (TCC→TCG)	SAUSA300_RS14745	integrase
Reads supporting (aligned to +/- strand): ref base C (20/38); new base G (17/120); total (37/158)
Fisher's exact test for biased strand distribution p-value = 5.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 209 total aligned reads displayed.

AGTATCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCAAAT  >  NC_007792/1722‑1915

aGTATCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCag                                                                                                <  1:3197881/100‑2 (MQ=40)

aGTATCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCag                                                                                                <  1:1227012/100‑2 (MQ=40)

   aTCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGcc                                                                                             <  1:3888465/100‑1 (MQ=40)

     cAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATgggct                                                                                           <  1:1000375/100‑5 (MQ=255)

      aaTACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCaggggcct                                                                                          <  1:1715129/100‑8 (MQ=255)

      aaTACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACg                                                                                          <  1:4883870/100‑1 (MQ=40)

              aaaGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACGAAATCAGc                                                                                  <  1:3950819/100‑1 (MQ=25)

                            caGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAg                                                                    >  1:3596100/2‑100 (MQ=25)

                                 aTTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTaaa                                                               >  1:919114/1‑100 (MQ=38)

                                     gctGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGaa                                                           >  1:2453066/3‑100 (MQ=255)

                                            ttACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTgct                                                    <  1:5028158/100‑3 (MQ=255)

                                                               aaCCGCAGGTCTCTTCAGATCT‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAg                                 >  1:3000123/1‑100 (MQ=12)

                                                                  cgcaGGTCTCTTCAGATCT‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCta                              <  1:1504807/97‑1 (MQ=12)

                                                                   gcaGGTCTCTTCAGATCT‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCtat                             <  1:4558901/98‑1 (MQ=18)

                                                                     tgtgctctTCCGATCT‑ACGGGTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCtaagt                           <  1:270982/96‑4 (MQ=11)

                                                                      gtgctctTCCGATCTGCCGGGTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCtatag                           <  1:4444384/97‑2 (MQ=11)

                                                                          tctTCAGATCT‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATAtttttt                      <  1:3423710/100‑1 (MQ=25)

                                                                                 ctCT‑GCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGc               <  1:3928420/99‑1 (MQ=255)

                                                                                 aTCT‑ATGGGTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGc               <  1:4724662/100‑1 (MQ=255)

                                                                                 aTCT‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGc               <  1:3648230/100‑1 (MQ=25)

                                                                                 aTCT‑ACGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGc               <  1:4223161/100‑1 (MQ=255)

                                                                                  tCT‑ACGGTTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCa              <  1:603283/100‑1 (MQ=255)

                                                                                  tCT‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCa              <  1:1337975/100‑1 (MQ=25)

                                                                                  tCT‑ACGGACTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCa              <  1:350312/100‑1 (MQ=255)

                                                                                   tctgggggttttCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCa              <  1:5009962/92‑1 (MQ=255)

                                                                                    t‑ACGGTTTTTCG‑CCATGCCGCGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:1458447/100‑1 (MQ=255)

                                                                                    t‑ACGGTTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:2257616/100‑1 (MQ=255)

                                                                                    t‑ACGGTTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4948345/100‑1 (MQ=255)

                                                                                    t‑ACGGCTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAGCACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:804651/100‑1 (MQ=255)

                                                                                    t‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4324784/100‑1 (MQ=25)

                                                                                    t‑ACGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:281432/100‑1 (MQ=25)

                                                                                    t‑ACGGACTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTCACAAc            <  1:4103548/100‑1 (MQ=11)

                                                                                    t‑ACAGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:236787/100‑1 (MQ=255)

                                                                                     ttccgATCTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:966690/96‑1 (MQ=255)

                                                                                     tgcGGGTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:3299137/98‑1 (MQ=255)

                                                                                     tatgtaTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4612139/95‑1 (MQ=255)

                                                                                     tatgggttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:2666766/94‑1 (MQ=255)

                                                                                     tatggATTTTCG‑CCGTGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4015152/97‑1 (MQ=255)

                                                                                     tatggATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4055316/97‑1 (MQ=255)

                                                                                     tatggATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4397824/97‑1 (MQ=255)

                                                                                     tatgaaTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:3363313/95‑1 (MQ=255)

                                                                                     tagggATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4909516/97‑1 (MQ=255)

                                                                                     tacgatttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:2034910/94‑1 (MQ=255)

                                                                                     ctgggATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4752323/97‑1 (MQ=255)

                                                                                     caCGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:3668325/99‑1 (MQ=25)

                                                                                     caCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:4862959/99‑1 (MQ=255)

                                                                                     acgcgATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:348318/96‑1 (MQ=255)

                                                                                     accGGTTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:2419429/98‑1 (MQ=255)

                                                                                     aaCGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAc            <  1:2212401/99‑1 (MQ=25)

                                                                                      aCGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTGCTTTGCAACa           >  1:493003/1‑100 (MQ=255)

                                                                                      aCGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:1648104/1‑100 (MQ=25)

                                                                                      aCGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:4187253/1‑100 (MQ=25)

                                                                                      aCGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTCGCAACa           >  1:90216/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGTAATTAGCATCATGATAGCAAGTTAAACGAACACTGACATGACATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:2627315/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:140979/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:2460869/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:1551240/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:2165926/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:3494096/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACa           >  1:4680309/1‑100 (MQ=255)

                                                                                      aCGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTCGCAACa           >  1:3122513/1‑100 (MQ=255)

                                                                                       tggATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:1167919/99‑1 (MQ=25)

                                                                                       tggATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:3163159/99‑1 (MQ=25)

                                                                                       gtgATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:3837309/98‑1 (MQ=255)

                                                                                       gggtttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:3283049/96‑1 (MQ=255)

                                                                                       cggtctttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:447939/95‑1 (MQ=255)

                                                                                       cGGTTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:4812979/100‑1 (MQ=255)

                                                                                       cGGCTTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:5101305/100‑1 (MQ=255)

                                                                                       cGGCTTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:2935940/100‑1 (MQ=255)

                                                                                       cGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:181004/100‑1 (MQ=25)

                                                                                       cGGATTTTCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAg          <  1:209217/100‑1 (MQ=25)

                                                                                       acgATTTTCCGCCATGCCACGAAATTAGCATCATGCTAGTAAGTTAATCGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAAcg           <  1:4707099/98‑2 (MQ=11)

                                                                                           ttttCG‑CTATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:1967971/100‑1 (MQ=255)

                                                                                           ttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:5194384/100‑1 (MQ=25)

                                                                                           ttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:852676/100‑1 (MQ=25)

                                                                                           ttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:1948624/100‑1 (MQ=25)

                                                                                           ttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:3411405/100‑1 (MQ=25)

                                                                                           ttttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:4516767/100‑1 (MQ=25)

                                                                                           ttttCG‑CCATGCCACGAAATCAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:3773723/100‑1 (MQ=255)

                                                                                           tcttCG‑CCGTGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:2697402/98‑1 (MQ=255)

                                                                                           ctttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:187518/99‑1 (MQ=25)

                                                                                           ctttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:810964/99‑1 (MQ=25)

                                                                                           ctctCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:669423/97‑1 (MQ=255)

                                                                                           atctCC‑CCGTGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAAcc      <  1:742624/97‑1 (MQ=255)

                                                                                            tttCG‑CCGTGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCa     <  1:1610037/100‑1 (MQ=255)

                                                                                            tttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCa     <  1:3012272/100‑1 (MQ=25)

                                                                                            tttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCa     >  1:2225363/1‑100 (MQ=25)

                                                                                            tttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCa     <  1:921177/100‑1 (MQ=25)

                                                                                            tttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCa     <  1:2290600/100‑1 (MQ=25)

                                                                                            tttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCa     <  1:759709/100‑1 (MQ=25)

                                                                                             ttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCat    <  1:2143132/100‑2 (MQ=25)

                                                                                             ttCG‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCat    <  1:1256807/100‑2 (MQ=25)

                                                                                             ttCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCag    <  1:2544736/100‑2 (MQ=255)

                                                                                             tgcc‑ccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCat    >  1:4140403/3‑99 (MQ=255)

                                                                                             tctc‑ccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCag    <  1:2922137/97‑2 (MQ=255)

                                                                                              tCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCagc   <  1:2077487/100‑3 (MQ=255)

                                                                                               tc‑ccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCAGAt  <  1:1751391/99‑1 (MQ=255)

                                                                                               cc‑ccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCagca  <  1:2766399/100‑4 (MQ=255)

                                                                                               cc‑ccATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTACAACAGatcggaag  >  1:1989832/1‑93 (MQ=255)

                                                                                               cc‑ccATGCCACGAAATTAGCATCATGCTAGAAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCagcc  >  1:4993488/1‑97 (MQ=255)

AGTATCAATACGGCAAAGAATACTTTAAAAGGTATTGAATGTATTTACGCTCTATATAAAAAGAACCGTAGGTCTCTTCAGATCT‑ACGGATTTTCC‑CCATGCCACGAAATTAGCATCATGCTAGCAAGTTAAACGAACACTGACATGATATATTAGTGATTAGCTATATTTTTTACTTTGCAACAGAACCAAAT  >  NC_007792/1722‑1915

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`