Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,825,558	G→A	G117S (GGC→AGC)	recG →	ATP‑dependent DNA helicase RecG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,825,558	0	G	A	98.0%	170.7 / NA	51	G117S (GGC→AGC)	recG	ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/29);  minor base C (1/0);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATGGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTCACAAGGAATGATGACGCTACCGGAAGCGTT  >  NC_000913/3825317‑3825802                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       cctctTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCaa                                                                                                                                                                                                                                               <  1:170850/250‑1 (MQ=255)       acacCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGcgcg                                                                                                                                                                                                                                         <  2:32657/251‑1 (MQ=255)                     gcTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGAGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATAAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCg                                                                                                                                                                                                                          >  1:312711/1‑251 (MQ=255)                                              ccccATCGGAGAACTACGGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTTCTGAACTGCAATATCTCTTTCGGCGGTCGGCGGATGATGTCCTGCCAGATCAGCGACGGTTCCGGCATTCTCACCATGCGCTTTTTCAATTTCAGCTCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCTAAGCAAAGCGCGGT‑AAAGATGGTGCGGAGATGATCCACCCGTAATACcgcg                                                                                                                                                                                                 <  2:46832/251‑1 (MQ=255)                                               cccATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGACCTGCACTATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCCGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCaacg                                                                                                                                                                                                                                             >  2:44875/1‑205 (MQ=255)                                               cccATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAg                                                                                                                                                                                                                                             <  1:44875/207‑1 (MQ=255)                                                 cagcGGAGATCTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTTAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCATCGACGGTTCCGGCTTCCTCACCATGCGCTTTTTCAATTTCTGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGc                                                                                                                                                                                              <  2:230126/248‑1 (MQ=255)                                                               tGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGGTCTCAGc                                                                                                                                                                                >  2:94006/1‑251 (MQ=255)                                                                  cGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATATCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGGTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGACAAATCGCCTGGCGGCGGGCCGCCGTGTCCTTGCTTATAGCGAAGCAAAGCGCGGG‑AAATATGGTGCGGGGATGATCCACCCGGAATACCGCGTTCAGGGGTATCTCAGCACg                                                                                                                                                                             >  2:231084/1‑251 (MQ=255)                                                                     gCGTTTATGCCACGGTGGAAGGCGAATTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAAGTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCa                                                                                                                                                                          <  2:123672/251‑1 (MQ=255)                                                                       gTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGGGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGcc                                                                                                                                                                           >  1:228629/1‑248 (MQ=255)                                                                             gCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTAcc                                                                                                                                                                  >  2:245218/1‑250 (MQ=255)                                                                               catggTGGAAGGCGAAGTGCTGAACTGCAATAGCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGTTCAGCGACGGGTCCGGCATTCTCACCATGCGCTTTTGCAATTTCAGCTCGGCAATGAAAAATTGCCTGGCGGCGGTCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAAGATGGTGCTGAGATGATCCAACCGGAATACCGCGTGCAGTGCAATCTCAGCACGCCAGAATTACAgg                                                                                                                                                                <  2:78896/248‑1 (MQ=255)                                                                                  ggTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGaaa                                                                                                                                                             <  1:94006/251‑1 (MQ=255)                                                                                   gTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGAAGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGAGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATAACGCGTGCAGGGCGCTCTCAGCACGCCAGAATTACAGGaaa                                                                                                                                                             >  1:110591/1‑250 (MQ=255)                                                                                        aGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAGCGCTCa                                                                                                                                                       <  1:86058/251‑1 (MQ=255)                                                                                            gaggTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCGCCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGcc                                                                                                                                                   <  2:225346/248‑1 (MQ=255)                                                                                                 gCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCGGCCAGATCAGCGACGTGTCCGGCATCCTAACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGtt                                                                                                                                               <  2:71849/250‑1 (MQ=255)                                                                                                             tatCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCCGCGACGGCTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCACTGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACgg                                                                                                                                  >  1:265237/1‑251 (MQ=255)                                                                                                                ctcCTTCGGCGGTCGGAGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACAGCGTGCAGGGCGATCTCAGCACGTCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGa                                                                                                                                 >  1:22412/1‑249 (MQ=255)                                                                                                                  ccTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGc                                                                                                                             <  2:252525/251‑1 (MQ=255)                                                                                                                       ggcggtcggcggATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTaa                                                                                                                         >  1:150726/1‑250 (MQ=255)                                                                                                                                cggATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCAc                                                                                                               <  1:231287/251‑1 (MQ=255)                                                                                                                                 ggATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACg                                                                                                              >  1:168756/1‑251 (MQ=255)                                                                                                                                 ggATGATGACCTGCCAGAGCAGCGTCGGTTCCGTAATCCTCACCATGCGTTTGTTCAATTTCAGCGCGGAAATGAAAAATAGCCTGGCGGCGGGCCGCTGTGTACGGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGTTGGGGAGATGTTCCACCCTGAAGACCGCGGGCAGGGCGATCTCAGCACGCCAGAAGTACGGGAAACGTTCACGTGGGTTGAACCAACAACGGAAGGCGTAAAGCAGGCCACg                                                                                                              <  1:294636/251‑1 (MQ=255)                                                                                                                                     gatgaCCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGAGCGGAGATGATCCACCCGGAATACCGCGTGCCGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCCCGCTGc                                                                                                          >  2:196405/1‑251 (MQ=255)                                                                                                                                                agaacAGCGACGGTTCAGGCATCCTCACCATGCGCATTTTCAATTGCAGCGCGGCAAGGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGAGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAAcc                                                                                               <  2:52289/247‑1 (MQ=255)                                                                                                                                                      gCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑g                                                                                        <  1:29004/251‑1 (MQ=255)                                                                                                                                                      gCGACGGGTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCGGCGTAAAGTAACCGACCA‑g                                                                                        <  2:150726/251‑1 (MQ=255)                                                                                                                                                              tgcGGCAACCACACCATGCGCTTTTTCAAATTCAGCGCGGTAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCTCGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTgg                                                                                 <  2:3415/248‑1 (MQ=255)                                                                                                                                                                 ggCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGAGGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCt                                                                             <  2:245214/251‑1 (MQ=255)                                                                                                                                                                      ccTCCCCATGTGCCTTTTCAATTGCATCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTTGACTTGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCAGAGATGATCCCCCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTTCGTAAATGAACCGACCA‑TGCGCTGGATCTGCTCg                                                                        <  2:251100/251‑1 (MQ=255)                                                                                                                                                                            cATGCGCTTTTTCAATTTCAGCACGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGATTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACAcc                                                                   >  2:116203/1‑250 (MQ=255)                                                                                                                                                                                          aTTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCAGGTTTATCCAACAACGGAAGGCGAAAAGAAAGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCACCATTgacga                                                     >  1:98539/1‑247 (MQ=255)                                                                                                                                                                                                  gcgGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTTCTCGACACCTGCTCCGTTGAAGAACTCCTgcc                                            <  2:177571/251‑1 (MQ=255)                                                                                                                                                                                                  gcgGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTgcc                                            >  1:119060/1‑251 (MQ=255)                                                                                                                                                                                                       aaTGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATCGCGAAGCAAAGCGCGGG‑AAATATGGCGCGGAGATGATCCACCCGGAATACCACGTGCAGGGCGATCTAAGCACGCCAGAATTCCAGGCAACGCTCACGGCGGATTATCCAACAACGGCAGGCGCAAAGCAGGACAAGCAGCGACAATTAAACG‑CCACGACGATGAGTCTGCGCTACCCCTGCGACATTGAAGAACTCCTgccgcc                                         >  2:172458/1‑249 (MQ=255)                                                                                                                                                                                                           aaaaaTAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACt                                   <  1:194483/251‑1 (MQ=255)                                                                                                                                                                                                            aaaaTAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTg                                  >  1:233068/1‑251 (MQ=255)                                                                                                                                                                                                               aTAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTACAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGCATACCGCGTGCAGGGCGATCTCACCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTc                                >  2:132646/1‑250 (MQ=255)                                                                                                                                                                                                                 aGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTATAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCGGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTcaca                             <  2:299074/251‑1 (MQ=255)                                                                                                                                                                                                                 aGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGG‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGCAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGCACTGTcac                              >  2:212781/1‑250 (MQ=255)                                                                                                                                                                                                                   ccTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTcaca                             >  1:83102/1‑249 (MQ=255)                                                                                                                                                                                                                    cTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTca                               <  1:50661/246‑1 (MQ=255)                                                                                                                                                                                                                    cTGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTca                               >  2:50661/1‑246 (MQ=255)                                                                                                                                                                                                                     tGGCGGCGGGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTCACAAgg                          >  1:120783/1‑250 (MQ=255)                                                                                                                                                                                                                        cggcggGCCGCCGTGTACTGGCTTATAGCGAAGCAAAGGGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTCACAAGGAatg                      <  2:228629/251‑1 (MQ=255)                                                                                                                                                                                                                            gggCCGCCGTGTACTGGCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTCACAAGGAAtgatga                  <  2:119060/251‑1 (MQ=255)                                                                                                                                                                                                                                    gtgtACTGGCTTATATCGAAGAAAAGCG‑GGTGAAATATGGTGCGGACAGGATCCTCGCGGAAGACCGAGTGCAGGGCTAACTCAGCACGCCAGAAGTACAGGAAACGCTCACGCCGGGTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGTGTAAATTAACCGACCA‑GGGGCTGGGTCGGCTCGACACCGGCGCCAATGAAGAACTTCATCCGCCGGAACTGTCGCAAGGAAATATGACGTTACCg          <  2:13032/251‑1 (MQ=255)                                                                                                                                                                                                                                          gggCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTCACAAGGAATGATGACGCTACCGGAAGCg    <  1:245218/250‑1 (MQ=255)                                                                                                                                                                                                                                            gCTTATAGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAAGACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTCACAAGGAATGATGACGCTACCGGAAGCGtt  <  2:234288/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       ACTCTTACACCTTCCCCTGCGCTACGAAGATCGCACCCATCTCTACCCCATCGGAGAACTACTGCCGGGCGTTTATGCCACGGTGGAAGGCGAAGTGCTGAACTGCAATATCTCCTTCGGCGGTCGGCGGATGATGACCTGCCAGATCAGCGACGGTTCCGGCATCCTCACCATGCGCTTTTTCAATTTCAGCGCGGCAATGAAAAATAGCCTGGCGGCGGGCCGCCGTGTACTGGCTTATGGCGAAGCAAAGCGCGGT‑AAATATGGTGCGGAGATGATCCACCCGGAATACCGCGTGCAGGGCGATCTCAGCACGCCAGAATTACAGGAAACGCTCACGCCGGTTTATCCAACAACGGAAGGCGTAAAGCAGGCCACGCTGCGTAAATTAACCGACCA‑GGCGCTGGATCTGCTCGACACCTGCGCCATTGAAGAACTCCTGCCGCCGGAACTGTCACAAGGAATGATGACGCTACCGGAAGCGTT  >  NC_000913/3825317‑3825802

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