Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,946,246	G→A	noncoding (2543/2904 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,946,246	0	G	A	100.0%	124.3 / NA	69	noncoding (2543/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/34);  total (35/34)

GGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCG  >  NC_000913/3946006‑3946478                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         ggtctgccatcatGAGGTTAGTGCAAGGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTGCATATCGACGGCGGTGTTGGGCACCTCGATGGCGTCTCATCACAGCCTGGAGCTGAAGGAGGTCCCATAGGTATGGCTg                                                                                                                                                                                                                                <  2:46599/244‑1 (MQ=11)  gTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGCTAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCCAAGGTATGGCTGt                                                                                                                                                                                                                               >  1:186673/1‑251 (MQ=14)   tCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGtt                                                                                                                                                                                                                              >  1:324517/1‑251 (MQ=14)                   tAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGGAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTgg                                                                                                                                                                                                               <  1:277665/250‑1 (MQ=14)                   tAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTgg                                                                                                                                                                                                               <  1:242127/250‑1 (MQ=14)                               aTAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                                                                                                                                                                                  <  2:93489/251‑1 (MQ=14)                                 aaGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                                                                                                                                                                                 <  1:294642/250‑1 (MQ=14)                                   gCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTTCCCATTTAAAGTGGTACGCGAGCTggggttt                                                                                                                                                                                              >  2:240460/1‑247 (MQ=9)                                   gCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                                                                                                                                                                               >  1:124659/1‑250 (MQ=14)                                     cAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGa                                                                                                                                                                                            <  1:75179/251‑1 (MQ=14)                                      aGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATAGCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGaa                                                                                                                                                                                           <  2:248117/251‑1 (MQ=14)                                                      tGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCAAAGAGTTCATATCGAAGGCGGAGATTGGCACCTCGATGTCGGCTCCTCACATCCTGGAGCTGAAGTAGGTCACACAGGTATGGCTGCTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                                                                                                                                                                           >  1:190087/1‑251 (MQ=11)                                                        aCGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGAGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCCTTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGGGACAGTTCgg                                                                                                                                                                         >  2:168865/1‑251 (MQ=9)                                                             gcgAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCt                                                                                                                                                                    <  1:32723/251‑1 (MQ=14)                                                              cgAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGAACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGATGGGTTTAGAACGTCGGGGGACAGTTCGGTCCCt                                                                                                                                                                    >  2:312884/1‑250 (MQ=9)                                                                 gCAGGTGCGAAAGCATGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGGCGTGAGACAGTTCGGTCCCTATCt                                                                                                                                                                <  2:231164/251‑1 (MQ=14)                                                                  cgggTGCGAAAGCAGGTCATATTGTCCCGGTGGTTCGGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACGCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCAGAGCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTTAAGTAGGTCCCAAAGTTATTGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGGTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTg                                                                                                                                                               <  2:190163/249‑1 (MQ=11)                                                                              cAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTg                                                                                                                                                   <  1:326237/251‑1 (MQ=14)                                                                                 gTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAATGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGgag                                                                                                                                                <  1:260397/251‑1 (MQ=14)                                                                                      aGTGATCCGGTGGTTCTGAATGGACGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCCTTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTg                                                                                                                                           >  2:376380/1‑251 (MQ=14)                                                                                      aGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTg                                                                                                                                           <  2:365926/251‑1 (MQ=14)                                                                                            ccGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAggggg                                                                                                                                     >  2:67611/1‑251 (MQ=14)                                                                                              ggtggttctgcctgacacggCCATCGCTCACCGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCAAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGGGAACTGAGGGGGgc                                                                                                                                   >  2:59854/19‑251 (MQ=11)                                                                                               gtggtTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgct                                                                                                                                  >  1:244038/1‑251 (MQ=14)                                                                                                tggtTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCCTTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgt                                                                                                                                                                                        >  1:400823/1‑196 (MQ=17)                                                                                                tggtTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCGCGATGTCGGCTCATCACATCCAGGAGCTGAAGTAGGTCCCAAAGGGAGGGCTGTTCGCCATTTAAAGTGCTACGCCAGCTGGGTTTAGAAcgt                                                                                                                                                                                        <  2:400823/196‑1 (MQ=11)                                                                                                     cTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCGGGAGAACTGAGGGGGGCTGCTccc                                                                                                                            >  2:379662/1‑250 (MQ=14)                                                                                                       gAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCt                                                                                                                            <  2:331660/249‑1 (MQ=14)                                                                                                              aGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACgagg                                                                                                                   >  1:272600/1‑250 (MQ=14)                                                                                                                    aTCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCAACGGCGGTGTTTGGCACCTCGATGTAAGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTAAGGATGTTCGCCATTTCAAGGGGTACGCGCGCTGGGTTTCGAACGTAGGGAGCCAGTGAGGTCCCTATCTGCCGCGGGCGCTGGGGAACTGAGGGGGGCTGCTCCTAGGACGAGAGGACCg                                                                                                             >  2:191321/1‑251 (MQ=11)                                                                                                                        cTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAgg                                                                                                         >  1:155455/1‑250 (MQ=14)                                                                                                                        cTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAgg                                                                                                         >  1:370676/1‑250 (MQ=14)                                                                                                                          cAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGGGGTAAGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGGGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTTGTACGAGAGGACCGGAGTgg                                                                                                       >  2:368246/1‑251 (MQ=9)                                                                                                                            aCGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTCTCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGAc                                                                                                     >  1:434835/1‑251 (MQ=14)                                                                                                                            aCGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGAc                                                                                                     <  1:158970/251‑1 (MQ=14)                                                                                                                            aCGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGAc                                                                                                     >  1:132162/1‑251 (MQ=14)                                                                                                                                      ggTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATTGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGGTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTgg                                                                                           <  1:59854/251‑1 (MQ=14)                                                                                                                                           tCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGGGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTc                                                                                      >  1:983/1‑251 (MQ=14)                                                                                                                                            ccGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCg                                                                                     <  2:15606/251‑1 (MQ=14)                                                                                                                                             cGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCgg                                                                                    >  1:169431/1‑251 (MQ=14)                                                                                                                                              ggggATAACAGGCTGATCCCGCCCAAGCGTTCATATCGACGGCTGTGGTTGGCACCGCGATGTCGGCTCATCCCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTTCGCTAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTg                                                                                                                                                          <  1:229218/180‑1 (MQ=11)                                                                                                                                                   tAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGGGGGGGGCTGCTCCTAGTACGGGAGGACCGGGGTGGACGCATCACTGGTGTTCGGGTTGTc                                                                              >  2:53017/1‑251 (MQ=9)                                                                                                                                                     aCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCAt                                                                            <  2:132162/251‑1 (MQ=14)                                                                                                                                                                  gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTg                                                               <  1:247208/251‑1 (MQ=14)                                                                                                                                                                    ccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGAACTGcc                                                             >  1:66122/1‑251 (MQ=14)                                                                                                                                                                     cAAGAGTTCATATCGACGGCGGTGTTTGGCATCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTGAGAACGTCGGGAGACAGTGCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGccc                                                            <  2:91323/251‑1 (MQ=9)                                                                                                                                                                     cAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGccc                                                            >  2:73433/1‑251 (MQ=14)                                                                                                                                                                      aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCg                                                           <  2:983/251‑1 (MQ=14)                                                                                                                                                                          gtgcagatCGCCGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCAATTTAAAGTGGTACGAGAGCTGGGTTTAGAACGTCGTGAGACAGTTCTGTCCCTATCTGCCGTGTGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAg                                                       <  2:3552/245‑1 (MQ=11)                                                                                                                                                                            tCATATCGACGGCGGTGTTTGGCTCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGc                                                      <  1:371465/250‑1 (MQ=14)                                                                                                                                                                             cATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCGCATCATATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGTCTGTTCGCCATTTAAATTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGGTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTATTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTa                                                    <  2:56502/251‑1 (MQ=11)                                                                                                                                                                             cATAGCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACGGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTGCGGGTTGTCATGCCAATGTCACTGCCCGGTAGCTa                                                    <  2:31507/251‑1 (MQ=9)                                                                                                                                                                              atatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCCAAGGTATGGCTGTTCGCCATATAACGTGGTACGCGAGCTGGGGTTAGAACGTCGTGAGACCATTCGGACCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGCGGACCGGAGTGGACGCATCACTGGAGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGcga                                                    >  2:183605/1‑248 (MQ=11)                                                                                                                                                                                atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGaa                                                                                                                                                                                           <  1:241147/113‑1 (MQ=18)                                                                                                                                                                                atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGaa                                                                                                                                                                                           >  2:241147/1‑113 (MQ=18)                                                                                                                                                                                      ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGc                                                                   >  1:278118/1‑227 (MQ=14)                                                                                                                                                                                      ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGc                                                                   <  2:278118/227‑1 (MQ=14)                                                                                                                                                                                           tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTCGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCCCTGCCCGGTAGCTAAATGCGGAag                                          >  2:279960/1‑247 (MQ=14)                                                                                                                                                                                             tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa                                                                                                           <  1:75400/180‑1 (MQ=17)                                                                                                                                                                                             tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa                                                                                                           >  2:75400/1‑180 (MQ=17)                                                                                                                                                                                             tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATa                                     >  1:133452/1‑250 (MQ=14)                                                                                                                                                                                               tGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGGGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCGATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGGGGAGGCATCACTGGTGTTGGGGTTGTCATGCCAATGGCACGGCCCGGTAGCTAAATGCGGAAGAGATAAg                                   <  2:68319/250‑1 (MQ=11)                                                                                                                                                                                                  cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCt                               >  2:214625/1‑251 (MQ=14)                                                                                                                                                                                                     cTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGa                             >  1:236877/1‑250 (MQ=14)                                                                                                                                                                                                       cGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGaa                            >  1:237727/1‑249 (MQ=14)                                                                                                                                                                                                                  catcaCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCGCTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGc                 <  1:168865/249‑1 (MQ=14)                                                                                                                                                                                                                  catcaCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGGTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGCAGAGATAAGTGCTGAAAGCATCGAAGCAc               >  2:135662/1‑251 (MQ=9)                                                                                                                                                                                                                            gggAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGcc     <  2:124659/250‑1 (MQ=14)                                                                                                                                                                                                                               aGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCg  <  2:133452/251‑1 (MQ=11)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         GGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCG  >  NC_000913/3946006‑3946478

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