Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,730,870	T→C	noncoding (288/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,730,870	0	T	C	100.0%	158.5 / NA	68	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/33);  total (35/33)

TCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCT  >  NC_000913/2730702‑2731114                                                                                                                                                                          |                                                                                                                                                                                                                                                       tCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑Agg                                                                                                                                                                     >  1:139515/1‑251 (MQ=14)      tcctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACCAGCTAATCCCATCTGGGCACCTCCGCTGGCAAGAGGCCCGA‑AGGTcccc                                                                                                                                                                >  2:431500/1‑251 (MQ=9)      tcctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTcc                                                                                                                                                                  <  1:526102/249‑1 (MQ=14)      tcctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGCcccc                                                                                                                                                                >  2:409825/1‑251 (MQ=14)          ccccgctgaaagtctggttcAACCCGAAGGCCCTCTTCATACACCCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATCTTCCCCACTGCTGCCTCCCGTAGGAGTCGGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTGCCCCACCAACAAGCTAACCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGGCCCCct                                                                                                                                                              <  2:287467/230‑1 (MQ=11)           cccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCtt                                                                                                                                                           >  1:539509/1‑251 (MQ=14)            ccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCttt                                                                                                                                                          >  2:731669/1‑251 (MQ=14)                agAAAGTTCATTACACCACGAAGGCCTGCTTTCTACACGCGGCATGGCGTCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCTTAGGAGTCTGTCCCGTGTCTCTGTTCCAGTTTGGCTGTTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCGTCTGGGCACATCCGATGGCAAGAGGCCCTA‑AGGTCCCCCTCTTTGG‑tg                                                                                                                                                     <  2:548737/250‑2 (MQ=11)                 gAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGTCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTGCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑t                                                                                                                                                      <  1:40019/249‑1 (MQ=14)                   aaGTACTATACAACCCGAAGGCCTTCTACATGCTCGCGTCAGGGCTGCATCAGGCTTGCGCCAATTGTGCAATACTCCCCACTGCTGCCTCCCGGAGGAGTCTGGACCGTGTCTCAGTTCGAGTGTGGCTGGTC‑ATGCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCGCATCCGAGGGCAAGAGTCCCGA‑AGGTCGCCCTCTTTGG‑TCTTg                                                                                                                                                  <  2:66256/251‑1 (MQ=11)                      tACTTTACAACCCGAAGGCCTTCTTCATCCACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCGGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCg                                                                                                                                                <  2:583959/250‑1 (MQ=14)                         tttACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCCGACCAGCCAGGGATCGTCGCCTTTGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGCTGGCAAGAGGCCCGACAGG‑CCCCCTCTTTGG‑TCTTGCGAccc                                                                                                                                            >  2:726901/1‑249 (MQ=11)                          ttACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACAtt                                                                                                                                           >  1:529966/1‑251 (MQ=17)                                 ccGAAGGCCTTCTTCATACCCGCGGCATGGCTGCATCAGGCTTTCGCCCATTGTGCAATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGt                                                                                                                                    <  2:674275/251‑1 (MQ=255)                                 ccGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGt                                                                                                                                    <  1:731669/251‑1 (MQ=18)                                   gAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTAt                                                                                                                                  >  1:182961/1‑251 (MQ=18)                                   gAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTAt                                                                                                                                  <  1:521271/251‑1 (MQ=18)                                       gCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGc                                                                                                                              <  2:121493/251‑1 (MQ=18)                                         cttcttCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCt                                                                                                                             >  2:257038/1‑250 (MQ=18)                                           tcttcATACCCGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGCGTCTGGAACGTGTCTCAGTGCCAGTGTGGCTGGTC‑ATCCTCTCAGGCCAGCCCGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGCCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTAc                                                                                                                           <  2:758531/250‑1 (MQ=11)                                            cttcATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTTAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCg                                                                                                                         >  2:44753/1‑251 (MQ=18)                                            cttcATACACGCGGCATGGCTGCATCACGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCCTCTGGGCACATCCGATGGCAAGAGGCCCGC‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTCTGCGGTATTAGCTACCg                                                                                                                         >  2:413180/1‑251 (MQ=255)                                                 tACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTcc                                                                                                                    >  1:402266/1‑251 (MQ=18)                                                  acacGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTcc                                                                                                                    >  1:430757/1‑250 (MQ=18)                                                   cacGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAAGCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                                                                                                                  <  1:673531/251‑1 (MQ=18)                                                    ccgcgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt                                                                                                                 <  1:549490/250‑1 (MQ=18)                                                     cgcgGCATGGCTGCTTCCGGCTTGCGCCCATTGTGCATTATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGTGC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACAGCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagta                                                                                                                <  2:179702/251‑1 (MQ=11)                                                     cgcgGCATGGCTGCAACAGGCTTGCGCCCATTGTGCAACATTCCCCACTGCTGCCTCCCGGAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGACAGGGATCGTCACCTTGGTGCGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACCTCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGGATCT‑GCGGCATTATGCGGTATTATCTACCGTTTcaagg                                                                                                                 >  2:705123/1‑246 (MQ=11)                                                      gcgGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag                                                                                                               <  1:198488/251‑1 (MQ=18)                                                         gCATGGCGGCATCAGGGTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTa                                                                                                            <  2:307733/251‑1 (MQ=18)                                                          cATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt                                                                                                           >  1:689017/1‑251 (MQ=18)                                                              gtcgcccgtcctgtatgcCCTTTGGGCAATAGTCCCCCCTGCGGCCTCACGTAGGAGTCTGGACCGTGTCCCAGTTCCAGTGGGGCGGGTC‑ATCCGCTCAGGCCAGCCAGGGATCGTCGCCTTGGGGAGCCGTTACCCCACCAACAATCTAATCCCATCTGGGCACATCCGATGGCAAGGTGACCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTTGTTATcccc                                                                                                       <  2:245842/235‑1 (MQ=11)                                                                    gcAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAt                                                                                                 <  2:182961/250‑1 (MQ=18)                                                                     cAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGAATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc                                                                                                >  2:470190/1‑251 (MQ=18)                                                                         cTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCctc                                                                                                                                                             <  2:439999/187‑1 (MQ=17)                                                                         cTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCctc                                                                                                                                                             >  1:439999/1‑187 (MQ=17)                                                                            gcgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGc                                                                                         >  2:439290/1‑251 (MQ=18)                                                                             cgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCt                                                                                        >  1:142072/1‑251 (MQ=18)                                                                              gcCCATTGTGCCATATTCCCCACTGCTGCCTCCCGTAGGCGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGGCATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc                                                                                                       >  1:245842/1‑235 (MQ=255)                                                                                ccATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTcc                                                                                     <  1:60974/251‑1 (MQ=21)                                                                                ccATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGACCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTcc                                                                                     >  1:246754/1‑251 (MQ=21)                                                                                             ttCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCCGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGCCATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTcc                                                                        >  1:4552/1‑250 (MQ=19)                                                                                               ccccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTccac                                                                      >  2:215398/1‑248 (MQ=21)                                                                                                       ctgcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGcc                                                              >  2:740943/1‑251 (MQ=21)                                                                                                         gcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCa                                                                                   <  1:476677/228‑1 (MQ=25)                                                                                                         gcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTACCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCa                                                                                   >  2:476677/1‑228 (MQ=25)                                                                                                             cccGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGt                                                        <  1:249273/251‑1 (MQ=21)                                                                                                             cccGTAGGAGCCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGCGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGCTGGCAAGAGGCCCGC‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTGGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTNCCCAGACCTTACTCACCCCTCCGCCACtagc                                                        >  2:362796/1‑248 (MQ=11)                                                                                                              ccGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGGCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCACGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTc                                                       >  1:502958/1‑251 (MQ=21)                                                                                                              ccGTAGGAGTATGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATACTCTAAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTc                                                       >  1:614988/1‑251 (MQ=255)                                                                                                               cGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCCAGAGGCCCGA‑AGGGCCCCCGCTTTGG‑TCTTGCGACATTATGCGGGATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTc                                                       >  2:167462/1‑250 (MQ=255)                                                                                                                  agtggtctggccctggcctcAGTTCCAGGGTGGC‑GGTCCACCCCCCCAGCCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCACCCCACCAAGCGAACCCCATCGGGGCACATCCGATGGCAAGAGGCCCGA‑AGGCCCCCCTCTTTGG‑TTTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCAAGACATTACTCACCCGTCCGCGACTCGTCAGc                                                    <  2:430757/233‑1 (MQ=11)                                                                                                                           gACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAgcagca                                           >  1:328129/1‑250 (MQ=37)                                                                                                                                      aGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCNCCTCGTCAGCGAAGCAGCACGCTGCTTCCtt                               >  2:395762/1‑250 (MQ=255)                                                                                                                                       gtgccAGGGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAg                                                                                  <  2:38635/196‑1 (MQ=255)                                                                                                                                       gTTCCTGTGTGGCTGCTC‑ACCCCCTCACACCAGCCAGGGATCGTCGCCTTGGTGCGCCGTTACCCCACCAACAATCAAATCCCATCTGGGCACACCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGTCATTATGCGGTATTAGCTACCGTTTCCAGTAGTTCTCCCCCTCCATCAGGCAGCTTCCCAGACATTCCCCTCCCGTCCGCCACTCGACAGCGAAGCAGCAAGCTGCTTCCTGt                              <  2:269332/251‑1 (MQ=14)                                                                                                                                       gTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAg                                                                                  >  1:38635/1‑199 (MQ=37)                                                                                                                                             ttgtgGCTGTTC‑ATCCTCTCAGTCAGGCCGGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCACCTACCTAATCCCTTCTGGGCCCACCCGATGGCAAGAGGCCGGA‑AGCTCCCCCGCTTTGG‑CCGTGCGACATTATGCGGGTTTAGCCTCCTTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCCCCCGCCCGCCACTCGTCAGCGAATCATCAAGCTGCTTCCTGTTACCg                         <  2:518288/249‑1 (MQ=2)                                                                                                                                                 ggCTGGTC‑ACCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGa                    <  1:441657/251‑1 (MQ=255)                                                                                                                                                     gggccacccTCTCAGACCAGCCAGGGACCGTCGCCTTGGTGAGCCTTTCCCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGGGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTCCTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACtt                 <  2:139515/244‑1 (MQ=14)                                                                                                                                                     ggTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTg                <  1:470190/251‑1 (MQ=38)                                                                                                                                                      gTC‑ATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGc               >  2:80905/1‑251 (MQ=38)                                                                                                                                                           cccTCCCAGACCAGCCAGGGTACGTCTCCTTGGTGAGCCGTACCCCCTCCAACAAGCTAACCCCATCTGGGCCCATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACCCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtgt           <  2:569924/250‑1 (MQ=14)                                                                                                                                                                ccAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTAAGCGAAGCAGCAAGCTGCTACCTGTTACCGTTCGACTTGCATGTGTTAg      <  1:431500/250‑1 (MQ=255)                                                                                                                                                                 cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAg      <  1:395762/250‑1 (MQ=39)                                                                                                                                                                  aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGc    <  1:409825/251‑1 (MQ=38)                                                                                                                                                                   gACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGcc   >  1:668296/1‑251 (MQ=38)                                                                                                                                                                    aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGACGCAGCAAGCTGACTCCTGTTACCGTTCGACTTGCATGTGttctgccg  >  2:4841/1‑245 (MQ=11)                                                                                                                                                                          |                                                                                                                                                                                                                                                       TCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTC‑ATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGA‑AGGTCCCCCTCTTTGG‑TCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCT  >  NC_000913/2730702‑2731114

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