Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,866	G→A	A568A (GCG→GCA)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,866	0	G	A	95.5%	185.7 / ‑1.8	68	A568A (GCG→GCA)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (27/37);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

CCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGAC  >  NC_000913/13622‑14108                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   ccgctgaatataaaaccttcgtttaaGAGCAAAAGATCACCAGCACGTATTCTGCTGATGTGAACTAGGATGTAATCTAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGTCCGGAAGAGTAAAGATCTGGTACAGAATCGCAGCCATGTCTACTATCTGCTGCATATCACCAGTAAGCATGTTGAAGAGCCATGCGTCTAACTGCCGGCGTACTACAAAACTGCTACCGAGTCTGCACTGAcg                                                                                                                                                                                                                                              <  2:317077/227‑2 (MQ=255)  cgcgAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCactgactg                                                                                                                                                                                                                                             >  1:68820/1‑251 (MQ=255)               aaCAGCGGTAAAGAGCAGAAGATCAACATCAAGGCTTCTTCTGGTCTGAACGCAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGc                                                                                                                                                                                                                                >  1:225056/1‑251 (MQ=255)                        aatgatcAGAAGAACACCAGCAAGGCTTCTTCTGGTCTGAAGGAAGAGGAAATCCAGAAAAGGGGACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGTATCGCAAGCAGGGTGACCATCTGCTGCCCAGCACCCGTAAGCAGGTTGAATAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATTGAGTCTGCACTGACTGCACTGGAAACTGCGCTGAAAgt                                                                                                                                                                                                                       <  2:506480/245‑2 (MQ=255)                         aaGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACAATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGGGACAACCTGCCGGATGACGACAAAACTGCTATCGAGTCTGCACTGACTTCACTGGAAACTGCTCTGAAAggg                                                                                                                                                                                                                      >  2:578830/1‑250 (MQ=255)                              cAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGa                                                                                                                                                                                                                 <  1:29959/251‑1 (MQ=255)                                  agaTCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACaa                                                                                                                                                                                                              >  1:456660/1‑250 (MQ=255)                                      cACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTTCGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGcc                                                                                                                                                                                                          <  1:472979/250‑1 (MQ=255)                                            cAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATc                                                                                                                                                                                                    <  1:269021/250‑1 (MQ=255)                                               ggCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCCGGTTGACGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGAGATCGAGTCTGCACTGACTGCCCTGGAAACTGCTCTGAAAGGCGAAGCCAAAGCACCTATCGacg                                                                                                                                                                                                >  1:386843/1‑249 (MQ=255)                                               ggCTTCTCCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCAGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGAGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGACTGACGACAAAACTGCTATCGAGTCTGCACTGACAGCACTGGACACGGCGCTGAAAGGTGAAGCCAAAGCCGcattcaa                                                                                                                                                                                                  >  2:296701/1‑243 (MQ=255)                                                     ttctGGTCTGAACGCAGATGAAATCCAGAAAATGGTACGCCACGCAGAAGCTAACGACGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCACCCAGGGCGCCCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTACCGGCTGCCGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTCTCGAAGCGaaa                                                                                                                                                                                           >  2:258306/1‑250 (MQ=255)                                                      tctGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGaaaa                                                                                                                                                                                          <  1:30160/250‑1 (MQ=255)                                                         ggTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATg                                                                                                                                                                                        <  2:362933/249‑1 (MQ=255)                                                           tCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGAAAAAACTACTATCGAGTCTGCACTGACTGCACTGGTAACTGCTCTGAAAGGTGAAGACAAAGACGCTATCGAAGCGAAAATGCAgg                                                                                                                                                                                    >  1:319104/1‑251 (MQ=255)                                                                aCGAAGAAGACATCCAGAAAAGTGTACGCGACGCACAAGCTAACGCCGAAGCTGCCCGTAAGGTTGAATAGCTGGTACAGCCGCGCAACCGGTGCGACCATCTGCTGCACAGCACACGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCGCGGAAAGGTTAAGACAAATCCGCTATCGAAGCGAAAATGCAGGATCTg                                                                                                                                                                               <  2:353354/251‑1 (MQ=255)                                                                  gAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGATTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGATGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGc                                                                                                                                                                             >  1:659647/1‑251 (MQ=255)                                                                    aGATGAAATCCAGAAAATGGAACGCGCCGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGCCGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGcac                                                                                                                                                                           >  1:39778/1‑251 (MQ=255)                                                                     gATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGcaca                                                                                                                                                                          <  1:513955/251‑1 (MQ=255)                                                                      aTGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGcac                                                                                                                                                                           >  1:419009/1‑249 (MQ=255)                                                                         aaaTCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGag                                                                                                                                                                                  >  1:48102/1‑238 (MQ=255)                                                                         aaaTCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGag                                                                                                                                                                                  <  2:48102/239‑2 (MQ=255)                                                                             ccAGAAAAGGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTccc                                                                                                                                                                  <  1:370012/251‑1 (MQ=255)                                                                                  aaaTGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGaa                                                                                                                                                              >  1:635245/1‑250 (MQ=255)                                                                                              aCGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATc                                                                                                                                                 <  1:478475/251‑1 (MQ=255)                                                                                                  aGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGccc                                                                                                                                             >  1:261512/1‑251 (MQ=255)                                                                                                             ccGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACAt                                                                                                                                  >  1:398963/1‑251 (MQ=255)                                                                                                              cGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACAtc                                                                                                                                 >  1:582836/1‑250 (MQ=255)                                                                                                                aaggtgacctgaaGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGcc                                                                                                                               <  2:189840/240‑1 (MQ=255)                                                                                                                 aGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCAGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGccc                                                                                                                              <  2:84141/251‑1 (MQ=255)                                                                                                                    tGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcagc                                                                                                                           >  1:324845/1‑251 (MQ=255)                                                                                                                     gACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTcc                                                                                                                                                                   <  1:589569/210‑1 (MQ=255)                                                                                                                      aCCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGtt                                                                                                                                                                      >  2:202782/1‑206 (MQ=255)                                                                                                                      aCCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGtt                                                                                                                                                                      <  1:202782/206‑1 (MQ=255)                                                                                                                            aggttGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCCGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACt                                                                                                                   <  2:194659/248‑1 (MQ=255)                                                                                                                               ttGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGc                                                                                                                 <  1:566718/250‑1 (MQ=255)                                                                                                                               ggtaaGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGCAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGc                                                                                                                 <  2:658765/247‑1 (MQ=255)                                                                                                                                            tcAGACTCGCATCCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCt                                                                                                    <  2:582836/249‑1 (MQ=255)                                                                                                                                                  ccGCAACCAGGGCGACCATCTGCTGCACTGCACCCGTAAGCAGGTTGAAGAAGCAGGCGCCAAACGGCCGTGTGACGACAAAACTGGTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCTAAAATGCAGGAGCTGGCGCAGGTGTCCCAGAAGCTGATGGAAACCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGGGCTGATGTTTCTGCaa                                                                                             <  2:39778/250‑1 (MQ=255)                                                                                                                                                       aCCAGGGCGACCATCTGCTTCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACTACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAaacaac                                                                                        <  1:3075/251‑1 (MQ=255)                                                                                                                                                            ggCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGa                                                                                    >  1:449286/1‑250 (MQ=255)                                                                                                                                                               gACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGa                                                                                <  1:254025/251‑1 (MQ=255)                                                                                                                                                                aCCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGaaa                                                                                                                                                             >  1:373150/1‑173 (MQ=255)                                                                                                                                                                  cATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATg                                                                              >  1:50874/1‑250 (MQ=255)                                                                                                                                                                          gcACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGAtgt                                                                       <  1:258306/249‑1 (MQ=255)                                                                                                                                                                             caGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                  <  1:173139/251‑1 (MQ=255)                                                                                                                                                                             caGCACCCGTAAGCAGGGTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAAATGCTATCGAGTCTGCACTTACTGCACGGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                  <  2:419009/251‑1 (MQ=255)                                                                                                                                                                                 aCCCGTAAGCAGGTTAAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGAAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGc                                                              <  2:171408/251‑1 (MQ=255)                                                                                                                                                                                  cccGAAAGCAGGTTGAAGTAGCATGCGACAAACAGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAATCGAAAATGCAGGAGCTGTCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCATCAACATGCCCAGCAGCAGGCTGCCGGTGCTGATGCTTCTGCAAACAACGCTAAAGATGGCGATGTTTTCTACGCt                                                             <  2:68820/251‑1 (MQ=255)                                                                                                                                                                                    cGTAAGCAGGTTGAAGAATCATGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTTCACGGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGa                                                           <  2:449286/251‑1 (MQ=255)                                                                                                                                                                                          cAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTg                                                      >  1:322042/1‑250 (MQ=255)                                                                                                                                                                                           aGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTgaa                                                    <  1:190698/251‑1 (MQ=255)                                                                                                                                                                                            ggTTGAAGAAGCAGGCGACAAACTGCCGGCGGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCGTCTGCAAACAACGCGAAAGATGACGAGGTTGTCGACGCTGAATTTgaag                                                   <  2:289884/251‑1 (MQ=255)                                                                                                                                                                                                     agCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGAATGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCAGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCACCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGATGATGCCTATGCAAACAACGCGAACGATGACGCTGGTGTCGTCCCTGAATTTGAAGAAGTcca                                            >  2:319034/1‑247 (MQ=255)                                                                                                                                                                                                          gCGAGAAACTGCCGGCTGACGACAAAACGGCTATCGAGTCTGCCCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTGCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACaaa                                     <  2:13123/251‑1 (MQ=255)                                                                                                                                                                                                          gCGACAAACTGCCGGCTGACGCCAAAACTGCTATCGAGTCTGCACTTACTGCACTGGAAACTGCGCTGAAAGGAGAAGACAAAGCCGGTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCTAAGACa                                       <  2:599766/249‑1 (MQ=255)                                                                                                                                                                                                                 aCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                                                                              >  1:292296/1‑187 (MQ=255)                                                                                                                                                                                                                 aCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCa                                                                                              <  2:292296/187‑1 (MQ=255)                                                                                                                                                                                                                 aCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAaataat                              >  1:281858/1‑251 (MQ=255)                                                                                                                                                                                                                        gCTGACGACAAAACTGCTATCGAGTCTGCACTGGCTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCta                       <  1:356104/251‑1 (MQ=255)                                                                                                                                                                                                                           gacgacAAAACTGGAATCGAGTCTTCACCGACTGCACTGGAAACTGCTCTGAAAGTTTAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAAGTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATaa                    <  2:580144/251‑1 (MQ=255)                                                                                                                                                                                                                           gacgacAAAACTGCTATCGAGCCTGCACTGACTGCACTGGAAACTGCACTGAAAGGTCAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCGGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGGTTCTGCAAACAACGCGAAAGATGACGATGTGGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAAGCGCCCTATaa                    <  2:456660/251‑1 (MQ=255)                                                                                                                                                                                                                            acgacAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGAAACCACCGCGAAAGAGGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAATTAATAGCCCTATaa                    >  2:411133/1‑250 (MQ=255)                                                                                                                                                                                                                            acgaaaaaccttcTATCTAGTCTGCACTTACTGCACTGGAAACTGCTTTGAAAGGTGAAGACAAAGCCGCTATCGTAGCGAAAATGCAGTAGCTTGCACAGGTTTCTCAGAAACTGATGGAAATCTCCCAGCAGCAACATGCCCAGCAGCAGACTTCCGGTGCTGATTCTTCTTCAAACAACTCGAAAGATGACGATGGTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATaaa                   <  2:659647/239‑1 (MQ=255)                                                                                                                                                                                                                                  aaaCTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                  >  2:105895/1‑198 (MQ=255)                                                                                                                                                                                                                                  aaaCTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCg                                                                  <  1:105895/198‑1 (MQ=255)                                                                                                                                                                                                                                       gCTATCGAGTCTGCACTGATTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAAtt          <  1:470773/249‑1 (MQ=255)                                                                                                                                                                                                                                         tATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATAc      >  1:432586/1‑251 (MQ=255)                                                                                                                                                                                                                                             gAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGac  >  1:436204/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGAC  >  NC_000913/13622‑14108

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