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	93.0%	152.3 / 3.0	58	A568A (GCG→GCA)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (22/31);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

TCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGAC  >  NC_000913/13621‑14108                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   tatgctaaaGATAAAAACAGCGGTAAAGAGCAGAAGACCACCATCAAGGCTTCTTCGGGTATGAACGAAGATGAAAACCAGAAAATGGGACGCGACGCAGAAGCTACCGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACGCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTCTCGAGTCTGCACTGa                                                                                                                                                                                                                                                <  2:280888/244‑1 (MQ=255)  ccGCGAAAGATAAAAACAGCTGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACt                                                                                                                                                                                                                                              <  1:24967/251‑1 (MQ=255)       aaaGATAAAAACAGCGGTAATGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACt                                                                                                                                                                                                                                         <  2:446386/251‑1 (MQ=255)                     cGGTAAAGAGCAGAAGATCACCATCAAGGATTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGGCCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCTTCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAACACTGATATCGAGTCTGCACTGACTGCCCTGGAAACTACTCTg                                                                                                                                                                                                                            >  2:162828/1‑250 (MQ=255)                         aaaGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAgg                                                                                                                                                                                                                       >  1:184135/1‑251 (MQ=255)                            gagCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGCCTGCACTGGAAACTGCTCTg                                                                                                                                                                                                                            >  2:434132/1‑243 (MQ=255)                            gagCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTg                                                                                                                                                                                                                            <  1:434132/243‑1 (MQ=255)                                agaagaTCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACACGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACAGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGa                                                                                                                                                                                                                 >  1:412787/1‑250 (MQ=255)                                   agaTCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACaaa                                                                                                                                                                                                             <  1:242887/251‑1 (MQ=255)                                     aTCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGACGAGCTGGTAAAGACTCGCAACCAGGCCGACCATCTGCTGCACAGCCCCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACAGCTATCGAGTATGCACTGACTGCACTGGAACCTGATCTGAACGGAGAAGACAAAGc                                                                                                                                                                                                           >  1:55601/1‑251 (MQ=255)                                         ccATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCt                                                                                                                                                                                                       >  1:389896/1‑251 (MQ=255)                                           aTCAAGGCTTCTTCGGGTCTGAACGAAGAGGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTAt                                                                                                                                                                                                     <  2:184135/251‑1 (MQ=255)                                              aaGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCg                                                                                                                                                                                                   >  2:180697/1‑250 (MQ=255)                                                          gggcTGAACGAAAATGAATACCAGAAAATGGAGCTCGACGCGTAAGCTAACACCGAAGCTTACCGGAATTTTGAAGCACGGGAACAGACTCGCAACCAGGGCCAACATCTGCTAAACAGCAGCCGGAAGCAGGTTGAAGAATCAGGCGACAAACTGCAGGATTACGACAAAACTGTTATCGAGTATGCACTGACTTCACTGAAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTGTCGAACCGAAAAtccc                                                                                                                                                                                      <  2:4742/248‑4 (MQ=255)                                                                   gAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGc                                                                                                                                                                             >  1:47066/1‑251 (MQ=255)                                                                                       ggTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCAATGACTGCACTGGAACCTGCTCTGAAAGGTGAAGACAAAGCCGCAATCGAAGCGAAAATGCAGGAGCTAGCACAGGTTTCCCAgcacctg                                                                                                                                                          >  2:102020/1‑244 (MQ=255)                                                                                               aCGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATc                                                                                                                                                 <  2:265055/251‑1 (MQ=255)                                                                                                cGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCg                                                                                                                                                <  2:369629/251‑1 (MQ=255)                                                                                                cGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCg                                                                                                                                                <  1:392273/251‑1 (MQ=255)                                                                                                  cAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGcc                                                                                                                                              <  2:156271/251‑1 (MQ=255)                                                                                                  cAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCAGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGcc                                                                                                                                              >  1:379888/1‑251 (MQ=255)                                                                                                   aGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGAGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGACGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGAAAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTACCAGAAACTGATGGAAATCGccc                                                                                                                                             >  1:457724/1‑251 (MQ=255)                                                                                                       gcaaaCGCCGAATCGGCCCGTAAGTTTGAAGAGCAGGTGAAGACTCGCAACCATGGCGGCAAGCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCTACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCcagca                                                                                                                                         <  2:265260/248‑1 (MQ=255)                                                                                                            cGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGGTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAAc                                                                                                                                    <  1:328123/251‑1 (MQ=255)                                                                                                                gAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGTACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGc                                                                                                                                                                                                        >  2:201130/1‑179 (MQ=255)                                                                                                                gAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGTACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGc                                                                                                                                                                                                        <  1:201130/179‑1 (MQ=255)                                                                                                                 aaGCTGACCGTTAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGTCAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCTCTGACTGCGCTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGcc                                                                                                                               <  1:497685/251‑1 (MQ=255)                                                                                                                   gccgcccttAAGTTTGAATCGCTGGTACAGACGCGTAACCAGGGCGACCATCTGATGCACAGCACCCGTAAGCAGTGTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCca                                                                                                                             <  2:437706/243‑1 (MQ=255)                                                                                                                      gACCGGAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcag                                                                                                                            <  2:186366/249‑1 (MQ=255)                                                                                                                                      agCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCt                                                                                                          <  1:452338/251‑1 (MQ=255)                                                                                                                                          ggTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGa                                                                                                        <  1:399043/249‑1 (MQ=255)                                                                                                                                             aCAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTAGCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCtt                                                                                                   <  2:170854/251‑1 (MQ=255)                                                                                                                                                                 accaccggctgcACAGCACCCGTAAGCATGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTACTAACGAGGCGGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCTAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAAGCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGGGCTGATGCTGCTGCAAACAACGCGAAAGag                                                                               <  2:373511/244‑2 (MQ=255)                                                                                                                                                                  ccATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTTCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATg                                                                              <  1:417147/251‑1 (MQ=255)                                                                                                                                                                  ccATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATg                                                                              >  2:337141/1‑251 (MQ=255)                                                                                                                                                                   cATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGa                                                                             <  1:91026/251‑1 (MQ=255)                                                                                                                                                                   cATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGATGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATg                                                                              >  1:473907/1‑250 (MQ=255)                                                                                                                                                                               aGCACACGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGAACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATAGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGGTGTCGa                                                                 <  1:367741/251‑1 (MQ=255)                                                                                                                                                                                 cACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCTCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGGCGACg                                                               <  1:102020/251‑1 (MQ=255)                                                                                                                                                                                       tAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAAt                                                         >  1:384199/1‑251 (MQ=255)                                                                                                                                                                                        aaGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGcc                                                                                                                >  1:181042/1‑195 (MQ=255)                                                                                                                                                                                        aaGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAAACGCCCAGCAGCAACATGCCCAGCAGCAGACTGcc                                                                                                                <  2:181042/195‑1 (MQ=255)                                                                                                                                                                                               ttGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGCAGCGAAAATGCAGGAGATGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCACAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTgaaga                                                  >  1:331602/1‑250 (MQ=255)                                                                                                                                                                                                tGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTgaagaa                                                 <  1:337141/250‑1 (MQ=255)                                                                                                                                                                                                 gaagaaGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGt                                               <  1:176207/251‑1 (MQ=255)                                                                                                                                                                                                  acgaagCCGGCGACCAACTGCAGGCTGACGACAAAACTGCTATCGAGTCTGCACTGAATGCACTGGAAACTGCTCTGAAAGGAGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCACAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGACGACGCTGACTTTggagaaaac                                              >  2:321839/3‑248 (MQ=255)                                                                                                                                                                                                             gACAAACTTCCGGCGGACTACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGGGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAAGCGCCCAGCAGCAACATGACCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACaaaaa                                   <  2:206228/251‑1 (MQ=255)                                                                                                                                                                                                             caCAAACTGCAGGCTGACGACAAAACTGCTATCGATTCTGCACTAACTGCACTGTAAACGACACAGAAAGGTGAAGACAAATCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTATGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAGTTTGAAGAAGTCAAAGACaaaaa                                   <  2:55601/250‑1 (MQ=255)                                                                                                                                                                                                               cAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAaat                                 >  1:432774/1‑251 (MQ=255)                                                                                                                                                                                                                     gCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGc                           <  1:184450/251‑1 (MQ=255)                                                                                                                                                                                                                                       tGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGCAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGATGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAAtt          >  2:363816/1‑250 (MQ=255)                                                                                                                                                                                                                                         cTATCGAGTCTGCACTGACTTCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAACGAAAATGCAGGGGCTGGCACAGGTTTCCAAGAAACTGATGGAAATCACCCAGCAGCAAAATGCCCAGCAGCAGACTGCAGATGCTGATGCTTCGACCAACAACGCGGAAGATGACGTTGTTGTCGATGCTGAGTTTGAAGAAGTCAAAGACAAAAAATAATCGCACTATAAAGGGGTAGTtat        >  2:33873/1‑250 (MQ=255)                                                                                                                                                                                                                                         cTATCGAGTCTGAACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACACAGCCGCTATCGAAGCGAAAATGCAGGCGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCCACATGCCCAGCAGCAGCCTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACgggaacttat        >  2:332111/1‑243 (MQ=255)                                                                                                                                                                                                                                           aTCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATCCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACt     <  2:4409/251‑1 (MQ=255)                                                                                                                                                                                                                                           aTCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGCAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGTTAATTATACt     >  2:46244/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACTCTGAATTTGAAGAAGGCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACt     <  1:332111/251‑1 (MQ=255)                                                                                                                                                                                                                                           aTCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATAc      <  1:427152/250‑1 (MQ=255)                                                                                                                                                                                                                                              gAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGac  <  1:63067/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   TCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGAC  >  NC_000913/13621‑14108

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