Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,459	C→T	V339I (GTT→ATT)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,459	0	C	T	93.8%	150.3 / NA	48	V339I (GTT→ATT)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/22);  minor base G (0/3);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGGACTGTATCGGCAGCGTGTCAGCGCCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCT  >  NC_000913/2446217‑2446705                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        ttCGGACTGTATCGGCAGCGTGTCAGCGCCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGcc                                                                                                                                                                                                                                                  >  1:121685/1‑250 (MQ=255)            gcGGCAGCGTGTCAGCGCCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTGGCCACGGCGTCGCTGCCAGGCGGGCGCGACGGGCAAGCAGAGCCATCAGCGCAATCGCGGTTTTAGTCATTTTTTACCAGCTTGGAATATCAGTCGTCACATCGGCATTTTGCGCCCGTTGCCGTAGCATGTGATCCATTAAAATGATCGCCAGCATCGCTTCt                                                                                                                                                                                                                                      <  1:175467/250‑1 (MQ=255)                  gCGTGTCAGCGCCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATc                                                                                                                                                                                                                                <  2:151762/251‑1 (MQ=255)                    gtgtCAGCGCCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCgg                                                                                                                                                                                                                              <  1:76297/251‑1 (MQ=255)                           cgcCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTgcg                                                                                                                                                                                                                       >  2:64114/1‑251 (MQ=255)                           cgcCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTgcg                                                                                                                                                                                                                       <  2:304367/251‑1 (MQ=255)                               gACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCgg                                                                                                                                                                                                                    >  1:307979/1‑250 (MQ=255)                               gACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCgg                                                                                                                                                                                                                    >  2:336931/1‑250 (MQ=255)                                                aaaaaCTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGGGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGAGCCAGTAAAAGGATCGCCAGCATCGGTTCTGCTATCGGCACTGCGCGGATCCCGACACAGGGATCg                                                                                                                                                                                                  <  1:214093/251‑1 (MQ=255)                                                 aaaaCTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGt                                                                                                                                                                                                 <  1:253150/251‑1 (MQ=255)                                                              gattgtgagttaaacGGCACAGGTTGGGTTATTTTTTGCCCCGGCGTCGCTGCCAGGCTGGAGCTACTGGCAAGCAGCGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGCCATTTTGGGCCCGTTGCCGTAACCGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGCTCGCGACACAGGGATCGTGACAGCCTTTGGt                                                                                                                                                                                    >  2:19109/15‑251 (MQ=255)                                                                        tacacggaaaaGGTTGGGTTATTTTTTGCCACGGTGTCGATGACAGGCTGGCGCTACTGGCAAGCAGAGCCAGAAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGACTATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCCTTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACCCAGGGCTCGTGCCGGCCTTTGGTGatcatatc                                                                                                                                                                           >  2:244866/11‑247 (MQ=255)                                                                         aCCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCa                                                                                                                                                                          <  1:203141/250‑1 (MQ=255)                                                                              gCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAActtct                                                                                                                                                                    <  1:70072/251‑1 (MQ=255)                                                                                acaGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACttcttc                                                                                                                                                                  >  1:83981/1‑251 (MQ=255)                                                                                   ggTTGGGTTTTTTGTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGTAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCATCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGCTCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGcc                                                                                                                                                               <  2:241035/251‑1 (MQ=255)                                                                                       gggTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAg                                                                                                                                                           <  2:34140/251‑1 (MQ=255)                                                                                                tttGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACGTCTTCGCCTAAGCGGTTAAt                                                                                                                                                   <  2:319023/250‑1 (MQ=255)                                                                                                 ttGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATg                                                                                                                                                  <  2:293167/250‑1 (MQ=255)                                                                                                   gCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTa                                                                                                                                               <  2:307979/251‑1 (MQ=255)                                                                                                    ccACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAAt                                                                                                                                                   >  1:68053/1‑246 (MQ=255)                                                                                                    ccACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAAt                                                                                                                                                   <  2:68053/246‑1 (MQ=255)                                                                                                           gTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACGGCGCGGATCCCGACACAGGGATCGTGACGGTCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCgg                                                                                                                                       <  1:373595/251‑1 (MQ=255)                                                                                                             cGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCa                                                                                                                                     <  2:172784/251‑1 (MQ=255)                                                                                                                 gCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACgg                                                                                                                                  <  2:309431/250‑1 (MQ=255)                                                                                                                 gCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACg                                                                                                                                   <  2:43338/249‑1 (MQ=255)                                                                                                                                   gCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGttt                                                                                                               <  1:159576/251‑1 (MQ=255)                                                                                                                                    cAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTc                                                                                                              >  1:161504/1‑251 (MQ=255)                                                                                                                                        cAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGAACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAgcg                                                                                                          >  1:269055/1‑251 (MQ=255)                                                                                                                                                      cAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATg                                                                                            >  2:264274/1‑251 (MQ=255)                                                                                                                                                        aTCGCGGTCTTATTCATTTTTCACTAGCATGGCATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTACCAGGTGCTACATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACACGGATCGTGACGCCCTTTGGTGATCATCCCAACTTCTCCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGCGGTCGCTTTCCGCCCCATATGGGCACTGAt                                                                                          >  2:359778/1‑251 (MQ=255)                                                                                                                                                          cgctgTTTTATTCATTTTTTGCCAGCGTTGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCACCGCTTCTGCGATCGGCACTGGGCGGATCCCGACACAGGGATCGTTACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCATCGCCATATTGGCAATGAttt                                                                                        <  2:39026/247‑1 (MQ=255)                                                                                                                                                                  tATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGCGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCg                                                                                >  2:136335/1‑251 (MQ=255)                                                                                                                                                                     tCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACTCAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAACGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCg‑ctg                                                                            >  2:264666/1‑251 (MQ=255)                                                                                                                                                                            ttACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTAGCCAACGCGGTTAATGGTACGACCCGGCAAGGTAATGCTGTAGGTCGGTTTCAGCGCCATCTGAGCACTGATTTGCTGCCCG‑CTGCTGATAc                                                                     >  2:362793/1‑251 (MQ=255)                                                                                                                                                                             taccctcttttcctatCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCACATTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTTTCTTACGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATAcc                                                                    >  1:188176/14‑251 (MQ=255)                                                                                                                                                                                               gTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGccgcccgc                                                  >  1:373852/1‑251 (MQ=255)                                                                                                                                                                                                  ttCACATCGGAATACTGCGCCCGTTGCCGTAACAGGTGCTCCATTAACATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATg                                               >  2:55944/1‑251 (MQ=255)                                                                                                                                                                                                     acaTCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACACGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATCTGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATGGtt                                            >  1:244998/1‑251 (MQ=255)                                                                                                                                                                                                         cGGCATTTTGCTTCCGTTTCCGTAACAGGTGATCCATTAAAATGATCGCCAGCACCGCTTCTGCGTTCGGCTCTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGTTCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTTGTGGTCGGTTTCAGCGCCATAGGGGTAATGATTTTCTGCCCG‑CGGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGct                                         <  1:177555/250‑1 (MQ=255)                                                                                                                                                                                                           gCATTTTGCGCCCGTTGCCGTAACGGGTGTTCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTCCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATGTTTGCTCTg                                      <  1:262891/251‑1 (MQ=255)                                                                                                                                                                                                                 tGCGCCCGTTGCCGTAACAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAAcc                                >  1:261670/1‑251 (MQ=255)                                                                                                                                                                                                                     cccGTTGCCGTAACAGGTGAGCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCt                            <  1:372889/251‑1 (MQ=255)                                                                                                                                                                                                                              gTAACATGTGATCCATTAAAAGGATCGCCAGCATCGCTTCTGCGATCGGCAATGCGGTTATCCCGACACAGGGATCGTTACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCTGTTAATTGTACGACCCGGCACTGTAATGCTGGAGGTCGTTTTCAACGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTTATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTAGAAACCGTCTTTGGTGAtt                   <  1:291125/251‑1 (MQ=255)                                                                                                                                                                                                                                aaCAGGTGATCCATTAAAATGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAACCCGTCTTTGGTGATTTc                 >  1:163403/1‑251 (MQ=255)                                                                                                                                                                                                                                       gATCCATTAAAATGATCGCCAGCATAGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATAATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCACCATATGAGCAATGATTTGCTGCTCG‑CTGCTGATACCGCCGAGACTGCCGCCCGACTGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCgg          >  2:100598/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTCCATTACAATGATCGCCAGCACAGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGCCCTTTGGTGATCATCTCAACTTCCTCGCCAACGCGGTTAATGGTACGACCCGGCACGGCAATGCTGCCGGCCGGTTTCAGCGCCATATGCGCACTGATTTGCTGCCCG‑CTGCTGATACCGCCGCGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTCTGGTGATTTCCTCGCGGt         >  2:231942/1‑251 (MQ=255)                                                                                                                                                                                                                                               aaaaGGATCGCCAGAATCGCTGCTGCGTTCGGCACGGCGCGGACCCCGACACAGGGATCGTGACTGCCTTTGGTGATCATCGCAACTTCTTCGCCAAAGCGGTTAAGGTTACGACCCGGCACGGGATTGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTGGC‑GCCCGGTTGTGGTTACCGCCGATAATGCCGCCCGCATGGTTGCTCTGGAAGCCGTCTTTGGTGATTTCATCGCGGTTCTGGCt  <  2:166592/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        TTCGGACTGTATCGGCAGCGTGTCAGCGCCGACAATACAGCCATTAGAAAAACTGCCTATCGATTGTGCGCTACCCGGCACAGGTTGGGTTATTTTTTGCCACGGCGTCGCTGCCAGGCTGGCGCTACTGGCAAGCAGAGCCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCG‑CTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCT  >  NC_000913/2446217‑2446705

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