Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,565,689 Δ1 bp coding (45/2400 nt) dosP ← oxygen‑sensing c‑di‑GMP phosphodiesterase DosP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,565,6860G.96.5% 107.9 / ‑1.2 30coding (48/2400 nt)dosPoxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base . (15/13);  minor base T (0/1);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATA  >  NC_000913/1565567‑1565814
                                                                                                                       |                                                                                                                                
ttATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGGTGTAAAACATCACTTCATCATTTTTGTTAATTAACACCGCACCCATCATATTTTTCTCCAGGGCGGGaaaaaaaagcc                                                                                                                     >  2:131296/1‑130 (MQ=255)
  aTTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCAt                                                                                                                  >  2:53023/1‑133 (MQ=255)
     gCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATcgg                                                                                                               <  1:120950/133‑1 (MQ=255)
     gCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATcgg                                                                                                               >  1:531877/1‑133 (MQ=255)
      ccAATGACTTCTTCACGTTTGTATACCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATcggc                                                                                                              <  2:105821/133‑1 (MQ=255)
                cttcACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGGTGAAAAACATCCCTTCATCATTTTCGTTAATTAACACCGCACCCCTCATATTTTGCTCAAGAGCGAGAACAAAATGCCCTCGGCGGCATTATcc                                                                                                    >  2:232436/1‑133 (MQ=255)
                        ttGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGt                                                                                            <  2:531877/133‑1 (MQ=255)
                        ttGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTCCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGt                                                                                            <  1:375372/133‑1 (MQ=255)
                           tATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAAAATCACTTCATCATTTTTGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAg                                                                                         >  2:413632/1‑133 (MQ=255)
                             tCCCCAGATCTTCTCTGCGGCGGGGTTGAATAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCt                                                                                       <  1:245147/133‑1 (MQ=255)
                                  agagCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATg                                                                                  >  1:438719/1‑133 (MQ=255)
                                         ctctGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACct                                                                           >  1:595950/1‑133 (MQ=255)
                                         ctctGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAAGGCGGGAAAAAAATGCCATCGGCGGGCTTATCCCCATCGGTTAGCTTCATGATTACct                                                                           >  2:60023/1‑133 (MQ=255)
                                         ctctGCGGCGGGGGTGAAAAACATCACTTCCTCATTTTCGTTAATTAACACCGCACCCACCATATTTTGCTCAAAGGCGGGAAAAAAATGCACTCGGCGGAATTCTCCACCTCGGTTCGCTTTATTATTACct                                                                           >  2:243165/1‑133 (MQ=255)
                                           ctGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACctct                                                                         <  1:60023/133‑1 (MQ=255)
                                                      ttGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGgcg                                                              <  2:438719/133‑1 (MQ=255)
                                                                         catTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTc                                           >  1:118238/1‑133 (MQ=255)
                                                                                   aattaaCACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGAGCATCTAGAGAGGGTTTTCCATAGTTCAac                                >  1:200246/1‑133 (MQ=255)
                                                                                      taaCACCGCACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAacac                              <  1:53023/133‑1 (MQ=255)
                                                                                              cACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTa                      >  1:127000/1‑133 (MQ=255)
                                                                                              cACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTa                      >  2:462792/1‑133 (MQ=255)
                                                                                              cACCCATCATATTTTGCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTACCATAGTTCAACAAGGTTTCTa                      >  1:246953/1‑133 (MQ=255)
                                                                                                        atTTTGCTCAAGGGCTGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt            <  1:232436/133‑1 (MQ=255)
                                                                                                        atTTTGCTCAAGGGCAGGAAAACAATGCCACCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt            <  1:183951/133‑1 (MQ=255)
                                                                                                          tttGCTCAAGGGCGGAAAAAAAATGCCATCTGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg          <  1:243165/133‑1 (MQ=255)
                                                                                                             gCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg       >  2:147700/1‑133 (MQ=255)
                                                                                                             gCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg       <  2:612442/133‑1 (MQ=255)
                                                                                                             gCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg       >  2:73437/1‑133 (MQ=255)
                                                                                                             gCTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGCCTTTCCAGTGTTGAACGCGGTTTCTGCCTCGTCTTTTGGCg       >  2:183951/1‑133 (MQ=255)
                                                                                                              cTCAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGa      >  1:22243/1‑133 (MQ=255)
                                                                                                                cAAGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGata    <  1:413632/133‑1 (MQ=255)
                                                                                                                  aGGGCGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGAtata  <  2:22243/133‑1 (MQ=255)
                                                                                                                       |                                                                                                                                
TTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATA  >  NC_000913/1565567‑1565814

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: