Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,429	0	T	G	25.0%	115.6 / 14.4	60	intergenic (‑128/+95)	thiM/rcnR	hydroxyethylthiazole kinase/DNA‑binding transcriptional repressor RcnR
Reads supporting (aligned to +/- strand):  ref base T (11/34);  new base G (15/0);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.45e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCAC  >  NC_000913/2185281‑2185576                                                                                                                                                     |                                                                                                                                                     aCCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGttt                                                                                                                                                     <  2:2018898/150‑1 (MQ=255)      gcagGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTg                                                                                                                                                <  1:617344/150‑1 (MQ=255)         ggTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCga                                                                                                                                             <  2:1962916/150‑1 (MQ=255)          gTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTCgat                                                                                                                                            >  1:1768931/1‑150 (MQ=255)           tCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgatg                                                                                                                                           <  2:736901/150‑1 (MQ=255)              aCTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGATa                                                                                                                                        >  1:1846277/1‑150 (MQ=255)               cTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTCGATGATaa                                                                                                                                       >  1:334163/1‑150 (MQ=255)                 tGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAAt                                                                                                                                     <  2:2098544/150‑1 (MQ=255)                 tGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAAt                                                                                                                                     <  2:2097795/150‑1 (MQ=255)                  gCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGGTTCGATTATAAAta                                                                                                                                    >  1:1184877/1‑149 (MQ=255)                    aTAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATaaaaagc                                                                                                                                  >  1:1531568/1‑146 (MQ=255)                      aGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCag                                                                                                                                <  1:1543363/150‑1 (MQ=255)                       gTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCaga                                                                                                                               <  1:368706/150‑1 (MQ=255)                       gTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGATAAAACGCaga                                                                                                                               >  1:509610/1‑150 (MQ=255)                             tCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGGGTGGTTGCGAGGATAAATAGCAGAggagga                                                                                                                         >  2:151942/1‑150 (MQ=255)                             tCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTACCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGTTAAATCGCAGAggcgga                                                                                                                         >  2:1113087/1‑150 (MQ=255)                                tGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATgg                                                                                                                      <  2:1357493/150‑1 (MQ=255)                                    aTAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AAt                                                                                                                 <  1:956070/150‑1 (MQ=255)                                       aCGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTc                                                                                                              <  1:215007/150‑1 (MQ=255)                                                   aaTGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATAGCAGAGGATGATGGT‑AATGTCCAGCGCACGCGt                                                                                                  >  2:262104/1‑150 (MQ=255)                                                   aaTGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGt                                                                                                  >  2:262297/1‑150 (MQ=255)                                                        ccTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGCTGTTAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTaa                                                                                             >  2:1511575/1‑150 (MQ=255)                                                        ccTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATTATAAAAAGCAGAGGAGGATGGT‑AATGTACAGCGCCCGCGTTTTaa                                                                                             >  1:1520964/1‑150 (MQ=255)                                                            gtggtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGa                                                                                         <  1:1184724/150‑1 (MQ=255)                                                            gtggtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGa                                                                                         <  1:872209/150‑1 (MQ=255)                                                            gtggtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTCGATTATAAATAGCAGAGGAGGATGGTAAATGACC‑GCGCACCCGGTGTAAACGa                                                                                         >  2:977233/1‑150 (MQ=255)                                                               gtCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGa                                                                                      <  2:791045/150‑1 (MQ=255)                                                                  cGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAGCGAGGaaaa                                                                                   <  2:1442653/150‑1 (MQ=255)                                                                  cGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGTTAAATAGCAGAGGGGGATGGT‑AATTTGCAGCGCACCCGGTGTTAACGAGGaaaa                                                                                   >  1:1076957/1‑150 (MQ=255)                                                                        ttCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGTTAAAACGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAAAGGCa                                                                             >  2:412566/1‑150 (MQ=255)                                                                          cccTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGGGATGATAAAAAGCAGAGGAGGAGGGT‑GATGTGCCGCGCACGCGTTGTAAAAGAGGAAAAACGGCAGt                                                                           >  1:697295/1‑150 (MQ=255)                                                                           ccTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGtt                                                                          <  2:205938/150‑1 (MQ=255)                                                                              aCGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTa                                                                       <  2:509610/150‑1 (MQ=255)                                                                                 cTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACaa                                                                    >  2:1993480/1‑150 (MQ=255)                                                                                  tGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAAt                                                                   <  2:948700/150‑1 (MQ=255)                                                                                  tGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTTCGATGATAAAACGCAGAGGAGGTTGGT‑AATTTCCAGCGCACGCGTTGTTAACCAGGAAAAAACGCAGTTTTAAAAt                                                                   >  2:675071/1‑150 (MQ=255)                                                                                  tGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATGATAAATAGCCGAGGAGGATGGT‑AATGTGCAGCGCCCGCGGTTTAAAAGAGGAAAAAAGGCAGTTTTaaaaa                                                                   >  1:475510/1‑149 (MQ=255)                                                                                   ggCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGATTATAAATAGCAGAGGATGATGGT‑AATGTGCAGCGCACGCGGTGTAAAAGAGGTAAAAAGGCAGTTTTaaaaaa                                                                  >  1:41470/1‑148 (MQ=255)                                                                                             aGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTtaataa                                                        <  2:1550266/150‑1 (MQ=255)                                                                                             aGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTtaataa                                                        <  2:463266/150‑1 (MQ=255)                                                                                              gATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATaaa                                                       <  1:1988107/150‑1 (MQ=255)                                                                                                 cAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAatta                                                    <  2:1749012/150‑1 (MQ=255)                                                                                                   ggTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAattatt                                                  <  2:1893529/150‑1 (MQ=255)                                                                                                   ggTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTAAAATCGCGTTAATAAattatt                                                  >  2:1883272/1‑150 (MQ=255)                                                                                                    gTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTAttt                                                 <  1:124039/150‑1 (MQ=255)                                                                                                       aTACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGTTTGCGATTATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGat                                              <  1:665967/150‑1 (MQ=255)                                                                                                       aTACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGTTGCGCTGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAAAGAGGAAAAACGGCAGTTTTACAAACGCGGTAAAAAATTATTTGat                                              >  1:1167596/1‑150 (MQ=255)                                                                                                           gggTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGatatat                                          <  1:1118565/150‑1 (MQ=255)                                                                                                                tttCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATc                                     <  2:1983329/150‑1 (MQ=255)                                                                                                                       gCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGTTAAATAGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATAGCGTTAATAAATTATTTTATATATGAATACAGCAcc                              >  2:1042236/1‑150 (MQ=255)                                                                                                                         cTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCtt                            <  2:547320/150‑1 (MQ=255)                                                                                                                            cACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAg                         <  1:251399/150‑1 (MQ=255)                                                                                                                                     gggCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATc                <  2:1908274/150‑1 (MQ=255)                                                                                                                                     gggCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATc                <  2:1907366/150‑1 (MQ=255)                                                                                                                                        cACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAg             <  2:1184877/150‑1 (MQ=255)                                                                                                                                         aCCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGa            <  2:475510/150‑1 (MQ=255)                                                                                                                                            ccGAGTCGTTTGGTTGCGATGATAAATCGCCGAGGAGGATGGT‑AATGTCCACCGCACGCGTTGTAAAAGAGGAAAAAAGGCAGTTTTAAAAACGCGGTAATAAAATATTTTATATATTAAACCAGCACCATAAGAAAGACAACCAGAtat         >  2:953631/1‑148 (MQ=255)                                                                                                                                              gAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAACACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGAtcttc       <  1:1648276/150‑1 (MQ=255)                                                                                                                                                 tCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCttcttc    <  2:1871333/150‑1 (MQ=255)                                                                                                                                                   gTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCAc  >  1:1508465/1‑150 (MQ=255)                                                                                                                                                   gTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCAc  >  1:467088/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCGCAGAGGAGGATGGT‑AATGTCCAGCGCACGCGTTGTAAACGAGGAAAAACGGCAGTTTTACAATCGCGTTAATAAATTATTTGATATATGAATCCAGCACCTTCAGAACGACATCCAGATCTTCTTCAC  >  NC_000913/2185281‑2185576

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