Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,844	0	A	C	51.5%	73.5 / 116.2	68	L728V (TTA→GTA)	arcB	aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (28/5);  new base C (0/35);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 2.38e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-10

TTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTC  >  NC_000913/3350696‑3350993                                                                                                                                                      |                                                                                                                                                      ttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAAc                                                                                                                                                      <  1:917086/150‑1 (MQ=255)        cTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTg                                                                                                                                               >  2:1010658/1‑150 (MQ=255)         ttttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAc                                                                                                                                             >  2:298390/1‑151 (MQ=255)          tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGa                                                                                                                                              <  1:285899/149‑1 (MQ=255)            tGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGAcccc                                                                                                                                          >  2:841641/1‑151 (MQ=255)               cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                                       <  1:18201/151‑1 (MQ=255)                        tttCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCttt                                                                                                                               >  1:137466/1‑150 (MQ=255)                              ctcTTCTACGTCGTGCCCCCACTCTTTTTTCATCTCTTCAACCCATCCCCCGCCTTTTTCTCCCCGGGCGGGAAGGTCAGG‑AGACTGAATTTGCTGACCCTGTTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAAtttt                                                                                                                         <  2:1043238/148‑1 (MQ=255)                                 ttCTCCGTCGTGCCCCCCCTCTTCTTTCATCTCTTCAATCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCGGG‑AGACTAATTTTGCTGCCCCATTTTTTGCAGATGGGGAACCCCCACGGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                       <  2:31768/149‑1 (MQ=255)                                 ttCTACGTCGTGACGCCCCTCTTCTTTCATCTCTCCAATCCATCCCCCGACGTTATCTCCCCGGGCCGGAAGGTCAGG‑AGACTAATTTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                       <  1:1024873/149‑1 (MQ=255)                                   cTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcc                                                                                                                   >  1:719728/1‑151 (MQ=255)                                    tACGTCGTGCCCCCACTCTTCTTTCATCTCTTAAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGG‑AGCCTGATTTTGCGGACCAAGTTTTTGCAGTGGGCGTACCCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                         <  2:1012550/144‑1 (MQ=255)                                    acgttcgtacccccACTCTTTTTTAATCTCTTAACCCAATCCCCCAACTTTATCTTCCCCGGCGGGAAGGTCAGG‑AGCCTGAATTTGCTCCCCTGGCTGTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAAtttt                                                                                                                         <  2:1062948/132‑1 (MQ=255)                                      cgtcgtGCCCCCCCTCTTCTTTCATCTCTCAAACCCATCCCCCGACGTTATCTCCCCAGGCCGGAGGGTCAGG‑AGACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCCCTGACCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                 <  1:990965/150‑1 (MQ=255)                                      cgccgtGCCCCCCCTCTTCTTTCATCTCTTAAATCGTTACCCCGCCTTTATCTCCCCAGGCCGGAAGTTAAGG‑ATACTGATTTTGCTGCCCAAGCTGTTGCAGAGGGCGAACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGtc                                                                                                                    <  1:687050/144‑1 (MQ=255)                                         cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                             >  2:928750/1‑151 (MQ=255)                                           tgccccccctctTTTTCCATCTCTCAAACCCATCCCCCGCCTTTATCTCCCCAGGCCGAAAGGTCAGG‑AAACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                              <  1:402062/140‑1 (MQ=255)                                            gccccccttcttttTTCCCCCCTCCAATCCATTCCCCGACGTTATCTCCCCAGGCGGGAAGGTCAGG‑AGACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                 <  1:125333/136‑1 (MQ=255)                                            gACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                          >  2:633692/1‑151 (MQ=255)                                                  aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                    >  2:503674/1‑151 (MQ=255)                                                      ttcttcCCTCTCTCCAACCCATCCCCCGCCTTTATCTCCCCAGGCCGGAGGGTCAGG‑AGCCTGATTGTGCTGCCCCGGTTTTTGCAGATGGCCAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtt                                                                                                  <  2:51930/149‑1 (MQ=255)                                                       tcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                               >  2:749447/1‑151 (MQ=255)                                                           tCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                           >  1:678964/1‑151 (MQ=255)                                                           tCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                           >  2:487657/1‑151 (MQ=255)                                                               ctctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                       >  2:361359/1‑151 (MQ=255)                                                                 ctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcgc                                                                                     >  2:1034684/1‑151 (MQ=255)                                                                     aacccttcccccGACTTTATCTCCCCGGGCGGGAAGGCCAGG‑AGACTAATTTTGCGGCCCCAGTTTTTGCAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGt                                                                                  <  2:313899/140‑1 (MQ=255)                                                                     aacccttcccccGACTTTATCTCCCCAGGCGGGAAGGTGGGG‑AGACTAAATTTGCTCACCCGGTTTTTGCAGTGGGCATACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                 <  2:478053/141‑1 (MQ=255)                                                                     aacccttcccccGACGTTATCTCCCCAGGCCGAAAGGTCAGG‑AGACTAATTTTGCTCCCCCAGCTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                        <  2:378613/134‑1 (MQ=255)                                                                     aacccATCCCCCGACTTTATCTTCCCAGGCGGGAAGGTCAGG‑AGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCATACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                        <  1:298446/141‑1 (MQ=255)                                                                        caatgcACCGACGTTATCGTCCCAGGCGGAGAGGGCAGG‑AGACGGAATTTGCGGGCCCAGGTTTTGCAGATGCCGTACCCCCACGGGCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                                 <  1:533941/127‑1 (MQ=255)                                                                           ttCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                           >  1:882680/1‑151 (MQ=255)                                                                            cccccGTCGTTTTCTTCCAAGCCCGGAAGG‑CAGGTGGACGGACTTTGCTAACCCGCTTGTTGTAGAGGGCGAAACCCCATCGTCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                          <  1:226251/148‑1 (MQ=255)                                                                             ccccGACTTTATCTCCCCAGGCGGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                           <  1:974047/147‑1 (MQ=255)                                                                             ccccGACTTTATCTCCCCAGGCCGGAGGGTCAGG‑AGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCTAACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                           <  2:350039/147‑1 (MQ=255)                                                                             ccccGACGTTATCCCCCCAGGCCGGAAGGCCAGG‑AGACTAAATTTGCTCCCCCGGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                          <  1:239025/148‑1 (MQ=255)                                                                                cGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTc                                                                      >  2:1078174/1‑151 (MQ=255)                                                                                      tATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                >  2:573865/1‑151 (MQ=255)                                                                                      tATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCAc                                                                >  2:234881/1‑151 (MQ=255)                                                                                         cTTCCCAGGCCGGACGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCt                                                             >  2:593561/1‑151 (MQ=255)                                                                                            cccAGGCCGGAAGGTCAGG‑AACCTAAATTTCCTCCCCCAGTTTTTGCAGATGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                          <  2:613700/151‑1 (MQ=255)                                                                                            ccagggCGGGAAGGTCGGG‑AGACTGAATTTGCTGACCCAGTTGTTGAAGAGGGCGAACCCCCACGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGa                                                           <  1:301891/146‑1 (MQ=255)                                                                                             ccAGGCCGAAGGGCAAGG‑AGCTGGATTTTGCTGCCCCATTTTTTGTAGATGGCTTACCCCCATTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACa                                                         <  2:1068470/151‑1 (MQ=255)                                                                                                    ggAAGGTCAGG‑AGACTGAATTTCCTGACCCAGTTTTTGAAGGTGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                   <  1:762737/150‑1 (MQ=255)                                                                                                    gaaaGGTCAGG‑AGACGAATTTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCg                                                  <  2:140061/149‑1 (MQ=255)                                                                                                      agggTCGGG‑AGACTAATTTTGCTCCCCCAGCTGTTGCAGAGGGCTAACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCgg                                                 <  1:18150/148‑1 (MQ=255)                                                                                                      aaGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGc                                                >  1:239378/1‑151 (MQ=255)                                                                                                         ggttagcaggcggattttGCTGACCCATCTGTTCGAGAGGGCGAACCCCCCCTGTCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcat                                              <  1:1057251/136‑1 (MQ=255)                                                                                                            aGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCatcatc                                          >  1:459323/1‑151 (MQ=255)                                                                                                              g‑agaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                        >  1:272758/1‑151 (MQ=255)                                                                                                                 gaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttct                                      >  2:555837/1‑151 (MQ=255)                                                                                                                  aCTGATTTTGCTGCCCCAGTTGTGGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTtctc                                     <  1:593686/151‑1 (MQ=255)                                                                                                                     gAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaa                                   >  1:1046908/1‑150 (MQ=255)                                                                                                                      attttCCTGACCCGGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCa                                    <  2:356069/146‑1 (MQ=255)                                                                                                                         ttGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                              >  2:535626/1‑151 (MQ=255)                                                                                                                           gCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                             <  2:213538/150‑1 (MQ=255)                                                                                                                           gCTGACCCAGCTTTTGCAGAGGGCGAACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                              <  2:418942/149‑1 (MQ=255)                                                                                                                           gCTCACCCAGTTGTTGGAGAGGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                             <  1:642230/150‑1 (MQ=255)                                                                                                                            ctcccccGGCTGTTGCAGATGGCGTACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                            <  2:838347/146‑1 (MQ=255)                                                                                                                             tgccccGGCTTTTGCAGATGGCTTAACCCCTCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                            <  1:488550/146‑1 (MQ=255)                                                                                                                             tcccccAGTTTTTGCAGATGGCTAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                              <  2:1118016/144‑1 (MQ=255)                                                                                                                               aCCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc                        >  1:801375/1‑151 (MQ=255)                                                                                                                                  cAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCg                     >  2:644270/1‑151 (MQ=255)                                                                                                                                   aGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                    <  1:344217/151‑1 (MQ=255)                                                                                                                                     cTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                  >  1:256128/1‑151 (MQ=255)                                                                                                                                          gCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg               >  2:386942/1‑149 (MQ=255)                                                                                                                                            aGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTg               <  1:742548/147‑1 (MQ=255)                                                                                                                                             gATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa           <  1:1010794/150‑1 (MQ=255)                                                                                                                                                     aaCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc  >  1:630439/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGG‑AGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTC  >  NC_000913/3350696‑3350993

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