Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,162,143 T→C 12.8% 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,1430TC12.8% 137.5 / 9.0 55intergenic (‑43/+34)PP_3663/pssAGGDEF domain‑containing protein/CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base C (4/3);  total (29/26)
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.94e-01

AAGGCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGAT  >  NC_002947/4161996‑4162290
                                                                                                                                                   |                                                                                                                                                   
aaGGCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTgg                                                                                                                                                   <  2:66870/150‑1 (MQ=255)
aaGGCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTgg                                                                                                                                                   <  1:601275/150‑1 (MQ=255)
  ggCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGc                                                                                                                                                  >  2:863287/1‑149 (MQ=255)
  ggCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGc                                                                                                                                                  <  1:863287/149‑1 (MQ=255)
   gCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCtt                                                                                                                                                >  1:226538/1‑150 (MQ=255)
     ccGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGc                                                                                                                                              <  1:974272/150‑1 (MQ=255)
     ccGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGc                                                                                                                                              >  1:377126/1‑150 (MQ=255)
             aTGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACa                                                                                                                                      <  1:987646/150‑1 (MQ=255)
             aTGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACa                                                                                                                                      <  1:1201810/150‑1 (MQ=255)
              tGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAg                                                                                                                                     <  1:124878/150‑1 (MQ=255)
                cGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAggg                                                                                                                                   <  2:226538/150‑1 (MQ=255)
                   tCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAg                                                                                                                                     >  1:1009948/1‑145 (MQ=255)
                   tCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAg                                                                                                                                     <  2:1009948/145‑1 (MQ=255)
                     gTGATCACGTTCAGAGGGGCCTTGCGCTGTTGCCTGGAGATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATGTCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCGGGCTTGCCGGCTACAGGGtgcgt                                                                                                                              >  2:407606/1‑146 (MQ=255)
                                        ccTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCAc                                                                                                           >  1:999228/1‑150 (MQ=255)
                                            gCACTGTTGCCTGCACATTGTAATAAAGATCCGGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGc                                                                                                       >  1:37563/1‑150 (MQ=255)
                                                 agatctggCACATTGCAATCAGGTCGCTGACAAAGTGCCTTTGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGAt                                                                                                  <  2:1373404/143‑1 (MQ=255)
                                                  ttGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATc                                                                                                 >  2:1016737/1‑150 (MQ=255)
                                                    gCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCaa                                                                                               >  1:742885/1‑150 (MQ=255)
                                                     ccTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACAg                                                                                                                                     <  2:504548/111‑1 (MQ=255)
                                                     ccTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCCGGCAAGCCAGCGACAg                                                                                                                                     >  1:504548/1‑111 (MQ=255)
                                                        gCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCagatcg                                                                                           >  1:1373404/1‑145 (MQ=255)
                                                          acaTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTc                                                                                         <  1:625133/150‑1 (MQ=255)
                                                            aTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGAt                                                                                                  <  2:1225504/139‑1 (MQ=255)
                                                            aTTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGAt                                                                                                  >  1:1225504/1‑139 (MQ=255)
                                                             ttGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                  >  1:856904/1‑106 (MQ=255)
                                                             ttGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                  <  2:856904/106‑1 (MQ=255)
                                                                       gCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAg                                                                            <  2:802667/150‑1 (MQ=255)
                                                                            tGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCaa                                                                       >  1:1115746/1‑150 (MQ=255)
                                                                                aaaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                  >  1:1188336/1‑87 (MQ=255)
                                                                                aaaGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGt                                                                                                                                  <  2:1188336/87‑1 (MQ=255)
                                                                                      ccGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACa                                                                                                                                      <  1:512919/77‑1 (MQ=25)
                                                                                      ccGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGACCCTGTCGCCGGCAAGCCAGCGACa                                                                                                                                      >  2:512919/1‑77 (MQ=25)
                                                                                              agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCg                                                                                                                                            >  2:182374/1‑63 (MQ=255)
                                                                                              agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCg                                                                                                                                            <  1:182374/63‑1 (MQ=255)
                                                                                              agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCGGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgcg                                                     >  2:1390346/1‑150 (MQ=255)
                                                                                              agcaATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTgcg                                                     <  2:742885/150‑1 (MQ=255)
                                                                                                caATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAgg                                                                                                                                    <  2:1142414/69‑1 (MQ=255)
                                                                                                caATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAgg                                                                                                                                    >  1:1142414/1‑69 (MQ=255)
                                                                                                   tatCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                               <  2:697432/87‑1 (MQ=255)
                                                                                                   tatCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCaga                                                                                                               >  1:697432/1‑87 (MQ=255)
                                                                                                    atCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCActtg                                               <  2:377126/150‑1 (MQ=255)
                                                                                                               aTCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTAAGTCCAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGtt                                    >  2:194794/1‑150 (MQ=255)
                                                                                                                    tGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGcc                               <  2:999228/150‑1 (MQ=255)
                                                                                                                        cgGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTtct                           >  2:187641/1‑150 (MQ=255)
                                                                                                                                accacTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAgc                   <  1:1016737/150‑1 (MQ=255)
                                                                                                                                         gCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGcc          >  1:1406969/1‑150 (MQ=255)
                                                                                                                                         gCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGcc          >  1:1177582/1‑150 (MQ=255)
                                                                                                                                           ccTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTc                                                                                         <  2:1267108/69‑1 (MQ=255)
                                                                                                                                           ccTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTc                                                                                         >  1:1267108/1‑69 (MQ=255)
                                                                                                                                            cTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGaa       >  2:180877/1‑150 (MQ=255)
                                                                                                                                            cTGTCGCCGGCAAGCCAGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGACCTTCTCCCCCAGCGCCTCTGCCGaa       >  2:39623/1‑150 (MQ=255)
                                                                                                                                            cTGTCGCCGGCAAGCCAGCGACAGGGCCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGGGCCTCTGCCGaa       >  2:441288/1‑150 (MQ=255)
                                                                                                                                             tGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAAc      >  1:1159259/1‑150 (MQ=255)
                                                                                                                                                 gCTGGCTTGCCGGCGACAGGCTCCGTACAGGCGCTGTCAGAGCTCCCGCTTGAGCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGAt  <  1:1161509/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGGCCCGACTGTATGCGATCGTGATCACGTTCAGAGGGGCCTTGCACTGTTGCCTGCACATTGTAATCCAGCCCCTGACAAAGGGCCGTGGGCAGCAATATCATAATGACATCAATGCGCGGTGCGCACCACTCTGGCCCTGTCGCTGGCTTGCCGGCGACAGGGTCAGTACAGGCGATGTCAGAGCACCCGCTTGATCAACGGCTCCAACCGGCTATAGCGCAACCGACGCAGAAACTTGCGCACTTGCTTGGGATGTTCGGCCTTCTCCCCCAGCGCCTCTGCCGAACCGAT  >  NC_002947/4161996‑4162290

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