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	98.4%	110.0 / NA	62	noncoding (2543/2904 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/37);  minor base C (1/0);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AATCCTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCG  >  NC_000913/3946000‑3946478                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         aaTCCTGGTCGGTCCTCAGGAGGTTAGTGCAATGGCATAAGCCGGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTGCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCGGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGGCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGt                                                                                                                                                                                                                                       <  2:54027/250‑1 (MQ=11)  aTCCTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCATGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTAt                                                                                                                                                                                                                                     >  1:47997/1‑251 (MQ=14)    ccTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGACCAAGAGTTCATATCGACGGCGGGGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATgg                                                                                                                                                                                                                                   >  1:65548/1‑251 (MQ=14)          cGGACATGAGGAGGGTAGTGCAATGGCATAAGCCAGCTTGCCGGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTTCTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCTACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCGCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTc                                                                                                                                                                                                                             <  2:230391/251‑1 (MQ=11)           ggACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCg                                                                                                                                                                                                                            <  1:255766/251‑1 (MQ=14)              cATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGcc                                                                                                                                                                                                                          >  1:354130/1‑250 (MQ=14)                  aggaggTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTa                                                                                                                                                                                                                     <  2:313890/251‑1 (MQ=14)                          aGTGCAATGGCATAAGCCAGCTTGACGGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGGGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGGCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTa                                                                                                                                                                                                             <  2:219921/251‑1 (MQ=9)                                         gCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTGCTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTTGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTGCGCGATGGAAAGTGGGACGCGAGCTGGGTTTa                                                                                                                                                                                              <  2:379442/251‑1 (MQ=11)                                           cAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGCAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGTCGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                                                                                                                                                                               >  1:273134/1‑248 (MQ=14)                                             gCTTGACTGCGCGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCGTAGTGATCCGGTGGTTCTGAATGGAAGGGCCCTCGCTCAACGGCTAAAAGGTACTCCGGGGATAAAAGGCTGATACCGCCCAAGAGTTCATCTCGCCGGCGGAGTTTGGCACCCCGATGTCGGCTCATCACATACTGGAGCTGAACTAGGTCCCAAACGTAAGGCTGTTCCCCATTTAAAGTGGTACGCGCGCTGGGTTTAgcac                                                                                                                                                                                          >  1:198656/1‑248 (MQ=11)                                                tGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                                                                                                                                                                                       >  1:133031/1‑251 (MQ=14)                                                 gACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                                                                                                                                                                                      >  1:12205/1‑251 (MQ=14)                                                           gTGACGGCGCGAGCCGGTGCGCAAGCAGGACATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCCAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGGACGCGCGATGGGTTTAGAACGGAGTGAGACAGtt                                                                                                                                                                            >  2:381088/1‑251 (MQ=11)                                                             gACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGTTTCTGAATGGAAGGGCCATTGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCTCCTCGATGTCGGCTCATCACATCCTGGAGCTTAAGTAGGACCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCGGGGTTTAGAACGTCGTGAGACAGTTCg                                                                                                                                                                          <  1:248619/251‑1 (MQ=11)                                                              aCGGCGCGAGCAGGGGCGAAAGCAGGTCATAGTGATCCGGTGGTTCGGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTTGGTTTAGAACGTCGTGAGACAGTTCgg                                                                                                                                                                         <  2:331555/251‑1 (MQ=9)                                                                     gAGCAGGTGCGAAAGCAGGTCATAGTGATCAGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGCGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTAt                                                                                                                                                                  >  2:74682/1‑251 (MQ=14)                                                                       gCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCt                                                                                                                                                                <  2:334379/251‑1 (MQ=14)                                                                         aGGTGCGAAAGCAGGTCATCGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGGGACAGTTCGGTCCCTATCTGc                                                                                                                                                              >  2:268076/1‑251 (MQ=14)                                                                           gTGCGAAAGCAGGTCATAGGTATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTAATCCGTGGATAACAGGCTGATACCGCCAAAGAGTTCATATCGACGGCGGTGTTTGGCACCGCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCTAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCg                                                                                                                                                            <  2:282946/251‑1 (MQ=11)                                                                              cGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTgg                                                                                                                                                         <  2:133031/251‑1 (MQ=14)                                                                                     agggcATAGTGATCTGGTGGTTCGGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGTCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGTGCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGTTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTgg                                                                                                                                                  <  2:40372/247‑1 (MQ=11)                                                                                       gTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTgg                                                                                                                                                  >  1:160937/1‑249 (MQ=14)                                                                                            aGTGATCCGGTGGTTCTGAAGGGAGGGGACAGCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGAGGTGTTTGGCACCTCGATGTCGGCTCATCACGTCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAGGTGGTACGCGAGCTGGGTGTAGAGCGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTTGAGAACTg                                                                                                                                           <  2:104470/251‑1 (MQ=11)                                                                                             gTGATCCGGTGGGTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACATGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTAGATGTCGGCTCATCACATCCTGGAGCTGAAGTATGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCACTATCTGCCGTGGGCGCTGGAGAACTGa                                                                                                                                          <  2:316536/251‑1 (MQ=255)                                                                                                aTCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGCTGTCGGCTCATCCCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgt                                                                                                                                                                                        >  2:58779/1‑202 (MQ=9)                                                                                                aTCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgt                                                                                                                                                                                        <  1:58779/202‑1 (MQ=14)                                                                                                     gtggtTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGGACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGGGGGGGgct                                                                                                                                  >  2:401638/1‑251 (MQ=14)                                                                                                       ggtTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCGCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGTCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGGTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGgctg                                                                                                                                 <  1:379933/250‑1 (MQ=9)                                                                                                        gtTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGAGACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCGCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTTAGAGAGGTCGGTCCCTATCTGCCTTGGGCGCTGGAGAACGGAGGGTGgctgcg                                                                                                                               <  2:172340/251‑2 (MQ=11)                                                                                                                   aaGGGCCATAGCTCACCGGATACAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATAGACGGCGCTGTTTGGCACCTCGATGTCGGCTCATAACATCCTGGAGCTGAAGAAGATCCCAACAGAATGTCTGATCGGCATTTAAAGTGGTCCGCGAGCTGGGTTTGGCACGTCGTGAGTCCGTTCGGGCCCTATCTGCCGTGGGGGCTGGAGAACTGAGGGGGGCTGCGCCtcggacg                                                                                                                      >  2:192162/1‑243 (MQ=11)                                                                                                                    aGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACg                                                                                                                      >  1:326049/1‑248 (MQ=14)                                                                                                                       gCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCCTTTAAAGTGGTACGCGAGCTGGGTTTAGAAAGTCGTGAGACATTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGGGGGGGGCTGCTCCTAGTCCGAGAggg                                                                                                                >  2:146090/1‑250 (MQ=255)                                                                                                                            cGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTTAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGa                                                                                                           <  1:381088/251‑1 (MQ=14)                                                                                                                                cAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTg                                                                                                        <  2:124297/250‑1 (MQ=14)                                                                                                                                    ggATAAAAGGGACTCCGGGGATAACAGGCTGATACCGCACAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGGACTGAGGGGGGCTGCTCCTAGTACGCGAGGACCGGAGGGGACGc                                                                                                   >  2:262666/1‑251 (MQ=255)                                                                                                                                     gATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCAAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAAGGGGGCTGCTCCTAGTACGAGAGGACCGGAGAGGACGc                                                                                                   >  1:201263/1‑250 (MQ=9)                                                                                                                                           aGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTg                                                                                            >  1:8823/1‑251 (MQ=14)                                                                                                                                                  ccGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACAGCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGtt                                                                                       <  1:201554/249‑1 (MQ=14)                                                                                                                                                  ccGGGGATAACAGGCTGATACCGCCCAAGAGTTCATAGCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGtt                                                                                       <  2:222132/249‑1 (MQ=14)                                                                                                                                                    ggggATAACAGGCTGAGACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCggg                                                                                   <  1:223519/251‑1 (MQ=14)                                                                                                                                                          aaCAGGCTGATTCCGCCCAAGATTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTTATCACATCCTGGAGCTGAAGTAGGTCCCAAAGTGATGGCTGTTCTCCATTTAAAGTGGTACGGGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCATTCTGCCGTGGGCGCTGGAGAACTGCGGGGGGCTGCTCCTAGTGCGAGAGGACCGTAGTGGACGCATAACTGGTGTTCGTGTTGTc                                                                              <  2:375121/250‑1 (MQ=11)                                                                                                                                                                cTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGcc                                                                                                                                                             >  1:99854/1‑165 (MQ=17)                                                                                                                                                                cTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGcc                                                                                                                                                             <  2:99854/165‑1 (MQ=17)                                                                                                                                                                  taTACCGCCCAAGTGAACATAACGCCCGCGGTGTTTGGAACCTAGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCACAAATGTATGGCTGTTCGACCTTTAAAGTGGTACGCGAGCTGGGGTTAGAACGTCGTGAGACAGTTCGGACCCTATCTGCCGCGGGAGCTGGAGGACTGAGGGGGGCGGCACCAAGTACGAGAGGAGCGGAGAGGCCGCATCACTGGGGGTCGGGTTGTCAtagcaagg                                                                     >  2:99972/2‑244 (MQ=255)                                                                                                                                                                            aaGAGTTCATATCGGCGGCGGTGGTTGGCACCTTGATGTCTGCGCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTATAGTGGTACGCGAGCGGGGTTTAGAACGGCGTGAGACAGTTCGGTCCCTATCGGCCGTGGGCGCTGGAGAACTGAGGGGGTCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATTTCACTGCCCg                                                           <  2:65548/251‑1 (MQ=11)                                                                                                                                                                            aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCg                                                           <  1:130377/251‑1 (MQ=14)                                                                                                                                                                            aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCg                                                           <  1:257645/251‑1 (MQ=14)                                                                                                                                                                               agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTa                                                        <  1:213349/251‑1 (MQ=14)                                                                                                                                                                                    atatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCTGGGAGGTGAAGTAGGTCGCAAAGGTATGGCTGGTCGCCAGTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGGCCCTATCTGCCGGGGGCGCGGGAGAACTGAGGGGGGCTGCTCCTATTACGAGGGGGCCGGAGTGGACGTATCTATGGTGGTGGTGTTGTCATGCCAATGTCGCTGCCCGGTAGCTaa                                                   <  2:336890/251‑1 (MQ=11)                                                                                                                                                                                            ggcggTGTTTGGCACCTCGCTGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGGGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCCCTGCACGGTAGCTAAATGCGGa                                            >  2:49502/1‑250 (MQ=9)                                                                                                                                                                                               ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAagag                                        <  1:78315/251‑1 (MQ=14)                                                                                                                                                                                                gtgtTTGGCACCTCGATGGCGGCTCATCACATCCTGGATCTGAAGTAGGTCCCAAAGGTATGGCGGTTCGCCATTTAAAGTGGTACGCGAGCTGGGGTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCAGCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAgaga                                       <  2:273134/251‑1 (MQ=255)                                                                                                                                                                                                      ggccccTCGATGTCGGCTCGTCACATCCTGGAGCCGAAGTAGGTCCCAAAGGTATGGCTGTTCGGCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTTAGACAGTTCGGTCCCGATCTGGCGTGGGCGCTGGAGAACTGAGGGGGGCGGCGCCAAGTACGAGAGGACCGGAGTGGACGCATCACTGGAGGTCGGGTTGTCATGCCAATGGCGCTGACCGGTAGTTAAATTCGGAAGAGTTAAGt                                  <  2:187472/246‑1 (MQ=11)                                                                                                                                                                                                        cACCTCGATGTCGGATCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGCATGGCTGTTCGCCATTTACAGCGGTACCCGCGCTGGGTTTAGAACGTCGTGAGACAGTTCGATCCCTATCTGCCCTGGGCGCACGAGAACTGAGGGGGGCTGCTACTCGTACGAGAGGACCGGCGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCCCACCCCGGTAGCCAAATGCGGAAGAGATAAGTGCt                               >  2:183267/1‑251 (MQ=11)                                                                                                                                                                                                            tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGAGCTGaa                            >  1:158652/1‑250 (MQ=14)                                                                                                                                                                                                             cGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGaa                            <  2:115692/249‑1 (MQ=14)                                                                                                                                                                                                                   cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCt                     <  2:326049/250‑1 (MQ=14)                                                                                                                                                                                                                     gcgcagaccaTCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTGAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACGGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAg                  <  2:354130/243‑1 (MQ=9)                                                                                                                                                                                                                           cacaTCCTGGAGGTTAAGTAGGTCCCATAGGTATGGCTGTGCGACATTTAAAGTGGTGCGCGAGCTGGGCTGAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTTCGAGAGGATCGGAGTGGACGCATCACTGGTGTTCGAGTTGGCATGCCAATGGCACTGCCCGGTAGCTAAATGCAGAAGAGATAAGTGCTGAAAGCATCTAAGCACGa             <  2:276881/250‑1 (MQ=11)                                                                                                                                                                                                                                ccTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTTGTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCGAAGCACGAAACTTg       <  2:272902/251‑1 (MQ=9)                                                                                                                                                                                                                                     aGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCg  >  1:254090/1‑251 (MQ=11)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         AATCCTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAAGCCAGCTTGACTGCGAGCGTGACGGCGCGAGCAGGTGCGAAAGCAGGTCATAGTGATCCGGTGGTTCTGAATGGAAGGGCCATCGCTCAACGGATAAAAGGTACTCCGGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGGGCGCTGGAGAACTGAGGGGGGCTGCTCCTAGTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCG  >  NC_000913/3946000‑3946478

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