Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,3680AT42.9% 80.8 / 27.6 56R69R (CGT→CGAyjiRputative DNA‑binding transcriptional regulator/putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (9/23);  new base T (24/0);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 9.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTG  >  NC_000913/4571219‑4571511
                                                                                                                                                     |                                                                                                                                                
ggCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTAc                                                                                                                                                 <  2:913108/151‑1 (MQ=255)
   gcTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGTTACGtt                                                                                                                                              >  1:1033980/1‑151 (MQ=255)
     ttCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCCGGGGTTTTCGtt                                                                                                                                              >  1:819785/1‑149 (MQ=255)
          cATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGAc                                                                                                                                       <  1:209690/151‑1 (MQ=255)
          cATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGAACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGAc                                                                                                                                       <  1:472317/151‑1 (MQ=255)
           aTCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCCTTTCGGTGTGCggc                                                                                                                                       >  1:361411/1‑148 (MQ=255)
            tCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTAGTTTTGCGGCaa                                                                                                                                     >  1:501242/1‑150 (MQ=255)
              cGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACaaaa                                                                                                                                   <  2:819744/151‑1 (MQ=255)
                 tgtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCATCCAGCATATCCAGCATCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCCCCGGACGCGTATTCGGCGGTACTGGTGGAGGGGCTTTACGTTGTGCGTCaaa                                                                                                                                    <  2:448332/147‑1 (MQ=255)
                  gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAAc                                                                                                                              <  2:794649/151‑1 (MQ=255)
                  gtgCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTTGTTGTGCGACAaaaaaa                                                                                                                               >  1:1019423/1‑150 (MQ=255)
                     cTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTTGCTGTGCGGCAAAAAAAacc                                                                                                                            >  2:533458/1‑150 (MQ=255)
                     cTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGTTTCGTTTTGCGACAAAAAAAacc                                                                                                                            >  2:340833/1‑150 (MQ=255)
                       tttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCg                                                                                                                           <  1:89585/149‑1 (MQ=255)
                       tttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGa                                                                                                                          <  1:473448/150‑1 (MQ=255)
                        ttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGGGGGGTTTCGGTGTGCGGCAAAAAAAacc                                                                                                                            >  2:385052/1‑147 (MQ=255)
                          gCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAcgc                                                                                                                      <  1:113665/151‑1 (MQ=255)
                          gCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAcg                                                                                                                       <  1:999600/150‑1 (MQ=255)
                             tAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGTTTTGCGGCAAAAAAAACCCaaaccgg                                                                                                                    >  2:908138/1‑145 (MQ=255)
                                cGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGCTTACGTTGTGGGGCAAAAAAAACCGAACGCgg                                                                                                                    >  1:354504/1‑147 (MQ=255)
                                      aaTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTTCGTTGTGGGGCAAAAAAAACCGAACGCGGCCGCgggg                                                                                                            >  1:583444/1‑147 (MQ=255)
                                       aTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTGCGGCTGGGCGTTTCGTTTTGCGGCAAAAAAACCCCCACGCCGCCGCgggggg                                                                                                          >  2:1026450/1‑146 (MQ=255)
                                       aTGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGCTTTTCGTTTTGCGGCAAAAAAAACCCAACGCGGCTGCGGGGtg                                                                                                          >  1:295234/1‑150 (MQ=255)
                                          gAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGTTTTCGGTGTGCGGCAAAAAAAACCCAACCCGGCCGCgggggg                                                                                                          >  1:304691/1‑143 (MQ=255)
                                           aaCTGTCGCTATGCGCCTCCAGCATATCCAGCGCCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAg                                                                                                     <  2:196866/151‑1 (MQ=255)
                                                             ccAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTTCGTTTTGCGGCCAAAAAAACCCAACGCGGCCGCGGGGGGGTGAGCTTTATCGTCTcc                                                                                        >  2:381659/1‑146 (MQ=255)
                                                                  atatCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGAt                                                                              <  1:434709/151‑1 (MQ=255)
                                                                   tatCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCCGGGGGGTTCGGTGTGGGACAAAAAAAACCCAAAGCGGCTGCGGGGGGATGGGCTTCCACCTCTcccaccgctg                                                                                >  2:61333/1‑140 (MQ=255)
                                                                    atCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGGTGTGGGGCAAAAAAAACCCCAAGCGGCCGCGGGGGGGTGGGCTTCATCGTCTCCCGCGTCCGAt                                                                              >  2:1036028/1‑149 (MQ=255)
                                                                      ccAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGGGGCCTGGCGTTTCGCTTTGGGGGCAAAAAAACCGAAAGCGGCGGCGGGGGGGTGAGCTTCATCGTCCCCCACGTCCGGTAct                                                                           >  2:480910/1‑149 (MQ=255)
                                                                        aGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTGTGGGACAAAAAAAACCCAAAGCGGGTGGGGGGGGGTGAGGTTTATTGTTTCCAGCGGCTGATaaggc                                                                         >  2:68609/1‑146 (MQ=255)
                                                                            ccTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctgct                                                                     <  1:263411/150‑1 (MQ=255)
                                                                                  cccACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGt                                                              <  2:473288/151‑1 (MQ=255)
                                                                                           tAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCCTTGTGCGGCAAAAAAACCCCAAACCGGCCGCGGGGGGATGAGGCTCATCGTCTCCCGCGGCTGATACGCCCGCTTCCCGGGGGTGAt                                                        >  1:983740/1‑148 (MQ=255)
                                                                                            aaTTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACt                                                     <  1:487555/150‑1 (MQ=255)
                                                                                             aTTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGa                                                   <  1:542770/151‑1 (MQ=255)
                                                                                                  cACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTGTGCGGCAAAAAAAACCCAACGCGGCGGGGGAGGGATGGGGTTCATCGTCTCCCGCGTCTGATACGCCTGCCGCCCGGGGGTGGTAATGACGcc                                               >  1:901000/1‑150 (MQ=255)
                                                                                                   aCCGGGCGCTGCACCGGACGCGTCATCGGCGGTTCTGGCGGGCGGGCTTTTCGGTGTGCGGCAAAAAAAACCCAAACCGGCCGCGGGGGGGTGAGCTTTCTCGTCTCCCGCGGCTGATACGCCTGCTGCAAGGGGCTGATACTGACGcc                                               >  1:189359/1‑149 (MQ=255)
                                                                                                      gggCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCt                                          <  2:442299/151‑1 (MQ=255)
                                                                                                       ggCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTTTGCGGCAAAAAAAACCCAAACCGGCCGCGGGGGGGTGAGGTTCATCCTCTCCCGCGGCTGATACGCCCCCTGCCCGGGGGCGGTACTGAAGcc                                               >  2:555235/1‑145 (MQ=255)
                                                                                                         cgcTGCACCGGACGCGTCATCGGCGGTACTGGCGGGTGGGGTTTTCGGTGTGGGACAAAAAAACCCCGAAGCGGGTGCGGGGGGGTGAGGCTCCTCGTCTCCCCCGTCTGATACCCCCGCTGCACGGGGCTGGTACTGACGCCCCGCTcc                                        >  1:887285/1‑150 (MQ=255)
                                                                                                             gCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACt                                   <  2:899613/151‑1 (MQ=255)
                                                                                                              cACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGTTTTCGTTGTGCGACAAAAAAAACCCAACGCGGCTGCGGGGGGGTTAGCTTCATCGTCTCCAGCGTCTGATACGCCCGCTGCACGGGGCTGATACTGACGCCCTGCTCCTGACtt                                  >  1:475291/1‑151 (MQ=255)
                                                                                                               aCCGGACGCGTCATCGGCGGTACTTGCGGCTGGGGTTTTCGGTGTGGGGCCAAATAAACCCGAAGCGGGTTCGGGGGGATGAGGTTTCTTGTCTCCCGCCTCTTATAACCCCGCTGCCCGGGGGGGGTACTGACGCCCTGCTCCTTACtt                                  >  1:105785/1‑150 (MQ=255)
                                                                                                                    aCGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTg                            <  1:183867/151‑1 (MQ=255)
                                                                                                                     cgcgTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTgc                           <  2:492182/151‑1 (MQ=255)
                                                                                                                     cgcgTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTgc                           <  1:266455/151‑1 (MQ=255)
                                                                                                                              ggcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACg                  >  1:914363/1‑151 (MQ=255)
                                                                                                                               gcggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg                 <  2:304675/151‑1 (MQ=255)
                                                                                                                                cggTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACgg                 >  1:957843/1‑150 (MQ=255)
                                                                                                                                   tACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAt             >  1:562834/1‑151 (MQ=255)
                                                                                                                                     cTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAttt           >  1:254272/1‑151 (MQ=255)
                                                                                                                                         cggcTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTcc       <  1:1022185/151‑1 (MQ=255)
                                                                                                                                          ggcTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTccc      >  2:954998/1‑151 (MQ=255)
                                                                                                                                            cTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCg    >  1:1011432/1‑151 (MQ=255)
                                                                                                                                              gggCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTg  >  1:136487/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTG  >  NC_000913/4571219‑4571511

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