Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,194,362	(G)6→7	coding (1499/1662 nt)	yqiK →	flotillin family inner membrane protein YqiK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,194,362	1	.	G	94.9%	166.0 / 2.3	59	coding (1499/1662 nt)	yqiK	flotillin family inner membrane protein YqiK
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (19/37);  total (20/39)
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.61e-01

AGCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑G‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTT  >  NC_000913/3194114‑3194610                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         aGCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑ggg‑gg‑ggc                                                                                                                                                                                                                                                         >  1:610478/1‑249 (MQ=255) aGCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAAAGATGCTATCAACGTACTTACTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCACTCGACGGTATCAAGATTATTCAGGTCGATGGTATGAATCGTGGAGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑ggg‑gg‑ggc                                                                                                                                                                                                                                                         >  1:433758/1‑249 (MQ=255)  gCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCACCGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑ggg‑gg‑ggcg                                                                                                                                                                                                                                                        >  1:517558/1‑248 (MQ=255)       aGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTTCTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTAGCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAAc                                                                                                                                                                                                                                                   <  1:31193/251‑1 (MQ=255)       aGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCAAAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATGCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGAGGGTCTGAATCGTGGCGGCGCTGCGGGTGATGGAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAAc                                                                                                                                                                                                                                                   <  1:509550/251‑1 (MQ=255)            cacGTAAAAAGGGGCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGAGGAACAAACCAGCCTTAAATTCAAACTGTCCTTGTTGCAGGCGCTGCCGGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTAGTCAGGTCGATGGTCGGAATCGTGGCGGCGCTGCGGGTGATGCAAACACCGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGc                                                                                                                                                                                                                                              <  2:369024/251‑1 (MQ=255)               gTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGa                                                                                                                                                                                                                                           <  1:104033/251‑1 (MQ=255)                 aaaaaGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAAc                                                                                                                                                                                                                                         >  1:316243/1‑251 (MQ=255)                  aaaaGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAAc                                                                                                                                                                                                                                         >  1:363921/1‑250 (MQ=255)                     aGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGc                                                                                                                                                                                                                                     <  1:38910/251‑1 (MQ=255)                        gTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCAtt                                                                                                                                                                                                                                  <  1:263960/251‑1 (MQ=255)                                   gcagaagcaCAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGctc                                                                                                                                                                                                                       <  1:234686/251‑1 (MQ=255)                                      gaagcaCAACGTGCGCTTAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGctctct                                                                                                                                                                                                                    <  1:422853/251‑1 (MQ=255)                                       aagcaCAACGTGCGCTGAACGATGCTATCAACGTCCTTTCTGATGAACAAACCAGCCTGAAATTCAAACTGGCCTTGTTGCAGGCGAGGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAAGCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGctctctc                                                                                                                                                                                                                   <  1:348475/251‑1 (MQ=255)                                        agcaCAACGTGCGCTGCACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGAAATGT‑GGG‑GG‑GGCGGAAACCTGGGGGAACAAGCATATTCGGCCGctctctct                                                                                                                                                                                                                  >  2:267576/1‑251 (MQ=255)                                        agcaCAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCATCCGcactctct                                                                                                                                                                                                                  >  1:75625/1‑251 (MQ=255)                                         gcaCAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCtt                                                                                                                                                                                                                 <  2:118855/251‑1 (MQ=255)                                         gcaCAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGAGGCAAACACGGGTAATGT‑GGG‑GGTG‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCtt                                                                                                                                                                                                                 <  1:258617/251‑1 (MQ=255)                                           acaACGTGCGGTGAACGATGCTATCAACGTACTGTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGGCGGTATCAAGATTATTCAGTTCGATGGTCTGAATCGTTGCGGCTGTGCGTGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCTGGACAAGCATTTACAGCCGCTCTCTCTTAc                                                                                                                                                                                                               <  2:574295/251‑1 (MQ=255)                                            caACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGCACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGATGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGACCAAGCATTATCAGGCGCGCGCTCTTAcc                                                                                                                                                                                                              >  2:425432/1‑251 (MQ=255)                                                   cgcTGAACGATGCTATAAACGTACTTTCTGATGACCAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGccccc                                                                                                                                                                                                        >  2:382700/1‑246 (MQ=255)                                                          cGATTCTATCAACGTACGTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCg                                                                                                                                                                                                <  1:345043/251‑1 (MQ=255)                                                            aTGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCt                                                                                                                                                                                              <  1:553247/251‑1 (MQ=255)                                                                       gTACTTTCTGATGAACAACCCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGAACCCCGGCACCGCTGATTGACTCCt                                                                                                                                                                                   >  2:117973/1‑251 (MQ=255)                                                                         aCTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGTTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCtt                                                                                                                                                                                  <  1:44559/250‑1 (MQ=255)                                                                         aCTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGCCTCCTTg                                                                                                                                                                                 >  2:449016/1‑251 (MQ=255)                                                                                    aaCAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCTGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGGGATGCAAACACGGGGACTGT‑GGG‑GG‑GGCGGAAACTTGGCGGAACAAGACTTATCAGCCGCTCTCTCTTCACGCAACCAGGCACCGCTGATTGACTCCTTGCTCAATGAAAt                                                                                                                                                                      >  2:57465/1‑251 (MQ=255)                                                                                       aaaCCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTgg                                                                                                                                                                   <  1:449016/251‑1 (MQ=255)                                                                                           cAGCCTTAAATTCAAACTGGCATTGGTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGACAATCGACGGTATCAAGATTCTTCAGGTCGATGGTCTGAATCGTGGCGACGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTATTACCGCCCACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCgat                                                                                                                                                               >  1:289739/1‑249 (MQ=255)                                                                                             gCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCCCACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTc                                                                                                                                                             >  1:298254/1‑251 (MQ=255)                                                                                                 tAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAg                                                                                                                                                           <  1:20472/249‑1 (MQ=255)                                                                                                      tCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCCTTGAGCCGATGAAGTCAATAGACGGTATCCAGATTATTCAGGCCGATGGTCTGACTCGTGGCGGCGCTGCGGGTGCTGCAACCACGGGTAATGT‑GGG‑GG‑GGCGGAAACCCGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCCATGAAATTGGCGTTTCAGGCGGCTc                                                                                                                                                    >  1:86577/1‑251 (MQ=255)                                                                                                            gggCCTTGTGGCATGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGG‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCGCAATGAAATTGGCGTTTCAGGCGGCTCACTggc                                                                                                                                              <  2:610478/250‑1 (MQ=255)                                                                                                             ggCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTggcg                                                                                                                                             >  1:420842/1‑251 (MQ=255)                                                                                                                cTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTggcggc                                                                                                                                           >  1:189524/1‑250 (MQ=255)                                                                                                                                      gTAATAGAGAAATCCGTGGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTcaa                                                                                                                    <  2:17093/251‑1 (MQ=255)                                                                                                                                         aTAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTcaaca                                                                                                                  >  1:122689/1‑250 (MQ=255)                                                                                                                                            gagaAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACATGCATTATCAGCCGCTCTCTCTGACCGCTCACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCTCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTc                                                                                                              <  1:47889/251‑1 (MQ=255)                                                                                                                                                   ccGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCg                                                                                                        <  1:382700/250‑1 (MQ=255)                                                                                                                                                    cGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACATGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGTCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCgaa                                                                                                      <  1:565337/251‑1 (MQ=255)                                                                                                                                                    cGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACGGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCgaa                                                                                                      <  1:402707/251‑1 (MQ=255)                                                                                                                                                            cGATCAAGGCAATCGACGGTAGCAAGATTATTCAGGTCGATGGGCTGAATCGTGGCGGCGCTGCGGGGTATGCAAACACGGGTAAGGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGTCTCACTGTCGGCATGGACTTCACCCTTACCCTCAACAACTCCCGTCGAAGAAAAAGc                                                                                              <  2:433758/251‑1 (MQ=255)                                                                                                                                                              aTGAAGTCATTCGCCGGTAACAAGATTAGTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGTAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGAGTGACTCCTTGCTCAATGAATTTGGCGTTTCAGGCGGCTCACTGTCGGCATTGACTTCACCCTTAACCTCAACAACGCCCGTCGAAGAAAAAGc                                                                                              <  2:316243/249‑1 (MQ=255)                                                                                                                                                                gAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCagag                                                                                          <  2:19514/251‑1 (MQ=255)                                                                                                                                                                 aaGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGt                                                                                         >  1:119913/1‑251 (MQ=255)                                                                                                                                                                 aaGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGt                                                                                         <  1:517338/251‑1 (MQ=255)                                                                                                                                                                       aTCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAAccc                                                                                   >  1:193926/1‑251 (MQ=255)                                                                                                                                                                          gACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCctct                                                                                <  1:267576/251‑1 (MQ=255)                                                                                                                                                                                  cAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGa                                                                         <  1:538476/250‑1 (MQ=255)                                                                                                                                                                                   aaGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGATGCGGGAGATGCACACACGGGTAATGTCGGC‑GG‑G‑CGGAAACCTGGCGGACCAAGCATTATCAGCCGCTCTCTCTTACCACACACAGGCACCGCTAATTGACTACTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCAATGGCGGCATTGAATTCACCCTTAACCTCAACAACTCCCGTCGAAGAACAAGCAGGGTCAACCCTCTTCCAGGaat                                                                       >  1:603036/1‑251 (MQ=255)                                                                                                                                                                                       ttattCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCCCTGGCGGCATTGACTTCAcc                                                                                                                               >  2:219701/1‑191 (MQ=255)                                                                                                                                                                                       ttattCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCAcc                                                                                                                               <  1:219701/191‑1 (MQ=255)                                                                                                                                                                                           tCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAGCACGGGTAAGGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCATGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACGCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCt                                                                <  2:99924/250‑1 (MQ=255)                                                                                                                                                                                                       gcTGAATCGTGGCGGCGCTGCGGGGGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATGAGCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACGTCACCCTTACCCTCAACAACTCCCGTCGAAGAAAAAGCAGAGGAAACCCTCTGCCATGAATAATCCCCGGCCCAAAAAGcc                                                   <  2:268455/250‑1 (MQ=255)                                                                                                                                                                                                         ggAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGG‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGGCTCCTTGCTCAATGAAAGTGTCGTGTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCgg                                                 <  2:517558/250‑1 (MQ=255)                                                                                                                                                                                                               gtggtgGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCATCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAATCCCTCTTCCAGTAATAATCCGTGGCCCAAAAAGCCGGGGAttt                                           <  2:95939/251‑1 (MQ=255)                                                                                                                                                                                                                      cgcTGCGGGTGATGCAAACCCGGGTAATGT‑GGGTGA‑G‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCGTGCTCAATGAAATGGGCGTTTCAGACGGCTCACTGGCGGCATTGGCTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTGAAGCCTCTTCCACGAACGATCCCTGGCCCAAAAAGCCGGGGATATTGTatat                                    <  2:89446/251‑1 (MQ=255)                                                                                                                                                                                                                      cgcTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTatat                                    <  1:27222/251‑1 (MQ=255)                                                                                                                                                                                                                                gaggCAAACACGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAATCAGAGTAAACCCTCTTCCATGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTcc                           <  2:363921/247‑1 (MQ=255)                                                                                                                                                                                                                                   gCAAACTCGGGTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTa                       <  2:189524/251‑1 (MQ=255)                                                                                                                                                                                                                                           gggTAATGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCTCACACAGGCCCCGCTGATTGACGCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTGCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTTCGTTCCTCTAAAAAGTGc               <  2:75625/251‑1 (MQ=255)                                                                                                                                                                                                                                               aaTGT‑GGG‑GG‑GGCGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAAt           >  1:517781/1‑251 (MQ=255)                                                                                                                                                                                                                                                    gggg‑gg‑g‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                 <  2:138200/191‑1 (MQ=255)                                                                                                                                                                                                                                                    gggg‑gg‑g‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATccc                                                                 >  1:138200/2‑192 (MQ=255)                                                                                                                                                                                                                                                         gg‑g‑CGGACACCTGGCGGTACACGCATTATCAGCCGCTCTCTCTTACCGCAAACAGGCACCGCTGATTGACTCCTTGCTTAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCCTACCCTCAACCACTCCCGTCGCAGACAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt  >  2:243055/1‑251 (MQ=38)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         AGCCGAAGCTACACGTAAAAAGGGTCTGGCGGAAGCAGAAGCACAACGTGCGCTGAACGATGCTATCAACGTACTTTCTGATGAACAAACCAGCCTTAAATTCAAACTGGCCTTGTTGCAGGCGCTGCCTGCGGTAATAGAGAAATCCGTTGAGCCGATGAAGTCAATCGACGGTATCAAGATTATTCAGGTCGATGGTCTGAATCGTGGCGGCGCTGCGGGTGATGCAAACACGGGTAATGT‑GGG‑GG‑G‑CGGAAACCTGGCGGAACAAGCATTATCAGCCGCTCTCTCTTACCGCACACAGGCACCGCTGATTGACTCCTTGCTCAATGAAATTGGCGTTTCAGGCGGCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTT  >  NC_000913/3194114‑3194610

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