Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 268,746 G→C 9.5% intergenic (‑77/+23) PP_0216 ← / ← desA sensory box/GGDEF family protein/delta‑9 fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947268,7460GC 9.5% 128.4 / 5.4 42intergenic (‑77/+23)PP_0216/desAsensory box/GGDEF family protein/delta‑9 fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base C (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTT  >  NC_002947/268599‑268894
                                                                                                                                                   |                                                                                                                                                    
gggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgggg                                                                                                                                                    <  1:67999/150‑1 (MQ=255)
           tGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCCCAGCGGGGCGTGcca                                                                                                                                           <  2:528135/148‑2 (MQ=255)
           tGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCCCAGCGGGGCGTGcca                                                                                                                                           >  1:528135/1‑147 (MQ=255)
           tGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAACGGCACGCCCCCCAGCGGGGCGTGccatt                                                                                                                                         <  2:1043327/150‑4 (MQ=255)
                  ttCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggc                                                                                                                                  >  1:1415741/1‑150 (MQ=255)
                   tCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                                 <  2:1041972/150‑1 (MQ=255)
                        gCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAg                                                                                                                            <  1:163895/150‑1 (MQ=255)
                         ccATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAgg                                                                                                                           <  2:324442/150‑1 (MQ=255)
                              ggCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCg                                                                                                                      <  2:1387839/149‑1 (MQ=255)
                                 tatgAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGCTTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCCCGCACCGCTGGGGGGGGTGCCGATCAACGGCGGTCAGGCGCCGGtg                                                                                                                   >  1:410219/4‑149 (MQ=255)
                                          gTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                                 <  2:799857/127‑1 (MQ=255)
                                          gTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                                 >  1:799857/1‑127 (MQ=255)
                                             tcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCCGCGAg                                                                                                       >  1:596379/1‑150 (MQ=255)
                                             tcAACGCGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAg                                                                                                       >  2:424856/1‑150 (MQ=255)
                                               aaCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggggg                                                                                                                                                   >  1:295301/1‑104 (MQ=255)
                                               aaCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggggg                                                                                                                                                   <  2:295301/104‑1 (MQ=255)
                                               aaCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTa                                                                                                     >  1:573692/1‑150 (MQ=255)
                                                aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                    <  1:355812/150‑1 (MQ=255)
                                                  gTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGgcgc                                                                                                                        <  1:953869/128‑1 (MQ=255)
                                                  gTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGgcgc                                                                                                                        >  2:953869/1‑128 (MQ=255)
                                                      aaGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTACGGATTAGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCGGGTCGGGCGCAGCGAGTAAGtagtg                                                                                              >  2:480168/1‑146 (MQ=255)
                                                       aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                             <  1:533155/150‑1 (MQ=255)
                                                        gCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCGGGTCGGGCGCGGGGAGTAAGTCTTGAg                                                                                            >  1:1252064/1‑150 (MQ=255)
                                                          gTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGGGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAGGGCGGTCAGGGGCCGGTCGGGCGCAGCGAGTAAGTCTTGAggt                                                                                          >  1:1176048/1‑148 (MQ=255)
                                                            tcGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTtgt                                                                                        <  2:423098/149‑1 (MQ=255)
                                                              ggCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGc                                                                                                                                             <  1:213850/95‑1 (MQ=255)
                                                              ggCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGc                                                                                                                                             >  2:213850/1‑95 (MQ=255)
                                                                   gCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGa                                                                                 >  2:674561/1‑150 (MQ=255)
                                                                      tCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggg                                                                                                                                                     >  1:217104/1‑79 (MQ=255)
                                                                      tCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggg                                                                                                                                                     <  2:217104/79‑1 (MQ=255)
                                                                                       gCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTc                                                                                                                   <  1:1167062/96‑1 (MQ=255)
                                                                                       gCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTc                                                                                                                   >  2:1167062/1‑96 (MQ=255)
                                                                                          atatTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgg                                                                                                                                                      <  2:1136881/58‑1 (MQ=255)
                                                                                          atatTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgg                                                                                                                                                      >  1:1136881/1‑58 (MQ=255)
                                                                                              tCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTgg                                                      >  2:269497/1‑150 (MQ=255)
                                                                                                cGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGtt                                                    >  2:961914/1‑150 (MQ=255)
                                                                                                  gATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTc                                                  <  1:269497/150‑1 (MQ=255)
                                                                                                                    ggggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTc                                <  1:1317588/150‑1 (MQ=255)
                                                                                                                         gAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCa                           <  2:573692/150‑1 (MQ=255)
                                                                                                                            gctCCCAATGGCCCGCGCCGTTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATgg                        <  2:1149025/148‑1 (MQ=255)
                                                                                                                            gTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATgg                        >  2:422710/1‑150 (MQ=255)
                                                                                                                            gTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATgg                        >  1:553382/1‑150 (MQ=255)
                                                                                                                                  aaTGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGGGGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCATCATGTATCCAGTCTTTCATGGCTGCCa                  >  1:1022793/1‑150 (MQ=255)
                                                                                                                                    tGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGc                <  1:961914/150‑1 (MQ=255)
                                                                                                                                      gCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTc                                                                                                 <  2:1121881/67‑1 (MQ=255)
                                                                                                                                      gCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTc                                                                                                 >  1:1121881/1‑67 (MQ=255)
                                                                                                                                       cACGCCCCCCAGCGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATg             >  2:784948/1‑150 (MQ=255)
                                                                                                                                        aCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTg                                                                 >  1:635187/1‑97 (MQ=255)
                                                                                                                                        aCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTg                                                                 <  2:635187/97‑1 (MQ=255)
                                                                                                                                                  ggggggCGTGCCGATCAACGGCGGTCAGGCGCCGTTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCAt                                          <  1:563596/110‑1 (MQ=255)
                                                                                                                                                  ggggggCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGTGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGtt  <  2:177326/150‑1 (MQ=255)
                                                                                                                                                  ggggggCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCAt                                          >  2:563596/1‑110 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTT  >  NC_002947/268599‑268894

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