Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,038,477	0	G	A	78.6%	90.5 / 19.1	43	A246A (GCG→GCA)	appC	cytochrome bd‑II subunit 1
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (14/19);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCTTGCAACGCGGACGTATGGCCTGGC  >  NC_000913/1038250‑1038715                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                ttCAGCGCGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGGGACCGGGGCCATGGTTATTATGGCGATCAGCGCCTGGGATTTATTGCGCGGACGGGAGGGAAATGTCGCGTTACCCTAGTTTGCAATCGGTTTCGTCTTCGGGATTCTGGCGGTTATCGGTATCCTGCCACTCGGAGACAGGTCTGCGTATGGAGT‑GCGCAAAAGACAACCGGTAAAAcaggg                                                                                                                                                                                                                         >  2:57402/1‑247 (MQ=255)   cAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTggcg                                                                                                                                                                                                                        <  1:11312/250‑1 (MQ=255)           tGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAAc                                                                                                                                                                                                                                      <  1:229879/228‑1 (MQ=255)           tGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTAAGTGACCGGGGCCATGTTTTTTATGGCGATCCGCGCCTGGGATTTACTGCGCGGGCGGGGGCGCAATGTCGCATTACGCTCGTTTTCCCTCGGGTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGGGACAGTTCTTCGTTTGACGTCGCCAAAG‑TACAAc                                                                                                                                                                                                                                      >  2:229879/1‑228 (MQ=255)                            aGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGGATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTAATCTGGCGATTATCGGGACCCTGCAACTCGGAGACAGTTCTGCGGATGAAGTCGCAAAAG‑GACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGGGGCaaaagg                                                                                                                                                                                               >  1:57711/1‑247 (MQ=255)                                           gtgCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTGCTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGGCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACctgcacctgcac                                                                                                                                                                                <  2:50284/250‑1 (MQ=255)                                           gggNACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTGCTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCtgcacctgcacc                                                                                                                                                                               <  2:444761/249‑1 (MQ=255)                                                                tACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGcc                                                                                                                                                          >  2:58651/1‑251 (MQ=255)                                                                tACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGcc                                                                                                                                                          >  1:359718/1‑251 (MQ=255)                                                                 aCGTGACCGTGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGcc                                                                                                                                                          <  1:362225/250‑1 (MQ=255)                                                                 aCGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCg                                                                                                                                                         <  2:150059/251‑1 (MQ=255)                                                                   gTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCgg                                                                                                                                                        >  2:423509/1‑250 (MQ=255)                                                                        cGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTggcg                                                                                                                                                                                                                        <  2:163273/181‑1 (MQ=255)                                                                        cGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTggcg                                                                                                                                                                                                                        >  1:163273/1‑181 (MQ=255)                                                                                       tatGGCGATTAGAGCTGGGTACTTACTGCGCGGCCGGGAGCGCAATGTCGCATGCCGCTCGTTTGCCATCGGTACCGTCTTCGGTAATCCGGCGATTATCGGTATCCTGCAACTCGGAGACAGTTATGCGTATTACGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGTGATGGAAGGGGAGTGGCAAACGGACCCTGCACCTGCTCCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAAc                                                                                                                                    <  1:157868/250‑1 (MQ=255)                                                                                        attgCGATAAGCGCCTGGTATTTACTGCGCGGACGGGCGCGCAATGTCGCATTACGCTCGTTTGCCATCTGTTCCTTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCTCCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACg                                                                                                                                   <  2:179378/247‑1 (MQ=255)                                                                                        atGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGc                                                                                                                                  <  1:343027/251‑1 (MQ=255)                                                                                          ggCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACNCCt                                                                                                                                >  2:437238/1‑251 (MQ=255)                                                                                          ggCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCt                                                                                                                                <  1:235143/251‑1 (MQ=255)                                                                                           gCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGTTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGTATCAAGAGCGTAACGCCtt                                                                                                                               <  1:230930/251‑1 (MQ=255)                                                                                            cGATCAGCGCCTGGTATTTACTGCGCTGACGGGAGCCCAATGGCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTTGCGATTATCGGTACCCTTCAACTCGGAGACATTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGTCCGGAACAGGATCAAGAGCGTAACGCCttt                                                                                                                              <  1:176309/251‑1 (MQ=255)                                                                                                     ccTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCaa                                                                                                                      >  1:132379/1‑250 (MQ=255)                                                                                                       tGGTATTTACTGCGCGGGCGGGAGCTCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGATCCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCaaa                                                                                                                     <  2:347770/249‑1 (MQ=255)                                                                                                        ggTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAAt                                                                                                                   <  2:84178/250‑1 (MQ=255)                                                                                                                    gcgGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGgcggcggttga                                                                                                                                                                                                                >  2:86539/1‑140 (MQ=255)                                                                                                                    gcgGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGa                                                                                                                                                                                                                <  1:86539/145‑1 (MQ=255)                                                                                                                     cgGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCa                                                                                                                       >  1:455627/1‑233 (MQ=255)                                                                                                                     cgGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGGAACGCCTTTGCCCTCa                                                                                                                       <  2:455627/233‑1 (MQ=255)                                                                                                                      ggACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAgg                                                                                                    >  1:94967/1‑251 (MQ=255)                                                                                                                                        gCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGACCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCactcactc                                                                                  <  2:105677/251‑1 (MQ=255)                                                                                                                                            tACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCAtt                                                                               >  2:197414/1‑250 (MQ=255)                                                                                                                                                    tttGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGGAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGCGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAAcc                                                                      >  1:19006/1‑251 (MQ=255)                                                                                                                                                        ccATCGGTTACGTCTTCGGTACTCTAGCGATTATCGGTAACCTGCAACTCGGACACAGTTATGCGTATGAAGTCGCACAAA‑TACAACCGGTAAAACTGGCGGCGATGGACGGGGAGTGGAAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGACTTTGCCCTCAATATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCACTCGATAAAcccatg                                                                  >  2:281664/1‑248 (MQ=255)                                                                                                                                                            cGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCgg                                                              <  1:346465/251‑1 (MQ=255)                                                                                                                                                                                   cGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGTCTGaa                                       <  2:57711/251‑1 (MQ=255)                                                                                                                                                                                            gTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCa                              <  2:308108/251‑1 (MQ=255)                                                                                                                                                                                              aCCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACg                            <  1:354366/251‑1 (MQ=255)                                                                                                                                                                                                 cTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCtt                         >  1:369518/1‑251 (MQ=255)                                                                                                                                                                                                   gCAACTCGGAGACAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTTCCGGGTCTTAATAATTTGATGGCTGAAGCCTACCCACGCtt                         <  1:402378/249‑1 (MQ=255)                                                                                                                                                                                                              aCAGTTCTGCGTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCTTGCAACGCGGACGt            >  2:337384/1‑251 (MQ=255)                                                                                                                                                                                                                        gTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCTTGCAACGCGGACGTATGGCCTGGc  >  2:72720/1‑251 (MQ=255)                                                                                                                                                                                                                        gTATGAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGAATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCTTGCAACGCGGACGTATGGCCTGGc  >  2:328530/1‑251 (MQ=255)                                                                                                                                                                                                                            gAAGTCGCAAAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGTAGTGGGAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGAt                                                                                                                                                <  1:73518/105‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                TTCAGCGAGCTGGTCTTTAATCCGGTCAGCCAGGTGAAATTTGTGCACACCGTAATGGCGGGCTACGTGACCGGGGCCATGTTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCGGAGACAGTTCTGCGTATGAAGTCGCGCAAG‑TACAACCGGTAAAACTGGCGGCGATGGAAGGGGAGTGGCAAACGGAACCTGCACCTGCACCGTTCCATGTGGTTGCCTGGCCGGAACAGGATCAAGAGCGTAACGCCTTTGCCCTCAAAATTCCCGCGCTGCTAGGGATCCTCGCCACTCACTCATTAGATAAACCCGTGCCGGGTCTGAAGAATTTGATGGCTGAAACCTACCCACGCTTGCAACGCGGACGTATGGCCTGGC  >  NC_000913/1038250‑1038715

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: ‑