Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,881,373 A→G 9.5% I146T (ATT→ACT)  dmlR ← transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,881,3730AG 9.5% 119.2 / 4.4 42I146T (ATT→ACT) dmlRtranscriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base A (18/20);  new base G (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

AATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCA  >  NC_000913/1881225‑1881520
                                                                                                                                                    |                                                                                                                                                   
aaTGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGaa                                                                                                                                                     <  3:303577/149‑1 (MQ=255)
    aCCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATc                                                                                                                                                 <  2:18220/149‑1 (MQ=255)
      cAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCga                                                                                                                                               <  4:43635/149‑1 (MQ=255)
          cACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATc                                                                                                                                           <  3:475456/149‑1 (MQ=255)
            cGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGa                                                                                                                                         <  4:422771/149‑1 (MQ=255)
                   gCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTc                                                                                                                                  >  4:258834/1‑149 (MQ=255)
                     aGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGAt                                                                                                                                            <  3:288166/137‑1 (MQ=255)
                     aGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGAt                                                                                                                                            >  4:288166/1‑137 (MQ=255)
                       aaaGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTg                                                                                                                              >  1:45449/1‑149 (MQ=255)
                       aaaGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTg                                                                                                                              >  2:584705/1‑149 (MQ=255)
                         aGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAc                                                                                                                            >  3:573218/1‑149 (MQ=255)
                          gATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGAcc                                                                                                                           <  2:198537/149‑1 (MQ=255)
                           aTCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCa                                                                                                                          >  3:200761/1‑149 (MQ=255)
                            tCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaa                                                                                                                         >  4:446123/1‑149 (MQ=255)
                                    aTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCa                                                                                                                 >  1:228442/1‑149 (MQ=255)
                                    aTGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAACACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCAAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCa                                                                                                                 <  4:360012/149‑1 (MQ=255)
                                              gcgcCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTc                                                                                                       <  1:497222/149‑1 (MQ=255)
                                                                              gcAGAGAATACGGTAATTGGTCGCCAGTTTGCGTGGGATTAAATTAGGGGCGAGGTCATCGCTAGTGCGAATATCGAGATCGACACCTTCATTGACCAAATGGAGCAGTCGGTCTGCCACATCAAAACGTAATTCTCGTGGTGGATACg                                                                       >  4:67386/1‑149 (MQ=255)
                                                                                   gaATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTc                                                                  <  2:157194/149‑1 (MQ=255)
                                                                                    aaTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaaa                                                                                                                        >  2:282342/1‑94 (MQ=255)
                                                                                    aaTACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCaaa                                                                                                                        <  1:282342/94‑1 (MQ=255)
                                                                                      tACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAGTATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGcc                                                               >  3:455376/1‑149 (MQ=255)
                                                                                      tACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAGTATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGcc                                                               >  2:409637/1‑149 (MQ=255)
                                                                                          gTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAccg                                                           >  4:347543/1‑149 (MQ=255)
                                                                                              ttGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCtag                                                       <  1:218985/149‑2 (MQ=255)
                                                                                                       aGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCgg                                              <  4:90675/149‑1 (MQ=255)
                                                                                                        gTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGa                                             <  3:22507/149‑1 (MQ=255)
                                                                                                           tGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAGTATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGcc                                          <  1:409637/149‑1 (MQ=255)
                                                                                                           tGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAGTATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGcc                                          <  4:455376/149‑1 (MQ=255)
                                                                                                                  gATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGcc                                          <  1:105633/142‑1 (MQ=255)
                                                                                                                  gATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGcc                                          >  2:105633/1‑142 (MQ=255)
                                                                                                                    ttAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgt                                 >  2:584550/1‑149 (MQ=255)
                                                                                                                     tAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtc                                <  4:171162/149‑1 (MQ=255)
                                                                                                                       aaTTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTgtcgt                              <  3:435973/149‑1 (MQ=255)
                                                                                                                         ttAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTcc                            >  2:62697/1‑149 (MQ=255)
                                                                                                                             gggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCaaaa                        >  3:13122/1‑149 (MQ=255)
                                                                                                                             gggCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCaaaa                        >  4:170852/1‑149 (MQ=255)
                                                                                                                                     tCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAgct                >  4:42742/1‑149 (MQ=255)
                                                                                                                                        tCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAgctgct             <  1:295028/149‑1 (MQ=255)
                                                                                                                                             aaTGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGataa        <  3:67386/149‑1 (MQ=255)
                                                                                                                                             aaTGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGataa        <  1:409008/149‑1 (MQ=255)
                                                                                                                                             aaTGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATCCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGataa        <  3:346865/149‑1 (MQ=255)
                                                                                                                                                  gAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACgc   >  1:228551/1‑149 (MQ=255)
                                                                                                                                                   aaTATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACgca  <  2:132222/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATGACCAGGCACGGCAATGCAGAAAGATCGGTTAAATGTTTCGGCGCGCCATGTTGCGCAATAAATTCCGGCGACGCGCAGAGAATACGGTAATTGGTCGCCAGTTTGCGGGCGATTAAATTAGGGGCGATGTCATCGCCAATGCGAATATCGAGATCGACACCTTCATTGACCAAATCGACCAGTCGGTCTTCCACATCAAAACGTAATTCCAGTTGCGGATACGCCTTCGCCAGCGCTAATAACGCCGGAGCCACCACCTGTCGTCCAAAACCAAAGCTGCTGATAATACGCA  >  NC_000913/1881225‑1881520

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