Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,565,696 (A)7→6 coding (38/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,6900A.96.5% 110.7 / ‑2.8 29coding (44/2400 nt)dosPoxygen‑sensing c‑di‑GMP phosphodiesterase DosP
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (13/15);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGCATATCAATGTTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCAT  >  NC_000913/1565552‑1565825
                                                                                                                                          |                                                                                                                                       
aTCAGCATATCAATGTTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGaa                                                                                                                                        <  1:98396/140‑1 (MQ=255)
   aGCATATCAATGTTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAGCACCGCACCCATCATATTTTGCTCAAGGGCGGGGaaaaa                                                                                                                                     <  2:507487/140‑1 (MQ=255)
         tCAATGTTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCa                                                                                                                              >  1:509622/1‑140 (MQ=255)
         tCAATGTTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCAACCATCATATTTAGCGCAAGGGCGGGGAAAAAATGCCa                                                                                                                              >  1:225819/1‑140 (MQ=255)
               ttATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATcggcg                                                                                                                        >  1:556260/1‑140 (MQ=255)
               ttATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAAAATCACTTCATCATTTTTGTTAATTAAGACCTCACCCATCATATTTTTTTCCAAGGCGGGGAAAAAAAGCCATcggcg                                                                                                                        >  2:156369/1‑140 (MQ=255)
                    gCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAAAATCACTTCATCATTTTTGTTAATTAACACCGCACCCATCATATTTTTCTCAAAGGCGGGGAAAAAAAGCCATCGGCGGCAtt                                                                                                                   >  2:130184/1‑140 (MQ=255)
                       aaTGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGGTGAAAAACATCACTTCATCATTTTTGTTAATTAACACCGCACCCATCATATTTTTCTCAAAGGCGGGGGAAAAATACCATCGGGGGCAttttt                                                                                                                >  2:142568/1‑137 (MQ=255)
                         tGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGGTGGAAAACATCACTTCCTCATTTTCGTTAATTAACACCGCACCCCTCATATTTTTCTCAAAGGCGGGGGAAAAATGTCATCTGCGGGATTATccc                                                                                                              >  2:599658/1‑139 (MQ=255)
                         tGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGGTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCg                                                                                                              >  1:624124/1‑140 (MQ=255)
                                 tcACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGtt                                                                                                      <  2:47941/140‑1 (MQ=255)
                                     gTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCt                                                                                                  <  2:627245/140‑1 (MQ=255)
                                       ttGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTc                                                                                                <  1:130184/140‑1 (MQ=255)
                                          tATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAAAATCACTTCATCATTTTCGTTAGTTTACACCGCACCCCTCATATTTTGCTCCAAGGGGGGGGAAAGATGCCTTAGGCGGGATTTTCCCACTTGGTTAGCTTccta                                                                                             >  2:312126/1‑137 (MQ=255)
                                                   agCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACctct                                                                                    <  2:624124/140‑1 (MQ=255)
                                                   agCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACctct                                                                                    <  2:505676/140‑1 (MQ=255)
                                                            gcggcgGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTgg                                                                           <  1:115317/140‑1 (MQ=255)
                                                            gcggcgGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTgg                                                                           <  2:556260/140‑1 (MQ=255)
                                                            gcggcgGGGTTCAAACACATCATTTCATCTTTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTgg                                                                           <  1:174221/140‑1 (MQ=255)
                                                             cggcggGGTTTAAAAAAAACACTTCATTATTTTCGTTTATTAACACCGCACCCATCATATTTTTCTCAAAGGCGGGGGAAAAAAGCCATCCGCGGCATTTTCCCCATCGGGTAGCTTTATTATTACCCCTGCATCCCggg                                                                          >  2:487381/1‑139 (MQ=255)
                                                              ggcggGGTTGAAAAACATCACTGCAGCATTTTCGTTACTTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGgcg                                                                         <  1:312126/140‑1 (MQ=255)
                                                                             cATCGCTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGc                                                          >  2:439770/1‑140 (MQ=255)
                                                                                   ttCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCa                                                    <  2:99967/140‑1 (MQ=255)
                                                                                        catTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGtt                                               >  1:63126/1‑140 (MQ=255)
                                                                                             tCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAaca                                          <  2:100615/140‑1 (MQ=255)
                                                                                             tCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAaca                                          >  1:503989/1‑140 (MQ=255)
                                                                                                             cACCCATCATATTTTGCTCAAGGGCGGGGACAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTc                          >  2:613238/1‑140 (MQ=255)
                                                                                                             cACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAATGACAGGCTTTCCAGAGATCAACACGGGTTCTACCTCGTc                          >  2:253780/1‑140 (MQ=255)
                                                                                                             cACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGCCTGGCTTACCATAGCTCAACATGGTTTCTACCTCGTc                          >  1:308492/1‑140 (MQ=255)
                                                                                                             cACCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTc                          >  1:423781/1‑140 (MQ=255)
                                                                                                              aCCCATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCt                         >  2:644495/1‑140 (MQ=255)
                                                                                                                ccATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt                       <  2:509622/140‑1 (MQ=255)
                                                                                                                ccATCATATTTTGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt                       <  1:142568/140‑1 (MQ=255)
                                                                                                                     atatTTTGCTCAAGGGCGGGGAATAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                  <  1:506571/140‑1 (MQ=255)
                                                                                                                          ttGCTCAAGGGCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGAtata             <  1:644495/140‑1 (MQ=255)
                                                                                                                                  gggCGGGGAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCGAGAGTTCATCACGGTTACTACCTCGTCTTTTGGCGATATACAGAGCtt     >  2:195592/1‑140 (MQ=255)
                                                                                                                                      cggggAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCTTGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCAt  >  2:334011/2‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                       
ATCAGCATATCAATGTTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCAT  >  NC_000913/1565552‑1565825

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