Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,162,099 T→C 7.5% intergenic (‑43/+34) PP_3663 ← / ← pssA GGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,162,0990TC 7.5% 147.6 / 2.1 53intergenic (‑43/+34)PP_3663/pssAGGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base T (29/20);  new base C (2/2);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAA  >  NC_002947/4161957‑4162241
                                                                                                                                              |                                                                                                                                              
ccGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTg                                                                                                                                          >  5:340862/1‑149 (MQ=255)
  gACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTgg                                                                                                                                              <  1:48842/143‑1 (MQ=255)
  gACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTgg                                                                                                                                              >  2:48842/1‑143 (MQ=255)
   aCTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCg                                                                                                                                       >  6:12146/1‑149 (MQ=255)
    cTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGcaagccag                                                                                                                                      >  2:105684/1‑142 (MQ=255)
        aTGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGAc                                                                                                                                  >  8:263227/1‑149 (MQ=255)
            gATCTTTATCACTTTCAGAGTGTCCTTGCACTGTTGCCTGAACATTGTAATCCAGCCCCTGACAAAGGGCCGGGTGCAGCAATATCATAATGACATCAATGCGCGGTCCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAgg                                                                                                                               <  8:125668/148‑1 (MQ=255)
             aTCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACAgggc                                                                                                                             <  1:105684/149‑2 (MQ=255)
              tCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGa                                                                                                                                   <  1:328010/142‑1 (MQ=255)
              tCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGa                                                                                                                                   >  2:328010/1‑142 (MQ=255)
                gTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCa                                                                                                                           <  1:128232/148‑1 (MQ=255)
                     cACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCAAAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAg                                                                                                                     <  7:263227/149‑1 (MQ=255)
                      aCGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAgg                                                                                                                    <  6:340862/149‑1 (MQ=255)
                        gTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGc                                                                                                                   >  3:229590/1‑148 (MQ=255)
                                gggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                          <  3:207245/149‑1 (MQ=255)
                                     ttGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCAcc                                                                                                     >  1:113317/1‑149 (MQ=255)
                                           tGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCag                                                                                                           >  7:216892/1‑137 (MQ=255)
                                           tGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCag                                                                                                           <  8:216892/137‑1 (MQ=255)
                                                   gaacaTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGACTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACgg                                                                                       >  4:76700/3‑149 (MQ=255)
                                                        ttGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTcc                                                                                   >  5:270916/1‑148 (MQ=255)
                                                              tCCAGCCCCTGACAAAGGGGCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGc                                                                            <  6:270916/149‑1 (MQ=255)
                                                                cAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAt                                                                                                                >  5:151260/1‑111 (MQ=255)
                                                                cAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAt                                                                                                                <  6:151260/111‑1 (MQ=255)
                                                                   ccccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAg                                                                       >  5:228762/1‑149 (MQ=255)
                                                                    cccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCta                                                                          <  3:78109/145‑1 (MQ=255)
                                                                    cccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCta                                                                          >  4:78109/1‑145 (MQ=255)
                                                                     ccTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCGAACCGGCTATAgcg                                                                     >  6:79832/1‑149 (MQ=255)
                                                                       tGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAgcgc                                                                    >  2:129819/1‑148 (MQ=255)
                                                                           aaaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAAcc                                                                >  6:230504/1‑148 (MQ=255)
                                                                               ggCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGc                                                           <  1:299123/149‑1 (MQ=255)
                                                                               ggCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAAGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACg                                                            >  1:279171/1‑148 (MQ=255)
                                                                                 ccGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACg                                                            >  3:28797/1‑146 (MQ=255)
                                                                                    tGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGaaa                                                      <  4:28797/149‑1 (MQ=255)
                                                                                       gcagcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTc                                                                                    <  2:152480/116‑1 (MQ=255)
                                                                                        cagcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTg                                                  >  5:267452/1‑149 (MQ=255)
                                                                                         agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgc                                                 >  3:47293/1‑149 (MQ=255)
                                                                                          gcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgcg                                                <  4:229590/149‑1 (MQ=255)
                                                                                           caATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgcg                                                >  8:175151/1‑148 (MQ=255)
                                                                                              tatCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCAct                                            >  1:188562/1‑149 (MQ=255)
                                                                                                tCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActt                                           <  7:175151/148‑1 (MQ=255)
                                                                                                tCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCAATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTAAAGCGCACCCAACGAACAAACTTACCCAcat                                           >  6:235751/1‑146 (MQ=255)
                                                                                                      tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCGTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTgg                                     >  1:70439/1‑148 (MQ=255)
                                                                                                          aTCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGt                                <  4:47293/149‑1 (MQ=255)
                                                                                                            cAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTc                              >  2:305013/1‑149 (MQ=255)
                                                                                                               tGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGc                           <  5:230504/149‑1 (MQ=255)
                                                                                                                  gcgGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCtt                        >  4:118763/1‑149 (MQ=255)
                                                                                                                      tGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTcc                    <  2:188562/149‑1 (MQ=255)
                                                                                                                      tGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTcc                    >  6:303695/1‑149 (MQ=255)
                                                                                                                             cacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTGAGAGCCCCCGCTTGCTCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGGACTTGCTTGGGCTGTTCGGCCTTCTCGCCCAgcg             >  7:287492/1‑149 (MQ=255)
                                                                                                                               ctctGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGcc           >  5:307584/1‑149 (MQ=255)
                                                                                                                                    gCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCctct        <  1:168850/147‑1 (MQ=255)
                                                                                                                                     cccTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGcc     <  6:228762/149‑1 (MQ=255)
                                                                                                                                        tGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGaa  >  6:21829/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAA  >  NC_002947/4161957‑4162241

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