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

ACTCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCT  >  NC_000913/2730700‑2731114                                                                                                                                                                           |                                                                                                                                                                                                                                                       aCTCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                                                                                                                       <  1:241406/251‑1 (MQ=11)    cccTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTGCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                                                                                                                                    <  1:62762/251‑1 (MQ=14)      cTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTcc                                                                                                                                                                  >  1:218671/1‑251 (MQ=14)            ccccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCCTCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCCCATCCGATGGCAAGAGGCCCGAAGGGCCCCctcc                                                                                                                                                            >  2:558564/1‑250 (MQ=9)            ccccGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                                                                                                                              <  1:480995/249‑1 (MQ=14)                   gAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGt                                                                                                                                                       <  1:177770/249‑1 (MQ=14)                   gAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCt                                                                                                                                                     <  1:32768/251‑1 (MQ=14)                    aaaGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCtt                                                                                                                                                    <  1:582603/251‑1 (MQ=14)                     aaGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTg                                                                                                                                                   <  1:4555/251‑1 (MQ=14)                       gTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCg                                                                                                                                                 <  1:288848/251‑1 (MQ=14)                         aCTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                                                                                                                                <  1:525153/250‑1 (MQ=14)                            ttACAGCCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                                                                                                                             <  2:12971/250‑1 (MQ=17)                                aaCCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GCGAGCCGTTACCCCACCACCAAGCTAATCCCATCTGGGCACATCAGATGGCAAGAGGCCCGANGGTCCCCCTCATTGGCCTTGCGACATTagcc                                                                                                                                        >  2:332706/1‑248 (MQ=11)                                  cccGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCCGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCAAGGGCTCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                                                                                                                      >  2:312157/1‑251 (MQ=18)                                  cccGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGACCCCCTCTTTGGTCTTGCGACATTATGCgg                                                                                                                                      >  2:132765/1‑251 (MQ=18)                                    cgtatgCCTTCTTCATCTACTCGGCATGTCTTGACCAGGCTTGCGCCCTGTGTGCAATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCGCAGACCAGCCAGGGATCGTCGCCTTG‑GTGGGCCGTTACCCCACCAACAAGCGAAGCCCATCTGGGCACATCCGATGGCAAGATGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATTCGGTa                                                                                                                                    <  2:654252/246‑1 (MQ=11)                                     gAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCCTCGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGGGTCTCAGTTCCAGGGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCCCCTTG‑GTGAGCCGTTACCCCACCAACAAGCCAATCCCATCTGGGCACATCCGACGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                                                                                                                   >  2:206746/1‑251 (MQ=11)                                     gAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                                                                                                                   >  1:156127/1‑251 (MQ=18)                                      aaGGCCTTCTTCATACACGCGGCATGTCTGCATCCGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGTCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCGTTGGTCTTGCGACATTATGCGGTAtt                                                                                                                                  <  2:442647/251‑1 (MQ=255)                                       aGGCCTTCTTCATACACGCGGCAGGGCTGCATCAGGCTTGCGCCCATGGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCGGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                                                                                                                                 <  2:186317/251‑1 (MQ=16)                                        ggCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGCCCCCCTCTTTGGTCTTGCGACATTATGCGGTAtt                                                                                                                                  >  1:562799/1‑249 (MQ=18)                                        ggCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GGGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGCTGGCCAGAGGCCCGAAGGTCCCCCTCTTTGGCCTTGCAACATTATGCGGCAtt                                                                                                                                  >  2:394539/1‑249 (MQ=11)                                         gCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                                                                                                                               <  1:725165/251‑1 (MQ=18)                                          ccTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGGGCAATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCGGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAACCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCt                                                                                                                              <  2:57642/251‑1 (MQ=255)                                           cgtcttcATACACGCGGCATGGCTGCATCAGNCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GAGAGCCGTTACCCCACCAACAAGCTAAGCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                                                                                                                             <  2:354385/249‑1 (MQ=16)                                             tcttctaacacGCGGCTTGGCTGCATCAGGCTTGCGCTCATTGTGCAATATGCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCGGTGTTGCTGGTCATCCTCTCAGACAAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAGGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATCAGCtcc                                                                                                                            <  2:709463/243‑3 (MQ=11)                                                tcATACACGCGGAATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGt                                                                                                                         <  2:133843/250‑1 (MQ=18)                                                  aTACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCCCATCCGCTGGCAAGAGGCCCCAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                                                                                                                      >  2:459219/1‑251 (MQ=16)                                                   tACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGGCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTcc                                                                                                                     >  2:94484/1‑251 (MQ=18)                                                          ggCATGGCTGCATCAGGCTTGCGCCCCTTGTGCAATATTCCCCCCTGCTGCCACCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGCTCGACCCCTTG‑GGGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGATGCCCGAAGGCCCCCCTCTTTGGTCTTGCGACCTTCTGCGGGATTAGCAACCGTTTCCAGTA‑Gtt                                                                                                             >  2:281521/1‑251 (MQ=11)                                                           gCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑Gtt                                                                                                             >  1:290652/1‑250 (MQ=18)                                                               ggCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCCACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                                                                                                                                   >  2:554243/1‑193 (MQ=17)                                                               ggCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                                                                                                                                   <  1:554243/193‑1 (MQ=17)                                                                   gCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATccccc                                                                                                      <  1:217557/249‑1 (MQ=18)                                                                      acAGGCTGGCGCCCATTGTGCAATATTCCCCACTGCGGCCTCCCGTAGGATTCTGGACCGGGTCTCAGTTCCAGGGTGGCTGGTCATCCTCTCAGACCAGCAAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGTGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCa                                                                                                  <  2:133183/249‑1 (MQ=11)                                                                       cAGGCTTGCGCCCATTGTTCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATc                                                                                                >  2:600716/1‑251 (MQ=18)                                                                       cAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATc                                                                                                <  2:224576/251‑1 (MQ=18)                                                                          gCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAgg                                                                                             <  1:538596/251‑1 (MQ=18)                                                                           cTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGc                                                                                            >  1:276324/1‑251 (MQ=18)                                                                             ttcgcCCATGGTGCAAGATTCCCCTCTGCTGCCTCCCGTAGGAGTCGGGACGGTGTCTCGGTGCCAGGGTGGCTGGTCCTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGTGTG‑GCCGTTACCCCCCCAACAAGCTAAGCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCa                                                                                           <  2:164380/248‑1 (MQ=11)                                                                                ccccATTGTGCAATATTCCCCCCTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTGATCCCCCTCCATCAGGCAGCtt                                                                                       <  2:510883/250‑1 (MQ=19)                                                                                cccaaTTGTGCAATATGCCCCACTGCTGCCTCCCGTAGGAGTCTGGACAGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGACGATACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGTAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATGAGATACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCtt                                                                                       <  2:524987/247‑1 (MQ=11)                                                                                   cATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTccc                                                                                    >  2:11752/1‑251 (MQ=21)                                                                                          cAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACAtt                                                                             >  1:135788/1‑251 (MQ=21)                                                                                           aaTATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑TTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCTCATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTa                                                                            <  1:91319/251‑1 (MQ=21)                                                                                            ctatTCCCCACTGCTGCCTCACGTAGGAGTCTGGACCGTGTCTCAGTTCCAGAGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTAGCCTTG‑GAGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGAACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCATTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCA‑TATGATATCCCCCTCCATCGGGCAGCTTCCCAGACAttcc                                                                           >  2:569647/2‑249 (MQ=14)                                                                                                 ccccACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCCTCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCTTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGCAATTACTCAcc                                                                      >  1:399264/1‑251 (MQ=255)                                                                                                 ccccACTGCTGCCTCCCGTAGGAGGCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCAcc                                                                      <  1:477495/251‑1 (MQ=21)                                                                                                 ccccACAGCTGCCTCCCGTAGGAGTCTGGACCGGGTCTCAGTTCCAGTGTAGCTGGTCATCCTCTCAGACCAGCCAGGGATCGACGCCTTG‑GTGAGCCGTTACCCCACCAAAAAGCTAATCCCATCTGGGCACATCCGATGGCCAGAGGCCCGAAGGTCCCCCTATTTGGTCGTGCGACATTATGCGGTATTAGCTACCGTTTCCAGCA‑GTTCTCCCCCTCCATCAGGCCGCTTCTCAGACATTACTAAcc                                                                      >  2:163316/1‑251 (MQ=11)                                                                                                  ccacctcctgctttccgtagtgggcTGGAACGTTTCTCAGTTCCAGTGTGGCTGGTTAGCCCCTCTGACCTGCCAGGGATCTTTCCCTCG‑GTGAGCCTTTACCCCACCAACAACCTAAGCCCCTTTGGGCGCCTCCGATGGCAAGAGGCCCTAAGTTCCCCCGCCTTGGTCGTGCGACATTATTCGGTATTAGCTACCGGTTCCAGTA‑GTTATTCCCCTCCATCAGGCAGCTTCCCAGACATTTCTCAccc                                                                     <  2:640636/227‑1 (MQ=255)                                                                                                      ctgctgCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTcc                                                                 >  1:215929/1‑251 (MQ=21)                                                                                                      cggctgcCTCCCGGAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCAGCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTcc                                                                 <  1:486832/249‑1 (MQ=255)                                                                                                               cccttAGGGGTCTGGACCGTGTCCCAGTTCCAGTGTTGCTGGTCATCCTCGCAGCCCAGCCAGGGTTCGCCGCCTTG‑GTGAGCCGTTACCCCCCCAACAAGCTAATCCCCTCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTCCCATTA‑GTTACCCCCCTCCATCAGGCGGCTTCCCAGACATTACTCTCCCGTCCGCCACTc                                                          <  2:281069/245‑1 (MQ=14)                                                                                                               cccGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTc                                                          <  1:644920/249‑1 (MQ=25)                                                                                                                        gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCCCAGCCCAGCAAGGGATCGTCGCCTTG‑GTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAagc                                               <  1:526491/251‑1 (MQ=255)                                                                                                                        gTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTAATCCTCTCAGACCAGCCAGGGATCGACGCCTTG‑GAGAGCCGTTACCCCCCCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGACTTGCGACATTATGCGGTATTAGCTCACGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTAATCACCCGTCCGCCACTCGTCAGCGAagc                                               >  1:31453/1‑251 (MQ=11)                                                                                                                          cTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAagc                                               >  1:317475/1‑249 (MQ=25)                                                                                                                            ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAgcagca                                           <  1:558564/251‑1 (MQ=37)                                                                                                                            ggACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAagcagc                                            >  2:690523/1‑250 (MQ=37)                                                                                                                                       cAGTTCCAGTTTGCCTGGTCACCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTcc                                 <  1:544344/250‑1 (MQ=255)                                                                                                                                       cAGTTCCAGTGTGGCTGGTCATCCTCTCAGCCCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCt                                <  1:488798/251‑1 (MQ=255)                                                                                                                                        aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTg                               >  1:600838/1‑251 (MQ=38)                                                                                                                                           tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGtt                             >  1:73227/1‑250 (MQ=39)                                                                                                                                             cagtgggtcgggTCATCCTCCCAGACCAGCCAGGTACCGGCGCCGTG‑GTGAGCCGTTACCCCACCAACAAGCTAAACCCATCTGGGCCCATCCGATGGCCAGAGGCCCGAAGTCCCCCCGCTTTGGTCTGGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGTAAGCTGCTTCCTGTTAcc                          <  2:614983/241‑1 (MQ=14)                                                                                                                                                ggtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGATCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGtt                       <  1:539298/250‑1 (MQ=255)                                                                                                                                                   ggCTGGTCACCCTCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGGGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAGGGCAGGAGGCCCGAAGGTCCCCCTCTTTGGTCTTTCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGa                    <  2:252852/251‑1 (MQ=14)                                                                                                                                                    gCTGGTCATCCCCTCAGACCAGCCAGGGATCGTCGCCTTG‑GTGGGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGAc                   <  1:541335/251‑1 (MQ=255)                                                                                                                                                           aTCCTCCCAGACCAGCCAGGTTTCGGCGCCTTG‑TTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACTTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTTTGGCGACATTATGCGGTATTAGCTACCTTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTCCTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCAtg            <  1:407893/251‑1 (MQ=14)                                                                                                                                                                     aCCAGCCAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCt  >  2:284708/1‑251 (MQ=38)                                                                                                                                                                           |                                                                                                                                                                                                                                                       ACTCCCTTCCTCCCCGCTGAAAGTACTTTACAACCCGAAGGCCTTCTTCATACACGCGGCATGGCTGCATCAGGCTTGCGCCCATTGTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTG‑GTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTA‑GTTATCCCCCTCCATCAGGCAGCTTCCCAGACATTACTCACCCGTCCGCCACTCGTCAGCGAAGCAGCAAGCTGCTTCCTGTTACCGTTCGACTTGCATGTGTTAGGCCT  >  NC_000913/2730700‑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: ‑