Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,696	C→T	A221T (GCC→ACC)	yaaW ←	putative enzyme‑specific chaperone YaaW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,696	0	C	T	97.9%	149.0 / NA	49	A221T (GCC→ACC)	yaaW	putative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/21);  minor base A (1/0);  total (27/22)
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

GGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATC  >  NC_000913/10455‑10933                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               ggCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGAGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGt                                                                                                                                                                                                                                      >  1:218686/1‑251 (MQ=255)     cAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCAcg                                                                                                                                                                                                                                  >  1:36221/1‑251 (MQ=255)        cGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGa                                                                                                                                                                                                                               <  1:4567/251‑1 (MQ=255)             ccGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCAGGAGCAGGCAGGCGAATTGCAGTACGGTTG‑GAATCGTCACTCGATAGgc                                                                                                                                                                                                                          >  2:550080/1‑251 (MQ=255)               ggATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGgcgc                                                                                                                                                                                                                        >  1:102537/1‑251 (MQ=255)                    aGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGAGTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGcc                                                                                                                                                                                                                    <  1:314595/250‑1 (MQ=255)                      ccAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGCATGGTGAAATAATTTCCCTGAATAACTGTAGTGATTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCt                                                                                                                                                                                                                 >  1:333514/1‑251 (MQ=255)                       cAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTg                                                                                                                                                                                                                <  1:249035/251‑1 (MQ=255)                        aaCTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGCATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTg                                                                                                                                                                                                                >  1:96488/1‑250 (MQ=255)                              gtgtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCtt                                                                                                                                                                                                         >  1:204479/1‑251 (MQ=255)                               ggtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCttt                                                                                                                                                                                                        <  1:309552/250‑1 (MQ=255)                                                        atcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGAATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTAc                                                                                                                                                                                >  2:518587/1‑250 (MQ=255)                                                              gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCt                                                                                                                                                                         >  1:556255/1‑251 (MQ=255)                                                              gaaggaatgGTAAATTAATTTCCCTGAATAACTGTAGTGGTTTCAAGGCTCGGCATAATAACCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTGTCCGGTATCAGACCTGAGTGGCGCTAACCAGCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCt                                                                                                                                                                         <  2:13858/247‑1 (MQ=255)                                                                         gAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgc                                                                                                                                                              >  1:386620/1‑251 (MQ=255)                                                                           aataatTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGAGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGGTTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAACCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCAccgcgc                                                                                                                                                            >  2:450704/1‑251 (MQ=255)                                                                                      tGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGACAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCTTTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCCCCGCGCAGCaa                                                                                                                                                  >  2:526644/1‑250 (MQ=255)                                                                                        aaTAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGACGCGCACCAAAcc                                                                                                                                               >  2:202641/1‑251 (MQ=255)                                                                                                            ggCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATtg                                                                                                                           >  1:449754/1‑251 (MQ=255)                                                                                                                cgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATtgctgc                                                                                                                       >  2:217323/1‑251 (MQ=255)                                                                                                                  gCATTATAATCAGTCAGTGGGGCAGTGTCTACGATCTTTGTAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCAGCCGGCGCAGGCATGCGATTTGCAGTACGGTTG‑GAAACGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAGCTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCTCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGt                                                                                                                     <  2:557711/251‑1 (MQ=255)                                                                                                                    ataataATCAGCCCGTGGGGCAGAGGCTAAGATCATTTGAGGGGAAAATGAAAATTTTCCGCGGTTTCCGGTATAAGACCTGAGTGGAGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAAACGTCACGCGATAGGAGCTGCCTCTGAAAGATTTCAACCCATTTCGCGCTGCCACTACAGGGCATCCCCGGCCCGCGCCGCGCAGCAAACCATGCCCACGTACGCGCATTGCTGCGTgg                                                                                                                   >  2:292794/1‑251 (MQ=255)                                                                                                                     taataaTCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTCTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTCACCCCa                                                                                                                                                                                                 >  2:48558/1‑172 (MQ=255)                                                                                                                     taataaTCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCa                                                                                                                                                                                                 <  1:48558/172‑1 (MQ=255)                                                                                                                                      ggCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGTCCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTc                                                                                                 <  2:173867/251‑1 (MQ=255)                                                                                                                                         aGTGTCTACGATCTTTTGAGGGGAAAATGAACATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGt                                                                                              >  1:431638/1‑251 (MQ=255)                                                                                                                                         aGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGt                                                                                              >  1:193590/1‑251 (MQ=255)                                                                                                                                              cTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCt                                                                                         >  1:353704/1‑251 (MQ=255)                                                                                                                                                tcgagcTTTTGAGGGGAAAATGAAATTTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCCAAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGGGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTg                                                                                        <  2:281273/245‑1 (MQ=255)                                                                                                                                                aCGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTgg                                                                                       <  1:66394/251‑1 (MQ=255)                                                                                                                                                 cGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCCGTTGGCTgg                                                                                       >  2:138666/1‑250 (MQ=255)                                                                                                                                                   aTCTTTTGAGGGGAAACTGAAAATTTTCCCCGGATTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGATG‑GAATCGTCACGCGCTAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGAGCGTAAAATGCGGGTCAGTTGGCTGGaa                                                                                     >  2:556991/1‑250 (MQ=255)                                                                                                                                                              ggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTGAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCGCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGac                                                                         <  1:519087/251‑1 (MQ=255)                                                                                                                                                                     gaaaatTTTCCCCGGTTTCCGGTAGCAGTCCTGAGTGGCGCTAACCATCCGGCGCAGGCATGCGATTTGCAGTACGGTTG‑GAAGCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCTAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGGAAAATGCGGGTCAGTCGGCTGGAAAGCAAATGCGACACACCtt                                                                  <  1:327197/251‑1 (MQ=255)                                                                                                                                                                       aaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGTTTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCtttt                                                                <  2:193590/251‑1 (MQ=255)                                                                                                                                                                         atTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGc                                                              <  1:14909/251‑1 (MQ=255)                                                                                                                                                                               cccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCaata                                                         >  1:68121/1‑250 (MQ=255)                                                                                                                                                                                  ggTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAAttt                                                     <  1:202641/251‑1 (MQ=255)                                                                                                                                                                                         ggTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTtc                                              <  1:383797/251‑1 (MQ=255)                                                                                                                                                                                          gTATCAGACCTGAGTGGCGCTATCCATCCGGCGCAGGTAGGCGATTTGTAGTACGG‑TGTGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAGTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCttgca                                             <  2:200450/251‑4 (MQ=255)                                                                                                                                                                                           tATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAGACCATGCCCAAGTACGCTCATTGCTGCGTGGGCGCGGAAAATGCGGGTCAGTTGGATGGAAGGCAACTGCGACCCACCTTTTg                                                               >  2:471215/1‑232 (MQ=255)                                                                                                                                                                                           tATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTg                                                               <  1:471215/232‑1 (MQ=255)                                                                                                                                                                                                       gTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCaataa                                                        >  1:56740/1‑227 (MQ=255)                                                                                                                                                                                                       gTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCaataa                                                        <  2:56740/227‑1 (MQ=255)                                                                                                                                                                                                             cTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTTCGCAAAATGCGGGTCAGTTTGCTGGAAAGCAAATGCGACTAACCATTTGCCCATAAATTGTCTTTCATCAGCAGCGGCagcagccc                          >  2:135599/1‑249 (MQ=255)                                                                                                                                                                                                                 cctcccGGCGCAGGCCGGCTAGTTGCAGTACGGTTG‑GAACCGTCGCGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAGTTAGTGCCGCACCTACAGGGCCTCCCATCCCCGCGCCGCGCAGCAAACCATGCCCAGGTACGCTCAGTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTGGGCTGGAAAGCAAATGCGACACGCCGTTTGCCAATAATTTGTCTTGCATCAGCAGCGTCAGCATCTCTTcc                      <  2:333514/247‑1 (MQ=255)                                                                                                                                                                                                                        gagcaggcagTCGATTTGCAGTGCGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTTCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTTAAATGCTGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCAt               <  2:349152/249‑1 (MQ=255)                                                                                                                                                                                                                                 ggcGATTTGCAGTACGGTTG‑TAATCGTCACGCGATGGGCGCTGCCGCTGAACTCTTTAACCCAATTTAGGGCCGCATCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGGCAGTTGGCTGGAAAGCAAATGCGACACCCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCAACCTgg      <  2:431576/251‑1 (MQ=255)                                                                                                                                                                                                                                     aTTTGCAGTACGGTTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCCTTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATc  >  1:543436/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               GGCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTG‑GAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATC  >  NC_000913/10455‑10933

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: ‑