Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5040AC49.1% 58.7 / 19.9 53T291P (ACC→CCC) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base C (0/26);  total (16/37)
Fisher's exact test for biased strand distribution p-value = 1.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCG  >  NC_000913/436355‑436648
                                                                                                                                                      |                                                                                                                                                
tttCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTa                                                                                                                                                  <  2:641735/150‑1 (MQ=255)
     gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGt                                                                                                                                             <  1:656039/150‑1 (MQ=255)
        gcgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCg                                                                                                                                         <  1:445046/151‑1 (MQ=255)
         cgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCgg                                                                                                                                        <  1:332759/151‑1 (MQ=255)
                 cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                  <  1:928075/149‑1 (MQ=255)
                   ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                              >  1:715286/1‑151 (MQ=255)
                   ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                              <  2:494037/151‑1 (MQ=255)
                    cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                  <  2:596027/146‑1 (MQ=255)
                         tGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGat                                                                                                                        >  2:665964/1‑151 (MQ=255)
                            aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATc                                                                                                                     >  2:426239/1‑151 (MQ=255)
                               ggAACAGGCGCGGCGCTGGGCCCTCTCGGGGGGTGAAGATTACGGGTTGTTTTTCCCTGTGCCGGAACTGAACCGTGGCGCGCTGGTTGTGGCTCCCGCACACGTGGGCGTACCTTTTCCTTGTATCGGGAAAATGACCGCCGATATCGaa                                                                                                                  <  1:542973/151‑1 (MQ=255)
                                    aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                             >  1:262969/1‑151 (MQ=255)
                                    aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                             >  2:899611/1‑151 (MQ=255)
                                                gggCGTTCTCGGGGGGAGAAGTTTACGGGTTGTTTTTCCCTGTGCGGAAACTGACCGGTGGCGCGTTGGATTTGGCTCTCGGACACGTGGGCTTCCCTTTTCCCTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTg                                                                                                          <  2:432959/142‑1 (MQ=255)
                                                gggCGCTTTCGGGGGGGGAAGATTACGATTTGTTTTTCACTGTGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCGTTTCCCTGTATCGGGAAAATAACCGCCGATATCGAAGGGCTTTGtttt                                                                                                      <  2:555336/146‑1 (MQ=255)
                                                gggCGCTCTCGGGGGGTGAAGATTACGGGTTGTGTTTCACTGTGCCGAAACTGACCCGTGGCGGGCGGGATGTGGCTCCCGGCCACCGGGGCGTCCTTTTTACCTGTATCGGGCAAATGCCCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                   <  1:262454/149‑1 (MQ=255)
                                                gggCGCTCTCGGGCGGTGAAGATTACGGGTTGTGTTTCACGGTGCCGAAACTGAACCGTGGCGCGCGGGATGTGGCTCCGGGCCCCCTGGGGGTACCTTTTCCCTGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                    <  2:154651/148‑1 (MQ=255)
                                                           gggggTGAAGTTTACGTGTTGTTTTTCCCTGGGCGGGAACTGACCCGTGGCGCGCGGGATGTGGCTCTCGGCCACCGGGGCGTCCTTTTTACTTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                        <  1:427419/146‑1 (MQ=255)
                                                           gggggTGAAGATTACGAGTTGTGTTTCATTGTCCCGGAACGAACCGGGGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCTTTTTACCTTTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                       <  2:745410/147‑1 (MQ=255)
                                                           gggggTGAAAATTAGGGGTTGTTTTTCCTTGTGCCGGAACTGACCCGGGGCGCGGTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTTTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                          <  2:721620/144‑1 (MQ=255)
                                                           ggCGGTGAAGTTTACGAGTTGTTTTTCACTGTGCCGGAACTGACCGGTGGCGCGTGGGATGTGGCTCTCGGACCCGGGGGCGTCCCGTTTCCTTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                        <  2:1011397/149‑1 (MQ=255)
                                                           ggCGGTAAAAATTACGGGTTGTGTTTCACTGTGCCGAAACTGACCGGTGGCGCGCTGGATGGGCCTCTCGGACCCGGGGGCGTCCCGTTTCCTTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                            <  2:21655/145‑1 (MQ=255)
                                                             cGGTGAAGATTACGTGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGTCGTACCGTTTACCTGGATAGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaagc                                                                                     >  1:728637/1‑148 (MQ=255)
                                                             cGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                    >  2:384605/1‑151 (MQ=255)
                                                              ggggAAGATTACGAGTTGTGTTCCGCGGTGCCGGACCTGACCCGTGGCGCGCTGGAGGGGGCTCTCGGCCCCCGGGGCGCACCTTTTCCCTTTAGCGGACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                   <  1:412604/148‑1 (MQ=255)
                                                              ggTGAAGTTTACGGTTTGTTTTTCACTGTGCCGGAACTGAACGGTGGCGCGTGGGATGTGGCTCGCGGACACCGGGGCGTCCCTTTTCCCTGTACCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                    <  2:550911/150‑1 (MQ=255)
                                                                tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                 >  1:538226/1‑151 (MQ=255)
                                                                  aaaattatgttttttgtTTCACTTTCCCGGATCTAAACCTGGCCGCGCTGTTTGTGCCTCTCTAACCCCTGGCCTTCCTTTTTCCTTTTATCGAACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCTAAcc                                                                                     <  1:398838/131‑1 (MQ=255)
                                                                  aaGATTAGGAGTTGTTTTTCACTGGCCCGGAACTGACCGGGGGCGCGGTGGTTGTGGCTCTCGGACACGGGGGCGTCCTGTTTACCTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGt                                                                                  <  1:33462/148‑1 (MQ=255)
                                                                  aaGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                               >  1:746385/1‑151 (MQ=255)
                                                                  aaGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                               <  2:311573/151‑1 (MQ=255)
                                                                   aGATTACGTGTTGTTTTTCACTGTGCCGAAACTGACCCTGGGCGCGCTGGTTTTGGCTCCCGCCCACCTGGGCGTCCCGTTTCCCTTTACCGGGCAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACa                                                                              <  1:70085/151‑1 (MQ=255)
                                                                        aCGGGTTGTTTTTCACTGTGCCGGACCTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                          <  2:196862/150‑1 (MQ=255)
                                                                             ttttttttCCCTGGGCGGAAACTAACCCGTGGCGCGCTGGTTGTGGCTCTCGGACCCCGGGGCGTACCTTTTACCTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                      <  2:316560/144‑1 (MQ=255)
                                                                             ttGTTTTTCACTGTGCGGGAACTGAACGGGGGGGCGGGGGATGGGGCTCTCGGACACGGGGGCGTCCCGTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                      <  2:34951/149‑1 (MQ=255)
                                                                                   ttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                              >  2:876111/1‑151 (MQ=255)
                                                                                        tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                         >  2:791730/1‑151 (MQ=255)
                                                                                        ggggcggaaaCTAACCCGTGGGGCGTGGGTTGTGCCTCTCGCACCCGTGGGCGTCCTTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa                                                          <  1:634341/142‑1 (MQ=255)
                                                                                                 aCTGACCCGGGGCGCGCTGGATGTGGCTCTGGGACCCGTGGGCGTACCGTTTCCCTGTATCGGACAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                <  2:157087/151‑1 (MQ=255)
                                                                                                   tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                              >  1:40525/1‑151 (MQ=255)
                                                                                                     acccGGGGCGGGTTGGATGTGGCTCTCGGCCCCCGGGGCGTCCCGTTTCCGTGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                             <  1:249130/148‑1 (MQ=255)
                                                                                                     aaCCGTGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCTTTTCCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                               <  2:860744/148‑1 (MQ=255)
                                                                                                      accgggggCGCGCTGTGTGTTGTTCGCGGACCCGGGGGGTTACTTTTTCCTTGTAGCGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                              <  1:673195/143‑1 (MQ=255)
                                                                                                       ccGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTACCGTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                          <  2:76401/151‑1 (MQ=255)
                                                                                                            ggggcgcGGGTTGGGGCTCTGGGCCCCGGGGGCGTCCGTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                       <  2:631068/147‑1 (MQ=255)
                                                                                                                cgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                  >  2:239527/1‑151 (MQ=255)
                                                                                                                     gtggggcctctGGGCCACCGGGGCGTGCGTTTTCCCGGTAGGGGGAAATTGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg                             <  1:562773/144‑1 (MQ=255)
                                                                                                                                cGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGCCCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                     >  1:225548/1‑148 (MQ=255)
                                                                                                                                 ggcccccTGGGCGTCCCTTTTCCTTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                     <  2:464265/142‑1 (MQ=255)
                                                                                                                                  accacgggggCGTACGTTTTCCCTTTAGCGGAAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                 <  1:784553/143‑1 (MQ=255)
                                                                                                                                      ctggggaaacctttttCCCTGTAGGGGAAAAAGGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAg              <  1:987633/136‑1 (MQ=255)
                                                                                                                                         gggCGTCCCTTTTCCTTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa               <  2:667289/145‑1 (MQ=255)
                                                                                                                                               aCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc   >  1:382203/1‑151 (MQ=255)
                                                                                                                                                ccGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg  >  2:692792/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TTTCTGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCG  >  NC_000913/436355‑436648

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