Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 268,744 T→A 8.2% intergenic (‑75/+25) 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,7440TA 8.2% 138.6 / 4.6 49intergenic (‑75/+25)PP_0216/desAsensory box/GGDEF family protein/delta‑9 fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base A (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATG  >  NC_002947/268599‑268883
                                                                                                                                                 |                                                                                                                                           
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                                                                                                             >  2:953869/1‑128 (MQ=255)
                                                  gTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGgcgc                                                                                                             <  1:953869/128‑1 (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                                                                                                                                  >  2:213850/1‑95 (MQ=255)
                                                              ggCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGc                                                                                                                                  <  1:213850/95‑1 (MQ=255)
                                                                   gCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGa                                                                      >  2:674561/1‑150 (MQ=255)
                                                                      tCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggg                                                                                                                                          <  2:217104/79‑1 (MQ=255)
                                                                      tCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTggg                                                                                                                                          >  1:217104/1‑79 (MQ=255)
                                                                                       gCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTc                                                                                                        >  2:1167062/1‑96 (MQ=255)
                                                                                       gCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTc                                                                                                        <  1:1167062/96‑1 (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             >  1:553382/1‑150 (MQ=255)
                                                                                                                            gTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATgg             >  2:422710/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                                                      <  2:635187/97‑1 (MQ=255)
                                                                                                                                        aCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTg                                                      >  1:635187/1‑97 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATG  >  NC_002947/268599‑268883

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