Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	495,680	C→T	F187F (TTC→TTT)	htpG →	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	495,680	0	C	T	81.4%	132.1 / 23.2	60	F187F (TTC→TTT)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (24/24);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

AATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAACCGTGTTGAA  >  NC_000913/495439‑495926                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       aaTCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGAc                                                                                                                                                                                                                                               <  1:145781/251‑1 (MQ=255)  aTCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACt                                                                                                                                                                                                                                              <  1:173384/251‑1 (MQ=255)    cccTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTgg                                                                                                                                                                                                                                            <  1:472362/251‑1 (MQ=255)     ccTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGgc                                                                                                                                                                                                                                           <  1:283799/251‑1 (MQ=255)       tGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGgcgc                                                                                                                                                                                                                                         >  1:376148/1‑251 (MQ=255)            tCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCgcgtgc                                                                                                                                                                                                                                     <  1:249656/250‑1 (MQ=255)                  cAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGCCAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAACGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCcat                                                                                                                                                                                                                              >  2:54351/1‑251 (MQ=255)                       gAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCatcatc                                                                                                                                                                                                                          >  2:391683/1‑250 (MQ=255)                       gAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCAtcatca                                                                                                                                                                                                                         <  2:331723/251‑1 (MQ=255)                             cagccagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAAt                                                                                                                                                                                                                   >  2:242512/1‑251 (MQ=255)                                          ggTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCatat                                                                                                                                                                                                      <  1:235149/251‑1 (MQ=255)                                             cAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATcgc                                                                                                                                                                                                   <  1:166940/251‑1 (MQ=255)                                                  tGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGc                                                                                                                                                                                              >  2:463503/1‑251 (MQ=255)                                                   ggTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGGTCCATCATCAGCAAATACTCCGACCATATCGCGCTGcc                                                                                                                                                                                             >  2:266311/1‑251 (MQ=255)                                                      gtTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGGGCGTGCGTGCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGt                                                                                                                                                                                          <  1:24618/251‑1 (MQ=255)                                                         ggTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCTTGGTACTGAAATCACCCTGCATCTTCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTag                                                                                                                                                                                        <  1:65655/250‑1 (MQ=255)                                                                aCTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATc                                                                                                                                                                                   <  1:54351/248‑1 (MQ=255)                                                                     gCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAAc                                                                                                                                                                            <  1:463503/250‑1 (MQ=255)                                                                         ttATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCgcgtgcgt                                                                                                                                                                                                                                   <  2:421157/191‑1 (MQ=255)                                                                         ttATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCgcgtgcgt                                                                                                                                                                                                                                   >  1:421157/1‑191 (MQ=255)                                                                               tGGCCGACAAAGTGACCGTGCGTACTCGCGCGCCAGGCGAAAAACCAGAAAATCGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACgcgaagcgaa                                                                                                                                                                  >  1:28733/1‑250 (MQ=255)                                                                                  ccGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGAc                                                                                                                                                              <  1:131573/251‑1 (MQ=255)                                                                                    gacaacgtgaccccgctaaCTCGCGCGGCAGGCGAAAAACCCGAGAATGGCGTCTTATGGGCATCGGCTGGCGAAGGAGAATACACACATT‑CCGACATCACCCAAGAAGATCGTGGTACTGAAATCACCCTGCCTTTGCGTATAGGCGAAGACGAGTTACTCGATGAATTGCCCTTGCGTTCCCTCAACCGCAAACACTCCGACCATATAGCGCTGCCGGTAGAGGGCGAAAAACGCGCAGAGAAAGACg                                                                                                                                                             >  2:266003/19‑250 (MQ=255)                                                                                                gTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGTTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTAt                                                                                                                                                <  2:201177/251‑1 (MQ=255)                                                                                                                 aGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGACGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGGTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGAGAAAACGTTATCTCCTGGGAGaaactc                                                                                                                                >  2:359107/1‑247 (MQ=255)                                                                                                                 aGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATc                                                                                                                                >  1:244332/1‑250 (MQ=255)                                                                                                                         aaCCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAgcg                                                                                                                       >  1:490870/1‑251 (MQ=255)                                                                                                                                       gTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGAc                                                                                                         >  2:327393/1‑251 (MQ=255)                                                                                                                                              gggAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTaa                                                                                                   <  1:120558/250‑1 (MQ=255)                                                                                                                                                         ggCGAAGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTGATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAAt                                                                                       <  2:490870/251‑1 (MQ=255)                                                                                                                                                             aaGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCAcc                                                                                   >  2:108079/1‑251 (MQ=255)                                                                                                                                                              aGGTGAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGAAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCg                                                                                  >  2:347552/1‑251 (MQ=255)                                                                                                                                                                  gAATACAC‑CGTTGCCGCCATCATCAAAGAAGATCGTGGTACTGACATCAACCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGCGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACt                                                                                                        >  2:269562/1‑225 (MQ=255)                                                                                                                                                                  gAATACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACAGCGAAACCGTTATCTCCTGGGAGAAAATCAGCAAAGCGCAGGCGCTGTGGACt                                                                                                        <  1:269562/225‑1 (MQ=255)                                                                                                                                                                    attcac‑CGTTGCCGACATCAGCAAAGAAGATCGGGGTACTGAAATCACCCTGCATCTGCGTGAAGGCAAAGACGAGTTTCTCGATGACTGGCGCCTGCGTTCCATCATCAGCAAATACTCCGACCATAGCGCGCTGCCGGTATATATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTAGCTCCTGGGAGAAAAGCAACAAAGCGCAGGCGCTGTTGTCTCTTAACAAGTCGGAAATCACCGATGAag                                                                            <  2:462394/248‑1 (MQ=255)                                                                                                                                                                     tACAC‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGa                                                                              >  2:114839/1‑248 (MQ=255)                                                                                                                                                                      acac‑CGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAag                                                                            >  1:3104/1‑249 (MQ=255)                                                                                                                                                                             ttGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCg                                                                                                                                                                                            <  1:158360/131‑1 (MQ=255)                                                                                                                                                                             ttGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCg                                                                                                                                                                                            >  2:158360/1‑131 (MQ=255)                                                                                                                                                                                    catcatCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCt                                                             <  1:65082/251‑1 (MQ=255)                                                                                                                                                                                      tcatcAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACAGATGAAGAGTACAAAGAGTTCTAc                                                           >  2:228145/5‑251 (MQ=255)                                                                                                                                                                                            aagaagATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAacac                                                     >  1:18425/1‑251 (MQ=255)                                                                                                                                                                                             agaagaTCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAAcaca                                                    <  1:242512/251‑1 (MQ=255)                                                                                                                                                                                                  aTCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGCTGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAATTCGGAAATCACCGATGAAGAGGACAAag                                                                   >  2:466245/1‑231 (MQ=255)                                                                                                                                                                                                  aTCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCTCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAag                                                                   <  1:466245/231‑1 (MQ=255)                                                                                                                                                                                                    cGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCa                                             <  1:391683/251‑1 (MQ=255)                                                                                                                                                                                                       ggTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACg                                           <  1:7403/250‑1 (MQ=255)                                                                                                                                                                                                               aaTCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATca                                                                                                                               <  2:221546/158‑1 (MQ=255)                                                                                                                                                                                                               aaTCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATca                                                                                                                               >  1:221546/1‑158 (MQ=255)                                                                                                                                                                                                                   aCCCTGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATcc                              >  1:148067/1‑251 (MQ=255)                                                                                                                                                                                                                    cccNGCATCTGCGTGAAGGCGAAGACGAGTTTCTCGATGCCTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCg                             <  2:452075/251‑1 (MQ=255)                                                                                                                                                                                                                      cTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCt                           >  1:63520/1‑251 (MQ=255)                                                                                                                                                                                                                         cATCTGCGTGAAGGCGAAGACGAGGTTCTCGATGACTGGCGCGTGCGTTCCATAATCAGCAAATACTCCGACCATATCGCGCTGCCGGAAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACACGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGAc                        >  1:302276/1‑251 (MQ=255)                                                                                                                                                                                                                             tGCGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTATAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTgg                    >  2:387649/1‑251 (MQ=255)                                                                                                                                                                                                                               cGTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGAACCGCTGACCTgg                    >  2:233703/1‑249 (MQ=255)                                                                                                                                                                                                                                gTGAAGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGc                 >  1:381814/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaGGCGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCcac              >  1:400964/1‑251 (MQ=255)                                                                                                                                                                                                                                       cGAAGACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAAcc          <  1:2666/251‑1 (MQ=255)                                                                                                                                                                                                                                           gACGAGTTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGTAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAACCgtgt      <  2:454370/251‑1 (MQ=255)                                                                                                                                                                                                                                               aGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAACCGTGTTGaa  <  2:165/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       AATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACAC‑CGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAGTTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAACCGTGTTGAA  >  NC_000913/495439‑495926

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