BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,351,228	G→C	13.9%	intergenic (‑189/+179)	fabI ← / ← ycjD	enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,351,228	0	G	C	13.9%	634.8 / 29.6	213	intergenic (‑189/+179)	fabI/ycjD	enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand): ref base G (88/92); new base C (28/4); total (117/96)
Fisher's exact test for biased strand distribution p-value = 3.33e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

TCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  NC_000913/1351080‑1351376

   ggCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                                                                                                                                      <  1:3791325/150‑1 (MQ=255)

   ggCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCAAACCCTCTCCCCACAggg                                                                                                                                                      >  2:2973916/1‑150 (MQ=255)

       aTTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAAGCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                                                                                                                                  <  2:3641857/150‑1 (MQ=255)

        ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAg                                                                                                                                                 <  2:3781044/150‑1 (MQ=255)

         tcgtcgTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgg                                                                                                                                                <  1:3595973/150‑1 (MQ=255)

         tcgtcgTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                <  1:2817770/150‑1 (MQ=255)

             cgtTTTAGTACACAGTACGAACAGATAAACGGTTATTAGAATCAACCTGGCTGTGAGTAGCTATAGGTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑c                                                                                                                                          <  2:3966791/150‑1 (MQ=255)

             cgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑c                                                                                                                                          <  2:892410/150‑1 (MQ=255)

             cgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCGCCCCTCTGGGGAGAgggtta                                                                                                                                            <  1:941040/150‑4 (MQ=255)

              gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACCCCCGAAACCTCTCCCCAGAGGGGCGAGGGG‑c‑cc                                                                                                                                         >  2:3532951/1‑150 (MQ=255)

               ttttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCg                                                                                                                                        >  1:1671835/1‑150 (MQ=255)

                tttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGt                                                                                                                                       <  2:2951440/150‑1 (MQ=255)

                tttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACAACCAAAACCACACCCCACATGGGCGAGGGG‑A‑ccct                                                                                                                                       >  2:219878/1‑148 (MQ=255)

                 ttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCACACCCAAAACCAAGCCCCAGAGGGGCGAGGGG‑G‑ccct                                                                                                                                       >  2:1362291/1‑148 (MQ=255)

                 ttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTACCCACACCCAAAACCACACCCCACAGGGGCGATGGG‑A‑CCga                                                                                                                                       >  2:1362352/1‑149 (MQ=255)

                  tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTAt                                                                                                                                     >  2:554179/1‑150 (MQ=255)

                    gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCTCACCCAAAACCACTCCCCAGAGGGGCGAGGGG‑G‑CCGTATtg                                                                                                                                   >  2:782478/1‑150 (MQ=255)

                      aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATtgt‑g                                                                                                                                <  1:2993992/150‑1 (MQ=255)

                        aCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCGAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCa                                                                                                                              >  1:1368540/1‑150 (MQ=255)

                           gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAACCCTCACCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAat                                                                                                                           >  1:4076377/1‑150 (MQ=255)

                            tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAata                                                                                                                          >  2:4069148/1‑150 (MQ=255)

                            tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCAAAACCGCTCCCCCGAGGGGCGAGGGG‑G‑CCGGGTTGT‑GGCAata                                                                                                                          >  1:2904906/1‑150 (MQ=255)

                            tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTACCCACACCCTAAACCAAACCCCACACGGGCGAGGGG‑A‑CCCTATTGT‑GCAAata                                                                                                                          >  1:2981853/1‑150 (MQ=255)

                               gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAACCCACTCCCCACAGGGGCGCGGGG‑A‑CCGGATTGT‑GGCACTATTg                                                                                                                       >  1:1295789/1‑150 (MQ=255)

                               gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAACCCACACCCCACAGGGGCGAGGGG‑A‑CCCTATTGT‑GCaaaaaatt                                                                                                                       >  2:113603/1‑145 (MQ=255)

                               gAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCAAACCCACACCCCCGAGGGGCGAGGGG‑A‑CCCGATTGT‑GCAAATAttt                                                                                                                       >  2:374627/1‑149 (MQ=255)

                                aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTACCCACACCCTAACCCTCACCCCAGAGGGGCGAGGGG‑A‑CCGGATTTT‑TGAAATATTG‑t                                                                                                                     >  1:3060637/1‑150 (MQ=39)

                                 aCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAACCCACTACCCACAGGGGCGAGGGG‑A‑CCGTATTGT‑TCAAAAATtt‑tt                                                                                                                    >  1:3917513/1‑150 (MQ=255)

                                  cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTa                                                                                                                   <  2:1896859/150‑1 (MQ=255)

                                     aTAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTAccc                                                                                                                <  2:940681/150‑1 (MQ=255)

                                     aTAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCACTCCCCGGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTAccc                                                                                                                >  1:2625500/1‑150 (MQ=255)

                                     aTAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAACCCACAACCCAGAGGGGCGAGGGG‑T‑CCGTATTGT‑GCAAATATTT‑TTAccc                                                                                                                >  1:1729464/1‑150 (MQ=255)

                                       aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCa                                                                                                              <  2:239606/150‑1 (MQ=255)

                                           ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAACCCAAACCCCAGAGGGGCGAGGGG‑G‑CCCGATTTT‑GCAAATAATGTTTAccccccaa                                                                                                           >  1:3849810/1‑146 (MQ=255)

                                            gTTATTATAATCAACCTGGCTGAGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTATCCCCACACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTCATGT‑GCAAATATTG‑TATCCCCCGccac                                                                                                         >  2:1265109/1‑150 (MQ=38)

                                               attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCaacaaa                                                                                                      <  1:2215017/150‑1 (MQ=255)

                                                ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCACTCCCCCGAGGGGCGAGGGG‑A‑CCCTATTGT‑GCAAAAATTG‑TTAACCCAGCaacaaac                                                                                                     >  2:416336/1‑150 (MQ=39)

                                                ttatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAACCCACACCCCACAGGGGCGAGGGG‑A‑CCGTAATGT‑GCAAATATTG‑TTAACCCAGCaacaaac                                                                                                     >  1:3419816/1‑150 (MQ=255)

                                                 tataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCaacaaaca                                                                                                    <  1:1249679/150‑1 (MQ=255)

                                                  ataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAg                                                                                                   <  1:2549908/150‑1 (MQ=255)

                                                  ataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCAAACCCTCACCCCACAGGGGCCAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCCGCAACAAACAg                                                                                                   >  2:243929/1‑150 (MQ=255)

                                                   taATCCACCTGGCTGTGAGTAGCGATAGGTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAgg                                                                                                  <  2:3483613/150‑1 (MQ=255)

                                                   taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAgg                                                                                                  <  2:576822/150‑1 (MQ=255)

                                                   taATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAAACCACAACCCCGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATAATG‑TTAACCCAGCAACAAAAAgg                                                                                                  >  2:811907/1‑150 (MQ=255)

                                                    aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCTCACCCAAAACCACACCCCAGAGGGGGGAGGGG‑G‑CCGGATTGT‑GCAAATATTT‑TTTCCCCCGCaaccaaaaagc                                                                                                 >  1:3320259/1‑145 (MQ=255)

                                                    aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGc                                                                                                 >  2:2371720/1‑150 (MQ=255)

                                                    aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCGTTTCCCCACACCCAAAACCGCACCCCACACGGGCGCGGGG‑G‑CCCGATTGG‑GCAAAAATTT‑TTAACCCAGCaaccaaaccgc                                                                                                 >  1:2146006/1‑145 (MQ=255)

                                                     aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCGAACCCACTCCCCACGGGGGCGCGGGG‑G‑CCCGATTGT‑GCAAATATTG‑TTACCCCACCAACAAACACGCt                                                                                                >  2:1086123/1‑150 (MQ=255)

                                                         aCCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATa                                                                                            <  2:1672195/150‑1 (MQ=255)

                                                           cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCACACCCAAAACCACGCCCCAGAGGGGGGGGGGG‑G‑CCCTATTGT‑GGAAAAATTT‑TTAAACCACCAAAAAAAAGGCTCaatca                                                                                          >  2:726040/1‑146 (MQ=255)

                                                                   gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCGCAACCAAACCCACAACCCCCAGGGGCCCGGGG‑ACCCACATTTT‑G‑AAATATTT‑TTAAACCACCAACAAACACGCGAAAAAAGCCCCaaa                                                                                  >  1:102902/1‑150 (MQ=255)

                                                                   gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTaa                                                                                  <  2:2915507/150‑1 (MQ=255)

                                                                     gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCGCCCCCAAACCCACACCACCGAGGGGCGAGGGGCC‑CCAAATTGT‑‑CAAAAAATT‑TTTCCCCACCAACAAAAAAGCTCAAACACCCCCAAAca                                                                                >  2:354192/1‑149 (MQ=255)

                                                                     gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTGCCCCCACCCAAACCCACACCCCACAGGGGCGAGGGG‑G‑CCGCAATGT‑GGCAAAATTT‑TTAACCCCCCACCAAACACGCGCAAACACCCCCAAAcc                                                                                >  1:1154883/1‑150 (MQ=255)

                                                                     gTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCGCACCCAAAACCAAAACCCAGAGGGGCGAGGGG‑A‑CCCTAATTT‑TCAAAAAATT‑TTAACCCACCAAAAAAAAGGCACAAACACCCCCaaaac                                                                                >  2:3876761/1‑148 (MQ=255)

                                                                        gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCACCAAACAGGCTCATACAGCCCCAAACCCtg                                                                             >  2:3075681/1‑149 (MQ=255)

                                                                         cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCAAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCGCCAACAGCCCCAAACCCtgt                                                                            >  2:2294420/1‑148 (MQ=38)

                                                                         cTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCAAACCCACTACCCCCAGCGGCCGCGGG‑A‑ACCTATTGG‑GCCAATATTT‑TTACCCCAGAAACAAACAGCCGCAAAACGCCCCGAAcccgta                                                                            >  2:26075/1‑147 (MQ=255)

                                                                         ataaaGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCTAACCAACACCCCAGAGGGGCGAGGGG‑A‑CCGTAATGT‑GCAAATATTT‑TTAACCCAGCAAAAAACAGGCACAAAAAGCCCCAAACACtat                                                                            >  2:2678523/5‑148 (MQ=255)

                                                                          tataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCCCCCCCAAAACCACACCCCCGAGGGGCGGGGGG‑C‑CCCGGTTGT‑GGAAAAGATG‑TTAACCCACCAACAAACACGCGCAGACACCCCCAAAaccaata                                                                           >  1:2294344/1‑146 (MQ=255)

                                                                           ataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTACCCCCACCCAAAACCACAACCCCCAGGGGCGCGGGG‑G‑CCCGATTTT‑GGAAAAATTT‑TTAACACACCAACAAAAACACGCAAACAGCCCCAAAcccaataa                                                                          >  1:3305462/1‑145 (MQ=255)

                                                                             aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                                                                        <  2:2748022/150‑1 (MQ=255)

                                                                             aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                                                                        <  1:851817/150‑1 (MQ=255)

                                                                             aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                                                                        <  1:1860106/150‑1 (MQ=255)

                                                                               ttGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCACTCCCCAGAGGGGCGAGGGG‑A‑CCCTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGc                                                                      >  2:3347991/1‑150 (MQ=255)

                                                                               ttGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCTAACCCACTCCCCAGAGGGGCGAGGGG‑C‑CCCTATTGT‑GCAAAAATTG‑TTACCCCAGCAACAAACAGGCTCAAAAAGCCCCAAACCCGTTCAAGGc                                                                      >  2:1625782/1‑150 (MQ=255)

                                                                                 gCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGa                                                                    <  1:4108994/150‑1 (MQ=255)

                                                                                 gCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGa                                                                    <  1:1416501/150‑1 (MQ=255)

                                                                                    aGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCACACCCAAAACCAAGACCCACACGGGCGAGGGG‑G‑CCCTATTTTGGCAAAATTTG‑TTA‑CCCAGCAAAAAAAAGGGGCCAAAAACCCCAAAACCCTTAATTGCGATgg                                                                 >  2:1051743/1‑150 (MQ=255)

                                                                                      gTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCt                                                               <  2:25157/150‑1 (MQ=255)

                                                                                          gACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTggg                                                            <  1:1362302/150‑1 (MQ=255)

                                                                                                         cGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGc                                            <  2:1421013/150‑1 (MQ=255)

                                                                                                          gggatgggggggCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGc                                            <  1:1971084/145‑1 (MQ=255)

                                                                                                            cA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGa                                         <  2:2879071/150‑1 (MQ=255)

                                                                                                             a‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGaa                                        <  1:2215201/150‑1 (MQ=255)

                                                                                                             a‑GGGGGGGCGCTTTTCCCCACACCCGAACCCACAACCCCGAGGGGCCAGGGG‑G‑CCCTATTTT‑GCCAATATTT‑TTAACCCAGCAACAAACACGCTAAAAAACCCCCAAAACCTTTAATTTCGATGGCTTGGACGGTTCAGACCGTGCCGaa                                        >  2:1940900/1‑150 (MQ=255)

                                                                                                                     gcgcTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATtctc                                 <  2:2711294/150‑1 (MQ=255)

                                                                                                                      cgcTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTcc                                >  2:2851330/1‑150 (MQ=255)

                                                                                                                       gcTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCa                               <  1:3663462/150‑1 (MQ=255)

                                                                                                                               ccTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                       >  2:3368626/1‑150 (MQ=255)

                                                                                                                               ccTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                       >  1:1343321/1‑150 (MQ=255)

                                                                                                                               ccTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtc                       <  2:2018366/150‑1 (MQ=255)

                                                                                                                                cTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGtct                      >  1:752544/1‑150 (MQ=255)

                                                                                                                                  cACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTcc                    <  2:3701838/150‑1 (MQ=255)

                                                                                                                                    cccTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCAt                  <  1:336095/150‑1 (MQ=255)

                                                                                                                                     ccTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATg                 >  1:3731945/1‑150 (MQ=255)

                                                                                                                                        aaCCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGttt              <  2:177928/150‑1 (MQ=255)

                                                                                                                                           ccTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTca           >  1:1710396/1‑150 (MQ=255)

                                                                                                                                              ctcCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcacca        <  1:2973916/150‑1 (MQ=255)

                                                                                                                                               tcCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTcaccac       >  2:3256250/1‑150 (MQ=255)

                                                                                                                                                 cccAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  1:1611598/1‑150 (MQ=255)

                                                                                                                                                 cccAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  2:621758/1‑150 (MQ=255)

                                                                                                                                                 cccAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  1:609551/1‑150 (MQ=255)

                                                                                                                                                 cccAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  1:1176105/1‑150 (MQ=255)

                                                                                                                                                 cccAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  1:3493047/1‑150 (MQ=255)

                                                                                                                                                 cccAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaa     >  1:1080477/1‑150 (MQ=255)

                                                                                                                                                  ccAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaac    >  2:185436/1‑150 (MQ=255)

                                                                                                                                                   cAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCacaaca   <  1:1510157/150‑1 (MQ=255)

                                                                                                                                                    agagGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAg  >  2:614735/1‑150 (MQ=255)

                                                                                                                                                    agagGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAg  <  1:4040618/150‑1 (MQ=255)

TCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGG‑A‑CCGTATTGT‑GCAAATATTG‑TTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTGCCGAATATTCTCCAGCACCGTCTCCATGTTTTCACCACAACAG  >  NC_000913/1351080‑1351376

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