Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	37,338	0	T	G	25.8%	119.4 / 10.0	66	T163P (ACC→CCC)	caiC	carnitine‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (12/37);  new base G (16/1);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 4.68e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCG  >  NC_000913/37088‑37585                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          aCCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGt                                                                                                                                                                                                                                                           <  2:122874/251‑1 (MQ=255)      caaaaaGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAg                                                                                                                                                                                                                                                      <  1:203621/248‑1 (MQ=255)         aaaGGTGGCCCCGGCAGAAAACGCCGCAATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTAAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgtt                                                                                                                                                                                                                                                   <  2:122805/251‑1 (MQ=255)          aaGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgttg                                                                                                                                                                                                                                                  <  1:244710/251‑1 (MQ=255)            ggTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgttgtt                                                                                                                                                                                                                                                <  1:118948/251‑1 (MQ=255)                 ccccGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACCAGGGGGCAGGGTGTTGAttt                                                                                                                                                                                                                                           >  1:228782/1‑251 (MQ=255)                  cccGGCATAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATGTGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAGGCGCAGGTTGTAATGGGGAATCACCACTCCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGAtttt                                                                                                                                                                                                                                          <  1:383216/251‑1 (MQ=255)                   ccGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCTGTGCATAGCACAAGGTGGCAGgtt                                                                                                                                                                                                                                                   <  1:150862/241‑1 (MQ=255)                                  ccctgcgccgcattccccTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTACCGCAGTTCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTTTATTGGGTAACCACCTCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGa                                                                                                                                                                                                                          <  1:83288/234‑1 (MQ=255)                                      tCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTCCCGGCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAAc                                                                                                                                                                                                                        <  2:154273/249‑1 (MQ=255)                                                acacTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTCGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGGAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCCCGCcatcct                                                                                                                                                                                                            >  1:304153/1‑249 (MQ=255)                                                        aaTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGa                                                                                                                                                                                                    <  1:312494/251‑1 (MQ=255)                                                             atatGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGc                                                                                                                                                                                               <  1:186871/251‑1 (MQ=255)                                                                 gAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATAACACAAGGGGGGAGGTTGTTGATTTTTTAGGTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGGGCcacc                                                                                                                                                                                           >  1:10968/1‑250 (MQ=255)                                                                   aaCGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGCGCACACTGCCAGGCGGAGTAATATCCCGCGAAGCGCAGGTTGTACTGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTCAACGAACTCCCCCCATCCTCAGCGGGGAGTGCcacctc                                                                                                                                                                                         >  1:201969/1‑248 (MQ=255)                                                                               aCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATCGCACAAGGGGGCAGGGTGTTGATTTTTCAGTTGAGGAAACGAACTCACGCCCTCCTCAGCGGGGAGTGCCACATCTGTCC‑GGGAAAt                                                                                                                                                                            >  2:361940/1‑251 (MQ=255)                                                                                 cGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAAtgt                                                                                                                                                                          <  1:352946/251‑1 (MQ=255)                                                                                    cAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGcc                                                                                                                                                                                              <  1:362900/229‑1 (MQ=255)                                                                                    cAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCCGGGTGTTGTTTTTTTAGTTGAGGAAACGAACTCACGCCCTCATCAGCGGGGAGTGcccgatcgggagagcgtcgggg                                                                                                                                                                          >  2:362900/1‑229 (MQ=255)                                                                                     aGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGGTTTTTCCGTTGAGTAAAAGAAATCACCCCATCATCAGCGGGGAGTGCCCCATCTTTCAAGGCAAA‑GTGccc                                                                                                                                                                      >  2:354339/1‑250 (MQ=255)                                                                                     aGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCCTAGCACAAGGGGGGGGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCCCGCCATCATCCGCGGGAGGTGCCCCATCTGTCA‑GGCaaagatgccg                                                                                                                                                                      >  2:96759/1‑244 (MQ=255)                                                                                             gTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTCCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAg                                                                                                                                                              <  1:334731/251‑1 (MQ=255)                                                                                                    tCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCCAAAGGTGGCCGGGTGTTGATTTTTCAGTTGAGGGAACGAACTCACGCCATCCTCAGCGGGAAATGGCCCATCTCTCA‑GGCAAATGGGCCGCCAtttggtgggatt                                                                                                                                                       >  2:407995/1‑241 (MQ=255)                                                                                                            ggcgcccggCCAGGCGGAGTAATACCCAGCGAAGCGGGGGTTGTAATGGGGGATCCCCACACTTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGGATTTCCGCCGTTTCGTCAGTCGATAGCGGCGGTGCATAGTACAAGGTGGCAGGTTGTTGATTTTTCAGTTGGGTAAACGAACTCACGCCATCATCAGCGTGAAGTGCCACATCTGTCA‑GGCAAATGGGCCTCAATTGATTGGGG‑TGTTCctgct                                                                                                                                              <  2:298248/243‑1 (MQ=255)                                                                                                                cacTGCCAGGCGGAGTAATCTCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAAt                                                                                                                                          >  1:349276/1‑251 (MQ=255)                                                                                                                cacTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGTTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAAt                                                                                                                                          <  1:312063/251‑1 (MQ=255)                                                                                                                           ggAGTAATATCCAGCGAAGCGCGGGTTGTAATGGGTAATCACCACTCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACa                                                                                                                               <  2:64952/251‑1 (MQ=255)                                                                                                                              gTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATATCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACCCCACCATCCGCGGGAAGTTGCACATCTGTCA‑GGGAAACATGCCGCAATTTCGGGGCA‑TCTTCCTTCTTAATCTTTTggttcat                                                                                                                              >  2:390675/1‑243 (MQ=255)                                                                                                                                aaTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCAACACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGAACAAGGTGGCAGGTTGTTGATTTTTTAGTTGAGTAAACGGACTCAAGCCATCATCAGCCGGAAGTGCCACATCTGTCA‑AGCCAATGTGCCGCAAATGAGTGGCA‑TCTTCCTGCTGACTCTGTTTATAAATAggg                                                                                                                          >  1:22022/1‑250 (MQ=255)                                                                                                                                  aacccAGCGAAGCGCAGGTGTTAATGGGTAACCACCACCCCTTCCGGTCGGGAGGGGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATa                                                                                                                        <  2:228782/248‑1 (MQ=255)                                                                                                                                   atCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGAt                                                                                                                         >  1:159591/1‑249 (MQ=255)                                                                                                                                     ccAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGaa                                                                                                                     <  1:375785/251‑1 (MQ=255)                                                                                                                                                   gTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGGCCAAGGGGGCAGGGTGGTGATTTTTCAGTTGAGTTAACGAACTCACGCCACCATCCGCGGGGAGGGCCACATCCTTCC‑GGCAAATGTGCCGCCATTTAGTGGGAGTCT‑CCTGCTGCAACTGTTTGTACATATGATAGAAATGCGCACCGGtac                                                                                                       >  2:396929/1‑249 (MQ=255)                                                                                                                                                        aaTGGGTAATCACCACACCTTCCGGACGGTAGGTGTGGCCGGAGTTGAAAAAAATTACCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACcagc                                                                                                  <  2:152899/251‑1 (MQ=255)                                                                                                                                                          tGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCCTAACACAAGGTGGCAGGTTGTTGCTTTTTCAGTTGAGTAAAAGAACCCACGCCCACCTCAACGGGAGGTGACACATCTGTAA‑GGCAAATGTGCCGCAACTGGGGGGCA‑TCTTCCTGCTGACGCTTGTGTTACATAGGATAGAATTgtggcctgggcccccgccg                                                                                                >  2:136673/1‑233 (MQ=255)                                                                                                                                                              tAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTa                                                                                                                                                                                                                               <  1:352485/122‑1 (MQ=255)                                                                                                                                                              tAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGGTGTTGGTTTTTCCGTTGGGTa                                                                                                                                                                                                                               >  2:352485/1‑122 (MQ=255)                                                                                                                                                               aaTCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTTAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGGGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAAATGGGCACTGGTCACCAGCAGGCACg                                                                                           >  1:175388/1‑251 (MQ=255)                                                                                                                                                                    ccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTgg                                                                                      <  1:177294/251‑1 (MQ=255)                                                                                                                                                                               ggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCCGCGGGAAGTGCCACATCTGTCC‑GGCAAAAGTGCCGCAATTTGGGGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGGTAGAATTTCGCCCTGGTCACCAGCAGGCACGCCCGGCTAtttttc                                                                             >  2:227520/1‑247 (MQ=255)                                                                                                                                                                                          ggtgCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCCTCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCCGCAGGCACGCCCGGCTATTTTGCAGGATccccccgc                                                                >  2:410166/1‑245 (MQ=255)                                                                                                                                                                                           gtgCCGGCGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGAACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCCTCATCAGCGGGGAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTTAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGGCACCAGCAGGCACGCCCGGCTTTTTTGGAGGGTCCACGCGCt                                                               >  2:169207/1‑251 (MQ=255)                                                                                                                                                                                              catgAGGTGAATATAATGTCCTCCGTCTCGTCAGTCGCGAGCTGCGTTGCATAGCACAAGGTGGCAGGTTGTTGATTTGTCAGTTGAGTAATCGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATTTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCATGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTc                                                            <  2:235993/248‑1 (MQ=255)                                                                                                                                                                                                ggAGGTGAAGAGAATTTCCGGCGTATCGTCAGTCGGTAGCGGCGGTGCATAACACAAGGTGGCAGGGTGTTGGTTTTTTAGTTGAGGAAACGAACTCCAGCCCTCATCAGCGGGACGTGGCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGGGGGA‑TCTTCCCTCTGAATCTGTTGAAACATAGGGTAGAATTGCGCCCTGGGCAACAGCAGGGAAGCCCGGCTAtcgatcggaagaggggtgggttggg                                                          >  2:295765/1‑227 (MQ=255)                                                                                                                                                                                                ggAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCGGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTAt                                                                                  <  1:295765/227‑1 (MQ=255)                                                                                                                                                                                                 gAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGGGGCA‑TCTTCCTGCTGAATCTGTTGATGCATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTAt                                                                                  <  1:192657/226‑1 (MQ=255)                                                                                                                                                                                                 gAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTc                                                         >  1:374698/1‑251 (MQ=255)                                                                                                                                                                                                   ggTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGGAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCC‑TCTTCCTGCTGGATCTTTTGATACATAGGATAGAATTGCGCACCGGTCACCAGCCGGCACGCCTTGCTATTTTTCCGGGTCCACCCGCttttccccc                                                       >  2:325778/1‑245 (MQ=255)                                                                                                                                                                                                     tGAAGAGAATTTCTGCCGTATCGTCAGTCGATAGCGGCGGCGCATCGCACAAGGTGGCAGGATGTTGATTTTTCAGGTGAGTAAACGAACTCACGCCATCCTCAGCGGGAAATGGCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGGGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCCGCCGGCACGCCTGGCTATTTTGCAGGAAACACGCGCTTTCCTCCcac                                                     >  1:275996/1‑251 (MQ=255)                                                                                                                                                                                                       aaGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGCGGATAGCACAAGGTGGCAGGTTGTTGATTTTTTAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCCCATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTTCTGCTGAATCTGTTGGTACATAGGATAGAATTGCGCACTGGTCCCCCGCAGGCACGCCCGGCTCTTTTGCCGGAGCCACCCGCTTTCCTccccca                                                    >  2:132429/1‑245 (MQ=255)                                                                                                                                                                                                            aaTTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggc                                              <  1:123503/251‑1 (MQ=255)                                                                                                                                                                                                             aTTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcg                                             <  1:366296/251‑1 (MQ=255)                                                                                                                                                                                                                tCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATATGTCA‑GGCAAATGTGCCGCAATTTAGTGGCA‑TCTTCCTGCTGAATCTGTTGATAAATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcgggg                                          >  1:81938/1‑250 (MQ=255)                                                                                                                                                                                                                cccgccgTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcgggc                                          <  2:374698/250‑1 (MQ=255)                                                                                                                                                                                                                 ccgccgTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcgggcg                                         <  1:400550/251‑1 (MQ=255)                                                                                                                                                                                                                    ccgTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCGCGCACAACAGGCGGGCGTTa                                      <  1:148305/251‑1 (MQ=255)                                                                                                                                                                                                                     cgTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCATGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCCAGTGAGTGTCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAACTGCGCACTGGTCACCAGCAGGCTCGCCTGGCTAGTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGAGCGTTaa                                     <  2:39985/251‑1 (MQ=255)                                                                                                                                                                                                                      gTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAAt                                    <  1:132429/251‑1 (MQ=255)                                                                                                                                                                                                                          cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGc                                <  1:328070/251‑1 (MQ=255)                                                                                                                                                                                                                          cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATGTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGc                                <  2:73550/251‑1 (MQ=255)                                                                                                                                                                                                                                   aTAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATc                       <  2:10968/251‑1 (MQ=255)                                                                                                                                                                                                                                     aGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGGTTAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCGTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATcgc                     <  1:288150/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggcggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGAACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGcc                  <  1:131555/251‑1 (MQ=255)                                                                                                                                                                                                                                            gTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAAtt              <  1:227520/251‑1 (MQ=255)                                                                                                                                                                                                                                                   gCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCa       <  1:178955/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cacaAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGcc        >  1:169351/1‑249 (MQ=255)                                                                                                                                                                                                                                                        aGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCg  >  1:91584/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          ACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCA‑TCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCG  >  NC_000913/37088‑37585

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