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	98.0%	112.6 / NA	50	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/26);  minor base A (0/1);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTG  >  NC_000913/2730701‑2731115                                                                                                                                                                          |                                                                                                                                                                                                                                                       cTCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCCTCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTTGCTGGTCATCCTCTCATACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACCAGCTAATCCCATCTGGACACATCCGATGGCAAGAGGCCCGaa                                                                                                                                                                       >  2:105407/1‑250 (MQ=255)   cccTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg                                                                                                                                                                      <  1:399068/249‑1 (MQ=11)    ccTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGTGCACCTCCGATGGCAAGAGTCCCGAAggcc                                                                                                                                                                   >  1:388533/1‑249 (MQ=9)     cTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGCTGGCAAGAGGCCCGAAGGTcc                                                                                                                                                                  >  1:509662/1‑251 (MQ=14)       tcctccCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc                                                                                                                                                                  <  1:87041/249‑1 (MQ=14)           ccccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCGACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct                                                                                                                                                            >  1:155561/1‑251 (MQ=14)               gCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGAGGCATGGCTGAATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGCCCCCCTCTTTgg                                                                                                                                                        >  1:310325/1‑251 (MQ=9)                  gaaggtacgttaccccccGAAGGCCTTCTTCATACACCCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATCCCCCACTTCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCATTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                                                                                                                        <  2:425096/233‑1 (MQ=11)                    aaGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTg                                                                                                                                                   >  1:91518/1‑251 (MQ=14)                      gTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                                                                                                                 <  1:365980/251‑1 (MQ=14)                           ttACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAtt                                                                                                                                            >  1:499090/1‑251 (MQ=17)                           ttACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAtt                                                                                                                                            >  1:467485/1‑251 (MQ=17)                             cccatccGAAGGCCTTCTTCATTCACTCGGCATGGCTGCATCAGGCTGGCGCCCATTTTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGGCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGACCGTCGCCTTAGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGTGCACCTCCGATGGCAAGGGGCCCGAATGTCCCCCTCTTTGGGCtttcgtcat                                                                                                                                             <  2:505524/243‑8 (MQ=11)                               aaCCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                                                                                                                        >  1:214766/1‑251 (MQ=18)                                    gAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                                                                                                                   >  1:281093/1‑251 (MQ=18)                                                   acacGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTcc                                                                                                                     >  1:306275/1‑250 (MQ=18)                                                     acGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagt                                                                                                                  >  2:628170/1‑251 (MQ=18)                                                      cgcgGCCTTGCTGCATCAGGCTTGCGCCCATTGTGCAAGATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGTGATCGTCTCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGCCAGTATGCGGTATTAGCTACCGGTTCCagtc                                                                                                                 <  2:45346/251‑2 (MQ=11)                                                           cATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt                                                                                                            <  2:276171/251‑1 (MQ=18)                                                            aTGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATc                                                                                                           <  1:5896/251‑1 (MQ=18)                                                             tGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTGCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc                                                                                                          >  2:599967/1‑251 (MQ=18)                                                              ggCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCCCCCGTAGGAGTCGGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTCATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATccc                                                                                                         <  2:430897/251‑1 (MQ=16)                                                               gCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc                                                                                                        >  1:120929/1‑251 (MQ=18)                                                                      cAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATc                                                                                                 <  2:332768/251‑1 (MQ=18)                                                                         gCTTGCGCCCATTGTGCAATATCCCCCACTGCTGCATCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGGCTTGCGACATTTTGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAgg                                                                                              >  2:530435/1‑251 (MQ=255)                                                                            tGCGCCCATTGTGCAATATTCCCCACGGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑ATCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGTCCCGAAGGTCCCCCTCTGTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTGCCAGTAGTTATCCCCCTCCATCAGGCa                                                                                            <  1:483396/250‑1 (MQ=255)                                                                             gcgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATTCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAg                                                                                           <  1:168576/250‑1 (MQ=18)                                                                             gcgcCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGc                                                                                          <  1:74868/251‑1 (MQ=18)                                                                                cccATTGTGCGATCTTCCCCAATGCTGCCTCCCGTAGGAGTCTAGACAGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCATGTCAGCCAGGGATCGTCGCCTTGGTG‑GGCCGTTACCCCACCAACAAGCTAATCTCATCTGGGCACAACCGAGGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATGATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAGCAGGCAGCTTc                                                                                       <  2:281495/251‑1 (MQ=11)                                                                                cccATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCTGGCAGCTTc                                                                                       <  1:620502/251‑1 (MQ=21)                                                                                 ccATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTcc                                                                                      >  2:580769/1‑251 (MQ=21)                                                                                  cATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGCCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGTGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTccc                                                                                     <  2:614533/251‑1 (MQ=21)                                                                                  cATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGTCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTccc                                                                                     <  1:319190/251‑1 (MQ=21)                                                                                   aTTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTccc                                                                                     >  1:32050/1‑250 (MQ=25)                                                                                        gCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCAGCAGGCAGCTTCCCAGACa                                                                                <  2:514752/250‑1 (MQ=25)                                                                                          atgaattccctcttctgcCTCCCGTAGGAGTCGGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACAAGCCAGGGATCGTCGCCTTGGTG‑CGCCTCTACCCCACCAACAAGCTAATCCCATCTGTGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCACTCCATCAGGCAGCTTCTCAGACATTa                                                                             <  2:57688/237‑1 (MQ=11)                                                                                            tatTCCCCACTGCTGCCTCACGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGATACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACAttgct                                                                           >  2:499756/1‑248 (MQ=255)                                                                                              ttCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTc                                                                          >  1:454585/1‑250 (MQ=25)                                                                                               tCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTccc                                                                        >  2:26947/1‑249 (MQ=21)                                                                                                          gcCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGcc                                                               <  1:219790/249‑1 (MQ=25)                                                                                                               ccGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCGGGTCAGCCTCTCAGACCAGCCAGGGATCGTCGCGTGGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTc                                                        <  1:329516/251‑1 (MQ=255)                                                                                                                gtgaggagtcgagtccGTGTCGCAGTTCCAGGGTGGCTGGTCATCCCCTCAGACCAGCAAGGGATCGTCGC‑TTGGGGTAGCCGTTCCCCCGCCAACAAGCTAATCCCATCTGGCCACTTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCCCCATCAGGCAGCTGCCCAGACATTACCCGACCGTCCGCCACTCGTCa                                                       <  2:22279/237‑1 (MQ=11)                                                                                                                   aGGAGTCTGTACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGc                                                     <  1:245148/250‑1 (MQ=25)                                                                                                                      aGTCTGGACCGTGTCTCAGTTCCATGGTGGCTGGCAATCCTCTCACCCCAGCCAGGGATCGGTGCCTTAGTG‑AGCCGTTCCCCCACCAACAACCTAATCCCATTCGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTCTTAGCTCCCGTTTCCAGTAGTGATCCCCCCCCCTCAGGCAGCTGCCCAGGCATTACTCACCCTTCCGCCACTCGTCAGCGAag                                                 <  2:253946/251‑1 (MQ=14)                                                                                                                            gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAgcagca                                            >  1:531153/1‑250 (MQ=37)                                                                                                                                      cAGTTCCAGTGTGGCTGGTCAGCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTcc                                  <  1:397575/250‑1 (MQ=255)                                                                                                                                                gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTGTCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTc                       <  1:271499/251‑1 (MQ=255)                                                                                                                                                gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGACGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTCGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCCAGATGCTTCTTGTGACCGtcg                       >  2:543019/1‑249 (MQ=14)                                                                                                                                                 tgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGTTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCg                      >  1:160192/1‑251 (MQ=255)                                                                                                                                                               cccAGACCAGCCAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAACCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTGGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCGAGCTGCTTCCTGTTACCGTTCGACGTGCATgtgt          <  2:499090/247‑1 (MQ=255)                                                                                                                                                                     ccagccagGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTg  <  1:87175/251‑1 (MQ=38)                                                                                                                                                                          |                                                                                                                                                                                                                                                       CTCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTG‑AGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCTG  >  NC_000913/2730701‑2731115

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