Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,457,177	0	T	G	26.8%	95.5 / 10.4	57	G661G (GGA→GGC)	fadJ	3‑hydroxyacyl‑CoA dehydrogenase FadJ
Reads supporting (aligned to +/- strand):  ref base T (23/18);  new base G (0/15);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTT  >  NC_000913/2456885‑2457437                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                     gCATTGTGCCTTATCCATTCACCGAATGAACCGCTGTTTTTTACCTCAATGGCGTAGGTATAGTCAATATGTGTTTATATTTTGCCCAACACTCCGCCATTCAGCGCGGAGTCATATAGCTTTGATCTTCTTATTGCAGGTCAGTTGCAGTCGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGCCCCCACGCTCGCCAGGGGTAAAAAGGGAACAATACTGCGCGGCCGGTCGTTGCTTTATTGCAACCACTTAGCCTGCGAAGAAAGAATCCGCATAGCGGA‑ACGGTCAACCgg                                                                                                                                                                                                                                                               >  2:89036/1‑299 (MQ=255) gCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGAACAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCACGCGCCGAGAGAATCGATCTAGCGGA‑ACGGTCCACCga                                                                                                                                                                                                                                                               >  1:10454/1‑300 (MQ=255)   aTTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCgaga                                                                                                                                                                                                                                                             <  1:388487/300‑1 (MQ=255)                    cACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTANAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAAt                                                                                                                                                                                                                                            >  1:568132/1‑300 (MQ=255)                           tGAAACGCTTTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATNTTGACCAACACTCCGCCATTCAGCGCTGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGNGAGAATCGATATAGCGGA‑ACGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAAt                                                                                                                                                                                                                                     <  2:725054/300‑1 (MQ=255)                           tGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATAGAGCGGA‑GCGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAAt                                                                                                                                                                                                                                     <  2:455842/300‑1 (MQ=255)                             aaaCGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAANGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGGCCANCGAGAAATGGCGGAAAACCAATGCCAAATAc                                                                                                                                                                                                                                   <  2:558083/300‑1 (MQ=255)                             aaaCGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATAc                                                                                                                                                                                                                                   <  1:116975/300‑1 (MQ=255)                                       ttttACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAata                                                                                                                                                                                                                         <  1:277349/300‑1 (MQ=255)                                               cAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTc                                                                                                                                                                                                                 <  1:731551/300‑1 (MQ=255)                                                          gTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCCTATAGCTCTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCACACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCACCGAGAGAATCGATATAGCGGC‑CCGGTCCACCGAGAAATGGCGGAGAAACAATGTCAAACACCACGCCAATATccacgccaaaaaagctcc                                                                                                                                                                                                      >  2:369678/1‑284 (MQ=255)                                                          gTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTAc                                                                                                                                                                                                      >  1:776634/1‑300 (MQ=255)                                                             taGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCNCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGa                                                                                                                                                                                                   <  1:570771/300‑1 (MQ=255)                                                             taGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACNTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGGTATGGCGGA‑GCGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGa                                                                                                                                                                                                   <  2:687174/300‑1 (MQ=255)                                                             taGTCAATAGGCGTTTACATTTTGACCAACACTCCTCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTACTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGGTAGAGCGGA‑ATGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGa                                                                                                                                                                                                   <  2:439344/300‑1 (MQ=255)                                                               gTCAATCTGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAACACGGGAACCATACTGCGTGGCAAGTCNTTGCATTATTGCAACCACTTCGCCCGCGCCTAGAGAATCGATATAGCGGA‑ACGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGNGCCAATATCCCCGTCACGCACGCTACGGATa                                                                                                                                                                                                 <  2:542721/300‑1 (MQ=255)                                                               gTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATa                                                                                                                                                                                                 >  1:465732/1‑300 (MQ=255)                                                                  aatgtGCGTGTACATTTCGACCAACACTCCGCTATTCAGCGCGGATTCATATAGCTTTGACCTTCTTACTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAATGGTTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGACAGAATCGATAGAGCGAA‑GCGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAcc                                                                                                                                                                                              <  2:24996/296‑1 (MQ=255)                                                                  aaTATGCGTCTACATTTTAACCAACACTCCGCCGTTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCNCGCAAGGGGTAAAACGGGAACCATACTGCGTTGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAAAGAATCGATAGCGATGA‑ACGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAAcc                                                                                                                                                                                              <  2:177985/300‑1 (MQ=255)                                                                                    gACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTTCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCACCTCGACCAAACACTCGCAAGGGGTAAAACGTNAACCATACTGCGTTGCAAGTCCTTGCATTACTGCAACCACTTCGCCCGCCCCTAGAGACTCGCTACCGCGGC‑ACCGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGt                                                                                                                                                                            <  2:194185/300‑1 (MQ=255)                                                                                    gACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATATCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGt                                                                                                                                                                            <  1:226238/300‑1 (MQ=255)                                                                                    gACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGt                                                                                                                                                                            >  1:470547/1‑300 (MQ=255)                                                                                      ccAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCACCTCGACCAAACGCTCGCAATGGGTAAAACGGGANCCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCTAGAGAATCGATAGAGCGAA‑ACGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTAc                                                                                                                                                                          <  2:190350/300‑1 (MQ=255)                                                                                           acTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCtt                                                                                                                                                                     <  1:688845/300‑1 (MQ=255)                                                                                              ccGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGC‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCAGTCACGCACGCTAAGGATAACCTGCTCATCAACACAAAGTACTGCttcat                                                                                                                                                                  >  1:216442/1‑300 (MQ=255)                                                                                                gCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttcattc                                                                                                                                                                <  1:446066/300‑1 (MQ=255)                                                                                                gCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCttcattc                                                                                                                                                                <  2:708676/300‑1 (MQ=255)                                                                                                   aTTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGc                                                                                                                                                             >  1:89854/1‑300 (MQ=255)                                                                                                                  atatAGCTTTGACCTTCTTATTGCAGTTCAGTTGCAGGTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCTAGAGACTCGGTATAGCTGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCacacac                                                                                                                                              <  2:749335/300‑1 (MQ=255)                                                                                                                  atatAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTTGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCTCCTAGATAATCGATATAGCTGA‑ACTGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCacacac                                                                                                                                              <  2:480328/300‑1 (MQ=255)                                                                                                                  atatAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCacacac                                                                                                                                              <  1:63178/300‑1 (MQ=255)                                                                                                                  atatAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACACGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGATAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCacacac                                                                                                                                              <  2:133892/300‑1 (MQ=255)                                                                                                                  acataGCTTTTACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCTACCAAACGCTCGCAATGGTTAAAACGGGAACCATACTGCTTTGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCTAAAGAATCGATACAGCGGA‑ACCGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCacacac                                                                                                                                              <  2:47879/298‑1 (MQ=255)                                                                                                                       gCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTc                                                                                                                                         >  1:586199/1‑300 (MQ=255)                                                                                                                          ttGACCTTCTTATTGCAGTTCAGTTGCAGTTGTCTTCCAAAAACTGTCCCCACGCGCGCCCACCGCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTTGCAATTCGTTGCATTATTGCAACCACTTCGCCCGCTCCTAGAGACTCTGTATAGATGAGTCGG‑CCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGc                                                                                                                                      <  2:400302/300‑1 (MQ=255)                                                                                                                          ttGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGc                                                                                                                                      <  1:675333/300‑1 (MQ=255)                                                                                                                            gACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCaa                                                                                                                                    <  1:233574/300‑1 (MQ=255)                                                                                                                                 tcttATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCtgc                                                                                                                               >  1:365546/1‑300 (MQ=255)                                                                                                                                   ttATTGNAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATNTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCtgcgg                                                                                                                             <  2:635091/300‑1 (MQ=255)                                                                                                                                             tCAGTTNCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAANGCTCGCAAGGGGTAAAACGTGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCTAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGAtt                                                                                                                   <  2:639823/300‑1 (MQ=255)                                                                                                                                                                aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACNGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCTAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGc                                                                                                <  2:579495/300‑1 (MQ=255)                                                                                                                                                                   aaCTTTCCTCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCTTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCTGA‑ACGGTCCACCTAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCaa                                                                                             <  2:39957/300‑1 (MQ=255)                                                                                                                                                                     cTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTTGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATc                                                                                           <  2:738467/300‑1 (MQ=255)                                                                                                                                                                                  cgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATgg                                                                              >  1:461000/1‑300 (MQ=255)                                                                                                                                                                                       aTCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCACGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGa                                                                         >  1:287424/1‑300 (MQ=255)                                                                                                                                                                                       aTCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCACGTCAAGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCATCGGGTAAATGGCGGGa                                                                         >  1:288627/1‑300 (MQ=255)                                                                                                                                                                                             aCCAAACGCTCGCAAGGGGTAAAACGTGAACCATACTGCTTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ATGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGAcc                                                                   <  2:458828/300‑1 (MQ=255)                                                                                                                                                                                                 aaCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCCATCAGCGGGTAAATGGCGGGATCGACCTGtt                                                               >  1:147360/1‑300 (MQ=255)                                                                                                                                                                                                       cGCAGGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCTAGAGAATCGATATAGCGGA‑ACGGGCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGc                                                         <  2:27726/300‑1 (MQ=255)                                                                                                                                                                                                            ggggTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTa                                                    >  1:670076/1‑300 (MQ=255)                                                                                                                                                                                                                        gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTTCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGa                                        >  1:67781/1‑300 (MQ=255)                                                                                                                                                                                                                        gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGa                                        >  1:140554/1‑300 (MQ=255)                                                                                                                                                                                                                        gAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAAACTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGa                                        >  1:193501/1‑300 (MQ=255)                                                                                                                                                                                                                        gAACCATACTGCGTGACAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTTCTTAATTCAGCATCAACATCACACACAGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGAGGGATCGACTTGTTTTTTGCTTTTACGCCCTTTCTGa                                        >  1:59199/1‑300 (MQ=255)                                                                                                                                                                                                                             aTACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACANACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGACCATa                                   >  1:561413/1‑300 (MQ=255)                                                                                                                                                                                                                                                      attGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGACCATAAAGATAGAAACCCCGGCCAtttttt          >  1:643532/1‑300 (MQ=255)                                                                                                                                                                                                                                                              cACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCttt  >  1:330617/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                     GCATTGTGCCTTATCCATTCACCGAATGAAACGCTGTTTTTTACCTCAATGGCGTAAGTATAGTCAATATGCGTTTACATTTTGACCAACACTCCGCCATTCAGCGCGGATTCATATAGCTTTGACCTTCTTATTGCAGGTCAGTTGCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGA‑ACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCCCCGTCACGCACGCTACGGATAACCTGCTCATCAACACAACGTACTGCTTCATTCAGCATCAACATCACACACCGTTCAGCAACCTGCGGTGCGGAGATTCGCCCCTGCCCTTGTGTGCCAATCAGCGGGTAAATGGCGGGATCGACCTGTTTTTTGCTTTTACGCCCTTTCTGACCATAAAGATAGAAACCCCGGCCATTTTTTCTGCCTTT  >  NC_000913/2456885‑2457437

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