Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 99,380 G→T A326A (GCG→GCT ftsW → putative lipid II flippase FtsW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091399,3800GT98.1% 185.1 / ‑3.0 53A326A (GCG→GCTftsWputative lipid II flippase FtsW
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (29/23);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCGATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGACCAAAGGT  >  NC_000913/99248‑99515
                                                                                                                                    |                                                                                                                                       
tAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGa                                                                                                                                  <  1:758932/140‑1 (MQ=255)
 aaCTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGAt                                                                                                                                 <  1:21542/140‑1 (MQ=255)
        tACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGt                                                                                                                          >  1:1136182/1‑140 (MQ=255)
         aCAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTa                                                                                                                         <  1:1183044/140‑1 (MQ=255)
         aCAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTa                                                                                                                         <  1:962701/140‑1 (MQ=255)
         aCAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTa                                                                                                                         <  2:12916/140‑1 (MQ=255)
         aCAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTa                                                                                                                         <  2:1208764/140‑1 (MQ=255)
         aCAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTa                                                                                                                         <  1:708170/140‑1 (MQ=255)
          cAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTaa                                                                                                                        <  2:634151/140‑1 (MQ=255)
                  ggAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAg                                                                                                                >  1:949999/1‑140 (MQ=255)
                        tCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTg                                                                                                          >  1:778142/1‑140 (MQ=255)
                                aaGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGt                                                                                                  >  2:531240/1‑140 (MQ=255)
                                     cacaCTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTc                                                                                             >  1:320740/1‑140 (MQ=255)
                                       cacTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCg                                                                                           <  1:671283/140‑1 (MQ=255)
                                             cTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTtc                                                                                     >  1:351706/1‑140 (MQ=255)
                                                 aTTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGc                                                                                 <  1:136967/140‑1 (MQ=255)
                                                 aTTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATCCTTCGTCGCTTTTCGCGCGATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGc                                                                                 <  2:1324040/140‑1 (MQ=255)
                                                  ttttCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGcc                                                                                >  2:1228692/1‑140 (MQ=255)
                                                      cGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGtt                                                                            >  1:1108332/1‑140 (MQ=255)
                                                      cGCCATTATCGGCGAAGAACTGGGGAATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGATTTTCCGGTTTTCTCGCCTGtt                                                                            >  1:1415277/1‑140 (MQ=255)
                                                          aTTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTAt                                                                        <  2:1136182/140‑1 (MQ=255)
                                                               cGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCa                                                                   <  2:98021/140‑1 (MQ=255)
                                                               cGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCa                                                                   <  1:352940/140‑1 (MQ=255)
                                                                 gCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATc                                                                 <  1:965653/140‑1 (MQ=255)
                                                                       aaCTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTa                                                                         <  2:279245/126‑1 (MQ=255)
                                                                       aaCTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTa                                                                         >  1:279245/1‑126 (MQ=255)
                                                                        aCTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTa                                                          >  1:1357226/1‑140 (MQ=255)
                                                                            gggTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCtt                                                      >  1:62764/1‑140 (MQ=255)
                                                                                 tGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCACTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAgg                                                 >  1:1272319/1‑140 (MQ=255)
                                                                                  gTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGgc                                                <  2:1357226/140‑1 (MQ=255)
                                                                                      gtgtGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGgcg                                               >  1:1125660/1‑137 (MQ=255)
                                                                                      gtgtGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGgcg                                               <  2:1125660/137‑1 (MQ=255)
                                                                                           gtgCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTaa                                       >  2:978077/1‑140 (MQ=255)
                                                                                            tgCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAAc                                      >  1:65041/1‑140 (MQ=255)
                                                                                            tgCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAAc                                      >  2:526090/1‑140 (MQ=255)
                                                                                                      tttAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGtt                                         <  1:1333633/127‑1 (MQ=255)
                                                                                                      tttAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGtt                                         >  2:1333633/1‑127 (MQ=255)
                                                                                                       ttAATGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGTGCTGGTTAACGTAGGCgcggc                           <  2:320740/140‑1 (MQ=255)
                                                                                                           tGGTATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCgcggcggcg                       >  2:203067/1‑140 (MQ=255)
                                                                                                              tATTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCgggg                    <  2:1415277/140‑1 (MQ=255)
                                                                                                               aTTCTTCGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGa                   >  2:1132315/1‑140 (MQ=255)
                                                                                                                   ttcGTCGCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCAGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGtt               >  1:1239251/1‑140 (MQ=255)
                                                                                                                     cgtcgCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTAc             >  1:1172957/1‑140 (MQ=255)
                                                                                                                      gtcgCTTTTCGCACTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTAcc            <  2:1272319/140‑1 (MQ=255)
                                                                                                                        cgCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTa              <  1:776932/136‑1 (MQ=255)
                                                                                                                        cgCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTa              >  2:776932/1‑136 (MQ=255)
                                                                                                                        cgCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGa          >  1:494796/1‑140 (MQ=255)
                                                                                                                        cgCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGa          >  1:373914/1‑140 (MQ=255)
                                                                                                                        cgCTTTTCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGa          >  1:1187634/1‑140 (MQ=255)
                                                                                                                           ttttCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGACCa       <  2:789404/140‑1 (MQ=255)
                                                                                                                             ttCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGACCaaa     >  2:487948/1‑140 (MQ=255)
                                                                                                                             ttCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGACCaaa     >  2:945742/1‑140 (MQ=255)
                                                                                                                             ttCGCGCTATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATCGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGACCaaa     >  2:1172095/1‑140 (MQ=255)
                                                                                                                                gcgctATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGACCAAAGGt  <  1:531240/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TAACTCGGTACAAAAACTGGAGTATCTGCCGGAAGCGCACACTGACTTTATTTTCGCCATTATCGGCGAAGAACTGGGGTATGTCGGTGTGGTGCTGGCACTTTTAATGGTATTCTTCGTCGCTTTTCGCGCGATGTCGATTGGCCGTAAAGCATTAGAAATTGACCACCGTTTTTCCGGTTTTCTCGCCTGTTCTATTGGCATCTGGTTTAGCTTCCAGGCGCTGGTTAACGTAGGCGCGGCGGCGGGGATGTTACCGACCAAAGGT  >  NC_000913/99248‑99515

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