Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,565,758 C→T *461* (TAG→TAA dosC ← diguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,565,7580CT97.4% 117.2 / ‑6.5 38*461* (TAG→TAAdosCdiguanylate cyclase, cold‑ and stationary phase‑induced oxygen‑dependent biofilm regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/19);  minor base G (0/1);  total (18/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 = 1.00e+00

ATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGA  >  NC_000913/1565638‑1565892
                                                                                                                        |                                                                                                                                      
atcatTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCag                                                                                                                      >  2:601683/1‑139 (MQ=255)
  catTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATGTAAAGACTGGCTTTCCagag                                                                                                                    <  2:225368/139‑1 (MQ=255)
         gTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAaca                                                                                                             <  1:273330/139‑1 (MQ=255)
                      gCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTc                                                                                                >  2:118306/1‑139 (MQ=255)
                       cACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCg                                                                                               >  2:799752/1‑139 (MQ=255)
                                 atcttcctccaaggCTGGGAAAAAAATGCCATCGGCGCCAATAGCGACATCGGTTAACTTCAGGATTACCTCTGCAGCCTGGGGTATTTAAAGAATGGCAGTCCAGAGTTCAACACGGTTTCTACCTCGTCCTTTGGCg                                                                                     <  1:102184/127‑1 (MQ=255)
                                 atTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg                                                                                     <  1:64463/139‑1 (MQ=255)
                                    ttGCTCAAGGGTGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGata                                                                                  >  2:469697/1‑139 (MQ=255)
                                      gCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGAtata                                                                                >  1:636110/1‑139 (MQ=255)
                                            gggCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGc                                                                          >  1:450418/1‑139 (MQ=255)
                                                 gggAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCAt                                                                     >  1:19648/1‑139 (MQ=255)
                                                           tGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCgg                                                                  >  1:231513/1‑132 (MQ=255)
                                                           tGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCgg                                                                  <  2:231513/132‑1 (MQ=255)
                                                                tCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACTTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTGTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGa                                                      <  2:452126/139‑1 (MQ=255)
                                                                  ggcggcATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAgg                                                    <  1:799752/139‑1 (MQ=255)
                                                                  ggcggcATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAgg                                                    <  2:390569/139‑1 (MQ=255)
                                                                        aTTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCa                                              <  2:19648/139‑1 (MQ=255)
                                                                                cATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAgg                                      >  1:673917/1‑139 (MQ=255)
                                                                                 aTCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                        >  1:89106/1‑88 (MQ=255)
                                                                                 aTCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg                                                                                        <  2:89106/88‑1 (MQ=255)
                                                                                 aTCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGa                                     >  1:236085/1‑139 (MQ=255)
                                                                                  tCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGAt                                    <  2:673917/139‑1 (MQ=255)
                                                                                   cGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATg                                   <  2:611073/139‑1 (MQ=255)
                                                                                      ttAGCTTCATGATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGAcc                                >  1:390556/1‑139 (MQ=255)
                                                                                                gATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGc                                                                                      >  1:707358/1‑75 (MQ=255)
                                                                                                gATTACCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGc                                                                                      <  2:707358/75‑1 (MQ=255)
                                                                                                    aCCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGAt                                    <  1:199055/120‑1 (MQ=255)
                                                                                                    aCCTCTGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGAt                                    >  2:199055/1‑120 (MQ=255)
                                                                                                         tGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGttt                                                                                                       <  1:254108/49‑1 (MQ=255)
                                                                                                         tGCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGttt                                                                                                       >  2:254108/1‑49 (MQ=255)
                                                                                                          gCATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGa            <  2:61632/139‑1 (MQ=255)
                                                                                                           cATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACaga                                                                            <  1:24875/74‑1 (MQ=255)
                                                                                                           cATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACaga                                                                            >  2:24875/1‑74 (MQ=255)
                                                                                                           cATCCTGGCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGAt           >  2:472259/1‑139 (MQ=255)
                                                                                                                 ggCGCATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAg     <  2:636110/139‑1 (MQ=255)
                                                                                                                  gcgcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCt                                                >  2:425346/1‑95 (MQ=255)
                                                                                                                  gcgcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCt                                                <  1:425346/95‑1 (MQ=255)
                                                                                                                    gcATTTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTga  <  1:601683/139‑1 (MQ=255)
                                                                                                                        |                                                                                                                                      
ATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATACAGAGCTTCATCGGCTATTTGAATGAGGCGCTCATAGTCAGGATGACCATTAAACATGGCGGCACCGATGGAAAGTGA  >  NC_000913/1565638‑1565892

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