Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,162,099 T→C 8.3% 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 8.3% 131.2 / 2.4 48intergenic (‑43/+34)PP_3663/pssAGGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base C (2/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

AAGGCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGA  >  NC_002947/4161952‑4162245
                                                                                                                                                   |                                                                                                                                                      
aaGGCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTg                                                                                                                                                       >  1:299400/1‑149 (MQ=255)
 aGGCCCGAATGTATGCGATAGTTATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTgg                                                                                                                                                      <  8:154821/149‑1 (MQ=255)
   gCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCt                                                                                                                                                    >  3:131733/1‑149 (MQ=255)
   gCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGc                                                                                                                                                     >  7:251860/1‑148 (MQ=255)
            tATGCGATCGTGATCACGATCAGAGGGGACTTGCACTGTTGCCAGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCACCAATATCATAATGACATCAATGCGCGGTGCGCACCACTATGGCCATGTCCATGGCTTGCCGGCGa                                                                                                                                           >  7:274751/1‑149 (MQ=255)
               gCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAg                                                                                                                                        >  3:246114/1‑149 (MQ=255)
               gCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACa                                                                                                                                         <  8:251860/148‑1 (MQ=255)
                cGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCCGCCCCGGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAgg                                                                                                                                       <  3:338923/149‑1 (MQ=255)
                cGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACAg                                                                                                                                        <  2:260928/148‑1 (MQ=255)
                       gATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTa                                                                                                                                >  1:413/1‑149 (MQ=255)
                         tCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTa                                                                                                                                >  7:70681/1‑147 (MQ=255)
                                cAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGAAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGAt                                                                                                                        <  7:144435/148‑1 (MQ=255)
                                 agagGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTCGGCAGCAATATCAAAATGACATCAATGCGCGGTGCGCACCACTCTGCCCCTGTCGCTGGCTTGCCGGCGACAGCGTCAGTACAGGCGATga                                                                                                                      >  4:36733/1‑148 (MQ=255)
                                  gagGGGCCTTGCACTTTTGCCTTCACATTTTAATCCAGAACCTGACAAAGGGCCTTGCGCAGCAATATCATAATGACATCAATGCGCGGAGGGCAACACTCTTGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTc                                                                                                                     <  5:273404/149‑1 (MQ=255)
                                    ggggCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGt                                                                                                                      >  8:178107/1‑146 (MQ=255)
                                               cTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCtt                                                                                                       <  3:31662/149‑1 (MQ=255)
                                                  ttGCCTGCACAATGTAATCAAGCCGCTGAGAATGGGCCGTGGGCAGCAATATCATATTGACATCAATGCGGGGTTCGCACCACTCTGGCCCTGTCGGTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGGACGC‑GCTTGa                                                                                                     <  8:174494/148‑1 (MQ=255)
                                                      cTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAAc                                                                                                >  4:92532/1‑149 (MQ=255)
                                                          acaTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCt                                                                                            <  3:176111/149‑1 (MQ=255)
                                                                  aTCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCgg                                                                                  <  8:70681/149‑1 (MQ=255)
                                                                           cTGACAAAGGGCCGTGGGCAGCCATATAATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGCGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAgcgc                                                                         >  7:237044/1‑149 (MQ=255)
                                                                             gACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCaa                                                                       >  3:303744/1‑149 (MQ=255)
                                                                                aaaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑cc                                                                    >  2:362884/1‑148 (MQ=255)
                                                                                 aaGGGCCGTGGTCAGCAATATCATAATGACATCTATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGa                                                                  <  6:312189/149‑1 (MQ=255)
                                                                                 aaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCTGCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCg                                                                   <  1:257312/149‑1 (MQ=255)
                                                                                 aaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGa                                                                  <  2:299400/149‑1 (MQ=255)
                                                                                     gCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCa                                                              <  2:136942/149‑1 (MQ=255)
                                                                                                    atCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGACCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCActt                                               >  2:170716/1‑149 (MQ=255)
                                                                                                          aTGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCActtgcttg                                          <  4:131733/148‑1 (MQ=255)
                                                                                                           tGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTgg                                         >  3:97629/1‑148 (MQ=255)
                                                                                                               aTCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGt                                    >  8:285248/1‑149 (MQ=255)
                                                                                                                 cAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCtata                                                                             >  1:188279/1‑107 (MQ=255)
                                                                                                                  aaTGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGTCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGATGCAGAAACTTGCGCACTTGCTTGGGATGTTCg                                 <  8:83298/149‑1 (MQ=255)
                                                                                                                  aaTGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTc                                  <  4:97629/148‑1 (MQ=255)
                                                                                                                     gcgcGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGcc                              >  1:355725/1‑149 (MQ=255)
                                                                                                                      cgcgGTGCGCACCACTCTTGCCCTGTCGCTGTCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATGAACGGCTCGA‑‑ACCGGGTATAGCGGAA‑GCGTCTCAGAAACTTGCTCACTTGCTTGGGATGTTCGGcc                              <  8:274751/148‑1 (MQ=255)
                                                                                                                         ggtcggcACCAGTCTGTCCCTGTCGCTGGCTTGCCTGCGACAGTGTCATTACAGGCGATGTCAGAGCACCC‑GCTTGTTCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACTCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTtct                          <  4:246114/144‑1 (MQ=255)
                                                                                                                         ggggcgcACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTtct                          <  1:206046/146‑1 (MQ=255)
                                                                                                                             cgcACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTccc                       <  4:303744/148‑1 (MQ=255)
                                                                                                                               caccacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTcccccc                    >  7:296137/1‑148 (MQ=255)
                                                                                                                                accacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAg                   >  5:194366/1‑149 (MQ=255)
                                                                                                                                 ccacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAgc                  <  7:178107/149‑1 (MQ=255)
                                                                                                                                  cacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAgca                 >  1:323941/1‑148 (MQ=255)
                                                                                                                                         gCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGc          <  7:285248/149‑1 (MQ=255)
                                                                                                                                          cccTGGCGCTGGCTTGCCGGCGTCAGGGGCAGTACACGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGGCTTCTCCCCCAGCGCCTCTGcc         <  1:367364/149‑1 (MQ=255)
                                                                                                                                           ccTGTCGCCGGCAAGCCAGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCATGGTGGGCTATAGCGCAATGCCACCCAGAAACATGCGCACTTGCTTGGGATTTTGGGCCTTATCCCCCAGCGccaatg           >  1:118546/1‑145 (MQ=255)
                                                                                                                                             ggTCGCTGGCTTTCCGGCGACATGGTCAGTTCATGCGATGTCATACCAACC‑GCTTTATCAACGGCTTCA‑‑AACCGCTATACCGCAA‑CCGAAGCATAAACTTTCGCACTTTCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGaa      <  8:296137/148‑1 (MQ=255)
                                                                                                                                                 gCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGa  >  7:77360/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAGGCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCC‑GCTTGATCAACGGCTCCA‑‑ACCGGCTATAGCGCAA‑CCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGA  >  NC_002947/4161952‑4162245

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