Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	109,960	G→A	R561H (CGT→CAT)	secA →	protein translocation ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	109,960	0	G	A	98.3%	192.6 / NA	61	R561H (CGT→CAT)	secA	protein translocation ATPase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (24/35);  minor base C (1/0);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATT  >  NC_000913/109715‑110201                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  gAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCAc                                                                                                                                                                                                                                               <  1:88321/250‑1 (MQ=255) gAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACg                                                                                                                                                                                                                                              <  1:207794/251‑1 (MQ=255)    cGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCGGCGGGGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAAGTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAAt                                                                                                                                                                                                                                           <  2:41932/251‑1 (MQ=255)          aTTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCcgtc                                                                                                                                                                                                                                     <  2:205167/251‑1 (MQ=255)              cACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTAt                                                                                                                                                                                                                                 <  2:472142/251‑1 (MQ=255)               aCGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCTTCGTATc                                                                                                                                                                                                                                <  1:471130/251‑1 (MQ=255)                   cAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATa                                                                                                                                                                                                                            <  2:147011/251‑1 (MQ=255)                                                                gACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGttct                                                                                                                                                                                <  1:77245/250‑1 (MQ=255)                                                                                  ggCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGGTCTTCCCGTTTCTTCCTGTCgg                                                                                                                                                             >  1:160125/1‑250 (MQ=255)                                                                                     gggTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGGACGTCACGATGCGGTACTGGAAGCCGGTGGCCTGCATATCATCGGTAACGAGCATCGCGAATCCCGTCGTATCGATAACCAGTTGAGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTACCGTTTCTACCTGTCGATgg                                                                                                                                                          >  2:309433/1‑251 (MQ=255)                                                                                      ggTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATgg                                                                                                                                                          >  1:57655/1‑250 (MQ=255)                                                                                            ggTACAGATATTGTGCTCGGTGGGAGCGGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAACTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTTCTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTGCTACCTGGCGATGGAAGATgc                                                                                                                                                   <  2:265724/251‑1 (MQ=255)                                                                                               aCAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGCAAAAATTAAAGCCGACTGGCAGGTAAGTCCCGATGCGGTACTTGAAGCAGGTGGCCTGCATATCATCGGTTCCGGGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCCTTCTGGTCGTCCGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGTTgcgtt                                                                                                                                                >  2:259554/1‑249 (MQ=255)                                                                                                   atatTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGCACGTCACGATGCGGTCCTGGAAGCAGGTGGCCTGCCTATCATCGGTACCGGGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTATGGTCGTCAGGGGGGTGCTGGTTCTTCCCGTTTCTACCTGTAGATGGAAGATGCGCTGATg                                                                                                                                            >  1:86146/1‑251 (MQ=255)                                                                                                       tgtgCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGCCCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCt                                                                                                                                                <  2:372051/243‑1 (MQ=255)                                                                                                       tgtgCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCt                                                                                                                                                >  1:372051/1‑243 (MQ=255)                                                                                                        gtgCTCGGTGGGAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTa                                                                                                                                        <  2:424434/250‑1 (MQ=255)                                                                                                             cGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTAtttt                                                                                                                                    <  1:517993/249‑1 (MQ=255)                                                                                                                    aGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACAGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGTATATTATCGGTACCGAGCATCACGAATCCCGTCGTCTCGCTAACCAGTTGCGCGGTCGGTCTGGTCGTCTGGGGGATGCGGGTTCTTCCCGTTTCTACCTGTCGATTGAAGATGCGCTGATGCGTATTTGTGCTTccga                                                                                                                           <  2:253762/251‑1 (MQ=255)                                                                                                                          caggcagAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTAGCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCTCGTTTCTACCTGTCGAGGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGt                                                                                                                     <  1:56796/251‑1 (MQ=255)                                                                                                                              cagAAGTTGCCGCGCTGGACAATCCGACTGCAGAGCAAATTGAAAAAATTACAGCCGACTGGCAGGTACGTCACGAATCGGTACTGGAAGAAGGTGGGCTGCATCTCATCGGTCACGAGCATCACTAACCCCGTCGCATACTTAAACACTTGCGCGGTCGTTATGGTCGTCAGGGGGATGCTGGTTCTTCCCGGTTCTACCTGTCGATGGAAGATGCGCTGATGAGTATTTTTGCTTccgaccgtatctcc                                                                                                                 >  2:183886/1‑244 (MQ=255)                                                                                                                               agAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCg                                                                                                                <  1:351208/251‑1 (MQ=255)                                                                                                                                gAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTAttt                                                                                                                                     <  1:354386/229‑1 (MQ=255)                                                                                                                                gAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGAATGGCAGGTACGTCACGATGCGGTACTGGAAACAGGTGGCCTGCATCTCATCGGTACCGAGCATCACGAATCCCGTCGGATCGGTAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTTTCGATGGAAGATGCGCTGATGCGTAttt                                                                                                                                     >  2:354386/1‑229 (MQ=255)                                                                                                                                          cgcTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGt                                                                                                     >  1:99629/1‑251 (MQ=255)                                                                                                                                             tGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTaa                                                                                                   >  1:178609/1‑250 (MQ=255)                                                                                                                                              ggAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTaa                                                                                                   <  2:160600/249‑1 (MQ=255)                                                                                                                                                 acctCCGACCGCAGAGCAAATTGCAACAATTAAAGCCGACTGGCAGGTACGTCACGATTCGGTACTGCAAGCAGGTGGCCTGCATATCCTACGTACCGAGCGTCACGCATCCCGTCGTATCACTAACCCGTTGCGCGGCCGTTCTGGTCGTCCGGGGGCTGCTGGTTCTTCCCCTTTCTACCTGTCTACGGCAGTTGCGCTGATGAGTATTTTTGCTTCCGACCGAGTCTCCGGCATGATGCGTaacatgg                                                                                              >  2:326952/4‑246 (MQ=255)                                                                                                                                                 aaaTCCGACCGCAGAGAAAATTGCAACCATTAAAGCCGACTGGCAGGTAAGTCACGATGCGGCAATGGAAGCACGCGGCCTCCATATCATCGGTACCGACACTCACAAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGGCGACAGGGGGATGCCGGTTCTTCCCGTTTCTACCTGTCCATGGACGATGCGCTGTTGCGTATTTTTGCTTCCGGCCGAGTATTCGGCCTGATGCggaacatgg                                                                                              >  2:334288/1‑243 (MQ=255)                                                                                                                                                   aTCCGACCGCAGAGCAAATTGAAAAAATTAACGCCGACTGGCAGGTACGTCACGATGCGGCACTGGAAGCAGGTGGCCTGCATATCATCGGGACCCAGCATCACGTATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTTTTTTTTCTTCCGACCGAGTATCCGGCCTTATGCGTAAACTgggg                                                                                            >  2:112474/1‑250 (MQ=255)                                                                                                                                                   aTCCGACAGCAGAGCAAAGTGAAAAAATTAAAGCCGACTGGCAGGTACGTAACAAGGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATAACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGGTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGATTATCCGGCATGATGCGTAAACTGGGt                                                                                            <  2:477174/251‑1 (MQ=255)                                                                                                                                                       gACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGAGGCGGTACTGGAAGCAGGGGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGCTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGGATTTTTTCTTCCGTCCGAGTATCCGGCATGATGCGTAAACTgggagtg                                                                                         >  2:123349/1‑246 (MQ=255)                                                                                                                                                                      tGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGcc                                                                          <  1:86643/250‑1 (MQ=255)                                                                                                                                                                      tGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTCGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATTGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGcc                                                                          <  1:334288/250‑1 (MQ=255)                                                                                                                                                                       gAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCAt                                                                        <  1:112474/251‑1 (MQ=255)                                                                                                                                                                        aaaaaaTTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGCTAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTACACTGGGTATGAAGCCAggagttcggc                                                                         >  2:124114/1‑241 (MQ=255)                                                                                                                                                                        aaaaaaTTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAgg                                                                                 <  1:124114/241‑1 (MQ=255)                                                                                                                                                                         aaaaaTTAAAGCCGACTGGCAGGTACGTCACGAGGCGGTACTGGAAGCAGGGGGCCTGCATATCATCGGTACCGAGCATCACGAAGCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCTAGTATCCGGCATTAGGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTg                                                                      <  1:218998/251‑1 (MQ=255)                                                                                                                                                                           aaagtAAAGCCGACTGGCAGGAACGACACGATGCGGTACTGGAAGCAAGTGGCCAGCATCTCATCGGTACCGAGCATCCCGAATCCCGCCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTAACTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCCTCCGACCGGGTATCCGGCATGCTGCGTAAACTGGGTATGGAGCAAGGCGCAGCCATTGa                                                                     >  2:384014/5‑250 (MQ=255)                                                                                                                                                                           aaaTTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTTCCTGTCGATGGAAGATGCGCTGATGCGTATTGTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGGATGAAGCCAGGCGAAGCCATTGaa                                                                    <  2:540056/251‑1 (MQ=255)                                                                                                                                                                            aaTTAAAGCCGACTGGCAGGTACGTCACGATGCTGTACTGGAAGCAGGTGGCCTGCATATCATCAGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAac                                                                   <  2:176973/251‑1 (MQ=255)                                                                                                                                                                                aaaGCCGACTGGCAGGAACGTCACGATGCGGTAATGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGGAAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGGTGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAacaacc                                                               >  2:37046/1‑248 (MQ=255)                                                                                                                                                                                     cGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTggg                                                          >  1:374125/1‑251 (MQ=255)                                                                                                                                                                                      gACTGGCAGGTGCGTCACGAGGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGt                                                         <  1:213102/251‑1 (MQ=255)                                                                                                                                                                                       aCTGGCAGGGACGTCACGATGCGGTACTGGAAGCAGGGGGCCCGCATATCATCGGTAACGAGCATCACGAATCCCGTCGTATCGATAAACAGTTGCGCGGGCGTTCTGGTCGTCCGGGGGGTGCTGGTTCTTCCCGGTTCTACCTGTCGATGGAAGATGCGCTGTTGGGTTTTTTTGCTTCCGACCGAGTATCCGGCATGATGCGGACACCGGGGATGAAGACAGGCGAAGCCATTGAACACCCGTggggg                                                        >  2:207309/1‑249 (MQ=255)                                                                                                                                                                                          ggCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCGGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGAGGCTGGTTCTTCCCGTTTCTACCTGTCGAGGGAAGATGCGCTGATGCGTATTTTTGCGTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGa                                                       <  2:376543/249‑1 (MQ=255)                                                                                                                                                                                           gCAGGTACGTGACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTTCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTGCCCGTTTCTACCTGTCGAGGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTa                                                    <  2:24200/251‑1 (MQ=255)                                                                                                                                                                                             aGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGGGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTaaa                                                  <  1:259554/251‑1 (MQ=255)                                                                                                                                                                                                     cACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTg                                           <  2:147572/250‑1 (MQ=255)                                                                                                                                                                                                        gATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGGGGGTGACTAAAGCGATTGCCa                                        >  2:423486/1‑250 (MQ=255)                                                                                                                                                                                                         aTGCGGTACTGGAAGCAGGTGGCCTGCATATCCTCGGTACCGAGCATCACGAATCCCGTCGTATCGATCACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAAc                                      >  2:59332/1‑251 (MQ=255)                                                                                                                                                                                                         aTGCGGTACTGGAACCAGGTGGCCAGCATATCATCGCTACCGAGCATCCCCAATCCCGTCGTAACGCCAACCAGATGCGCGGTCGTTCAGGTCGACAGGCAGATGCAGGTTCTTCCCGCTTCTACCTGTCGATGGAAGATGCGCTGATGCGCATTTTTCATTCCGACCGAGAATCAGGCATGATGCGTAAACTGGGAATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGCTTGCCaa                                       >  2:45781/1‑250 (MQ=255)                                                                                                                                                                                                                aCTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGc                               <  2:521343/251‑1 (MQ=255)                                                                                                                                                                                                                      aGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTaa                           >  1:158862/1‑249 (MQ=255)                                                                                                                                                                                                                      aGCAGGGGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACGAAAGCGATTGCCAACGCCCAGCGTAAAg                         <  2:269141/251‑1 (MQ=255)                                                                                                                                                                                                                         aGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTg                      <  1:59332/251‑1 (MQ=255)                                                                                                                                                                                                                            tGGCCTGCCTATCATCGGTACCGTGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGAACGTTCTGGTCGTAAGGGGGGTGCTGGTTATTCCCGTTTATAACTGTCGATGGAAGATGCGCTGATGAGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGCAGCAATGGAACTCCCGGGGGTGACTACAGCGATTCCCacggccccacgacaagttgaa                    >  2:537852/1‑230 (MQ=255)                                                                                                                                                                                                                             ggCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAg                  <  1:52419/251‑1 (MQ=255)                                                                                                                                                                                                                                      atcatcGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACtt         <  1:542143/251‑1 (MQ=255)                                                                                                                                                                                                                                       tcatcGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGAGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTc        >  2:37266/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACAtt  >  1:181401/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GAACGCCAAATTCCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATT  >  NC_000913/109715‑110201

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