Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA41.5% 66.0 / 27.7 54Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (27/4);  new base A (0/22);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 3.23e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC  >  NC_000913/436366‑436662
                                                                                                                                                            |                                                                                                                                                    
cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                      >  1:418175/1‑151 (MQ=255)
   tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGAAACCTGGGCGTACCGTTTACCTGTATCGGGCaa                                                                                                                                                   >  2:220474/1‑150 (MQ=255)
      cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGa                                                                                                                                              <  2:28845/151‑1 (MQ=255)
        ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                                            <  2:53046/151‑1 (MQ=255)
            gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgcc                                                                                                                                         <  1:813288/150‑1 (MQ=255)
                   cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCg                                                                                                                                  >  2:740228/1‑150 (MQ=255)
                    ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGa                                                                                                                                 >  1:654476/1‑150 (MQ=255)
                    ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGGACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGa                                                                                                                                 >  2:692103/1‑150 (MQ=255)
                          ggggcggcggggggCGCTCTCGGGGGGAGAAGATTGCGGTTTGTTTTTCCCTGTCCCGAAACTGACCCGGGGCGCGTTGGGGGTGCCTCCCGCCCCCTGGGGGGTCCTTTTTACCTTTTGGGGAAAAATGACCGCCGATATCGAAGGGCtt                                                                                                                          <  2:1180320/140‑1 (MQ=255)
                           gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCttt                                                                                                                         >  1:402729/1‑151 (MQ=255)
                                 cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                    >  1:979128/1‑150 (MQ=255)
                                 cgcggGGCGCCCTCGGGGGGTGAAGATTCGGGTTTGTTTTTCCCTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCCCGGACCCTGGGGCGTCCTTTTTTCTTGTACCGGAAAATTGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                    <  2:1106924/150‑1 (MQ=255)
                                  ggtgGGCCTTCTCGGGCGGTGAAGTTTCCGATTTGTGTTTCCCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCCCTCGGCCACCTGGGGGTCCCTTTTCCTTTTATGGGGAAAAAGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                  <  1:817194/149‑1 (MQ=255)
                                     gggcctctgcggggggTGAAGATTAGGGTTTTTGTTTCCCTGGGCCGGAACTGAACGGTGGCGCGTTGGTTGTGCCCCGCGCACCCGGGGGGTTCCTTTTTCCCTTTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                 <  2:1078737/135‑1 (MQ=255)
                                     gggCGCTCTGTGGCGGGAAAGTTTACGTGTTTTTTTTCTGTGTGCGGGAACTGACCGGGGGCGCGCTGGGTGTGGCTCTCGGACCCTGGGGCGTCCTGTTTCCCTGTGGCGGGAAAATGACCGCCGATATCGAAGGGCttt                                                                                                                         <  2:254210/141‑1 (MQ=255)
                                        cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                             >  1:1168517/1‑150 (MQ=255)
                                              ctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                        >  1:804978/1‑149 (MQ=255)
                                                ggggttGAAGATTCGGATTTGTTTTGTATGGTCCCGGAACTGCTCCGGGGCGCGCGGGTTGTCTCCTTCGGCGACGGGGCTTTCCGTTTTCCCTTTGTCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                     <  2:784792/145‑1 (MQ=255)
                                                ggCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGACCCGGGGCGCGTTGGTTGTGGCTCTCGGACACTGGGGCGTACCGTTTCCTTTTATGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                        <  1:1032022/147‑1 (MQ=255)
                                                ggCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                     >  2:203604/1‑150 (MQ=255)
                                                   ggtgaagtttcgggtttGTGTTTCACTGTGCGGGAACTGACCGGGGGCGCGTGGGATGGGTCTCTCGGACCCCGGGGCGTCCTTTTTCCTTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                                 <  2:815475/136‑1 (MQ=255)
                                                   ggTAAAAATTTCGATTTGTTTTTCACTGTGCGAGAACTGACCTGCGGCCCGTGGGGTGGGGCTCGCGGCCACGTGGGCGTACTTTTTACCTGTACCGGGAAAGTCACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                  <  2:349134/150‑1 (MQ=255)
                                                        aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                             >  2:791576/1‑150 (MQ=255)
                                                            taggAGTTTTTTTTCCCTGGGCCGGAACTGACCGGGGGCGCGTTGGTTGTGGCTCTGGGACCCCGGGGCGTCCCTTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                         <  1:1050362/147‑1 (MQ=255)
                                                            tacgttttGTGTTTCCCTGTGCCGGAACGGACCGGGGGCGCGTTGGTGGTGCCTCCGGGACACCTGGGCTGCCCTTTTCCCTTTGCGGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                         <  2:435075/144‑1 (MQ=255)
                                                            tACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGCCACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                          >  2:2904/1‑149 (MQ=255)
                                                              cGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                                      >  2:140274/1‑151 (MQ=255)
                                                                ggttgttTTTCCTTGTGCCGAAACGGACCCGTGGCGCGGTGGGTGGGGCTCCGGGAACCTGGGGCGTCCCTTTTTCCTGTACCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                        <  1:634064/146‑1 (MQ=255)
                                                                   tgtgTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                  >  1:415923/1‑150 (MQ=255)
                                                                    tttttttcccgggCCGGAAACTGACCCGGGGCGCGTTGGTTGTGCCTCTCGGCCCCCGGGGCGTCCCTTTTCCTTGTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                  <  1:1086262/138‑1 (MQ=255)
                                                                          ctcTGTGCGGGAACGAACCGGGGGGCGGTGGGGTGTGCCTCCGGGCCCCTGGGGGTTCCCTTTTACCTGTATGGGAAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                           <  2:189651/148‑1 (MQ=255)
                                                                          cACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                          >  2:788386/1‑151 (MQ=255)
                                                                               tgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                      >  1:380946/1‑150 (MQ=255)
                                                                                    gaaaCTGACCCTTGGCGGGTTGGTTGTGTCTCCGGAACCCCTGGGCGTCCCTTTTCCTTTTCCGGGGAAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                   <  1:396090/147‑1 (MQ=255)
                                                                                      accggcccggTGGCGCGTTGGATGTGGCTCCCGGCCACCTGGGCGTACCTTTTCCCTGTATGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                 <  2:262218/140‑1 (MQ=255)
                                                                                      aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                               >  1:468079/1‑151 (MQ=255)
                                                                                         tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                            >  2:200081/1‑151 (MQ=255)
                                                                                          gAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                           >  2:261005/1‑151 (MQ=255)
                                                                                              cgtgcggcgcTGGAGGTGCTCCTCGTACCCTGGGGCGTACTTTTTCCTTTTGTCGGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                           <  2:1222531/141‑1 (MQ=255)
                                                                                                 ggCGGGTTGGGTGTGCCTCTCGTCCACCTGGGTTTCCCGTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                            <  1:703850/143‑1 (MQ=255)
                                                                                                     cgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGTCGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                >  2:322051/1‑151 (MQ=255)
                                                                                                      ggtGGATGTGTCTCTCGGCCCCTGGGGCGTCCTTTTTCCTTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                <  2:243903/148‑1 (MQ=255)
                                                                                                      gcTGGATGTGGCTATCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                                >  2:651959/1‑150 (MQ=255)
                                                                                                        tGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                             >  1:1188573/1‑151 (MQ=255)
                                                                                                           gttggtgcctccGGACCCCTGGGGCGTCCCTTTTACCTTTACGGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                              <  2:95360/138‑1 (MQ=255)
                                                                                                               tgGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaa                                        >  2:937685/1‑150 (MQ=255)
                                                                                                                ggCTCTGGGCCCCCGGGGCGTCCCTTTTCCTTTTAGGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaaa                                       <  1:284147/150‑1 (MQ=255)
                                                                                                                    ctcGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                  >  1:851242/1‑151 (MQ=255)
                                                                                                                             cgggggCGTCCCTTTTCCTTGTAGGGGAAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa                           <  1:1132430/147‑1 (MQ=255)
                                                                                                                                             tttCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg                <  1:920961/148‑1 (MQ=255)
                                                                                                                                              ttACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc            >  1:433743/1‑151 (MQ=255)
                                                                                                                                                  cTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt        >  1:1273077/1‑151 (MQ=255)
                                                                                                                                                     tATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt     >  1:637848/1‑151 (MQ=255)
                                                                                                                                                        cGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGc  >  2:1178118/1‑151 (MQ=255)
                                                                                                                                                         gggAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt    <  1:140276/148‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                                    
CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGC  >  NC_000913/436366‑436662

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