Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,881,483 Δ1 bp intergenic (‑126/‑482) dnaA ← / → rpmH DnaA/50S ribosomal subunit protein L34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,881,4830G.97.9% 621.8 / 3.4 142intergenic (‑126/‑482)dnaA/rpmHDnaA/50S ribosomal subunit protein L34
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base . (64/75);  total (67/75)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 142 total aligned reads displayed.

TGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAGTCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCA  >  NC_000913/3881413‑3881553
                                                                      |                                                                      
tgttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGGACCGCTGTCTGCGGTTATATGCCCGATCAAGa                                                                  <  2:4300373/76‑1 (MQ=255)
tgttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGa                                                                  <  2:1409902/76‑1 (MQ=255)
 gttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGAt                                                                 <  2:5470996/76‑1 (MQ=255)
  ttgACTGACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                <  2:2202258/76‑1 (MQ=255)
  ttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                >  1:3065044/1‑76 (MQ=255)
  ttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                >  2:3410270/1‑76 (MQ=255)
  ttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                <  1:828494/76‑1 (MQ=255)
  ttgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                <  2:2978323/76‑1 (MQ=255)
  gtgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATc                                                                <  2:1642109/75‑1 (MQ=255)
   tgACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATcc                                                               <  2:2685213/76‑1 (MQ=255)
    gACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCt                                                              >  1:3745520/1‑76 (MQ=255)
       gtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCa                                                           >  1:3373966/1‑76 (MQ=255)
       gtacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCa                                                           >  2:2052002/1‑76 (MQ=255)
        tacgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaa                                                          >  1:172748/1‑76 (MQ=255)
          cgtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                        <  2:3673159/76‑1 (MQ=255)
          cggcGAAAAGACCCTGATTGAGGGTTACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCaaaa                                                        <  2:164355/73‑1 (MQ=255)
           gtCGAAAAGAGCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                       <  2:5663850/76‑1 (MQ=255)
           gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                       >  2:3512218/1‑76 (MQ=255)
           gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                       <  2:3038650/76‑1 (MQ=255)
           gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                       <  1:3959457/76‑1 (MQ=255)
           gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                       <  1:2197463/76‑1 (MQ=255)
           gtCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAAc                                                       <  2:5413448/76‑1 (MQ=255)
            tCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACg                                                      <  2:3061875/76‑1 (MQ=255)
             cGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGa                                                     <  1:4680268/76‑1 (MQ=255)
              gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAt                                                    <  2:4670466/76‑1 (MQ=255)
              gAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGAt                                                    <  2:2738612/76‑1 (MQ=255)
                aaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCg                                                  >  2:5364731/1‑76 (MQ=255)
                aaaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCg                                                  >  2:491022/1‑76 (MQ=255)
                 aaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCgg                                                 >  2:1842797/1‑76 (MQ=255)
                 aaGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCgg                                                 <  1:4916555/76‑1 (MQ=255)
                  aGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCggg                                                >  1:1920806/1‑76 (MQ=255)
                  aGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCggg                                                >  1:2333195/1‑76 (MQ=255)
                      ccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcg                                            >  1:872666/1‑76 (MQ=255)
                      ccTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcg                                            <  1:2182634/76‑1 (MQ=255)
                       cTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGCTATATGCCGGATCAAGATCCTGCAAAACGATCGGGCCcgc                                           <  2:2458110/76‑1 (MQ=255)
                        tGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcgcg                                          >  1:398116/1‑76 (MQ=255)
                        tGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACcgcg                                          >  2:2025166/1‑76 (MQ=255)
                         tcTTGAGTGTGACGCACGAACCACTGGCTGCGGTTATATGCCCGATCAATATCCGGCAAAACGATCGGGACCGCgg                                         <  1:4060042/75‑1 (MQ=255)
                          cTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                        <  1:95819/76‑1 (MQ=255)
                          cTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                        >  1:4000685/1‑76 (MQ=255)
                          cTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGa                                        <  1:2528902/76‑1 (MQ=255)
                            tGAGTGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATc                                      <  2:3680028/76‑1 (MQ=255)
                            tGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGCTCGGGACCGCGGATc                                      >  2:4979126/1‑76 (MQ=255)
                            tGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATc                                      <  1:5525259/76‑1 (MQ=255)
                             gggggTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCa                                     <  2:4331368/74‑1 (MQ=255)
                             gAGGGTGACGCACGAACCACTGTCGGCGGTTAGATGCCCGATCAAGATCCTGCAAAACGATCGGGCCCGCGGATCa                                     >  2:1776439/1‑76 (MQ=255)
                              aGGGTGACGCTCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                    <  1:1570166/76‑1 (MQ=255)
                              aGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCAt                                    >  1:1034541/1‑76 (MQ=255)
                               gggTGACGCACGAACCGCTGTCTGGGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                   <  2:3576010/76‑1 (MQ=255)
                               gggTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                   >  2:4030716/1‑76 (MQ=255)
                               gggTGACGCACGAACCACTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATa                                   <  1:986096/76‑1 (MQ=255)
                                ggTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAg                                  <  2:4037508/76‑1 (MQ=255)
                                 gTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGc                                 <  1:5227493/76‑1 (MQ=255)
                                  tGACGCCCGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGcc                                <  1:665536/76‑1 (MQ=255)
                                   gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCt                               <  1:774460/76‑1 (MQ=255)
                                   gACGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCt                               <  1:1926323/76‑1 (MQ=255)
                                    gcgcgcgAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTa                              <  2:1986008/75‑1 (MQ=255)
                                    gcGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTa                              <  2:4244934/75‑1 (MQ=255)
                                    aCGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTa                              >  2:4758387/1‑76 (MQ=255)
                                     cGCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaa                             >  1:2720061/1‑76 (MQ=255)
                                      gCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaaa                            >  1:1254577/1‑76 (MQ=255)
                                      gCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaaa                            >  2:2434389/1‑76 (MQ=255)
                                      gCACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaaa                            >  1:3499623/1‑76 (MQ=255)
                                       cACGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                           >  2:1277391/1‑76 (MQ=255)
                                         cGAACCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTg                         >  1:2546907/1‑76 (MQ=255)
                                            aCCGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTgcgc                      <  2:751162/76‑1 (MQ=255)
                                              cGCTGTCTGCGGTTATATGCCCTATCAAGAGCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCaa                    <  1:3465669/76‑1 (MQ=255)
                                              cGCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCaa                    <  1:4503978/76‑1 (MQ=255)
                                               gCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAag                   >  2:2655577/1‑76 (MQ=255)
                                               gCTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAag                   >  1:4425944/1‑76 (MQ=255)
                                                cTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAaga                  <  2:4952758/76‑1 (MQ=255)
                                                cTGTCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAaga                  <  2:4386876/76‑1 (MQ=255)
                                                   tCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAt               >  2:2648814/1‑76 (MQ=255)
                                                   tCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAt               <  1:2868356/76‑1 (MQ=255)
                                                   tCTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAt               >  2:4141642/1‑76 (MQ=255)
                                                   tCTGCGGTTATATGAACGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAt               >  1:1871472/1‑76 (MQ=255)
                                                    cTGCGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtc              <  2:2629968/76‑1 (MQ=255)
                                                      gCGGTGATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtctt            <  2:5228829/76‑1 (MQ=255)
                                                       cGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaaa                            <  1:1334624/59‑1 (MQ=255)
                                                       cGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTaaa                            >  2:1334624/1‑59 (MQ=255)
                                                       cGGTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttc           >  2:181364/1‑76 (MQ=255)
                                                        ggTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct          >  2:2309708/1‑76 (MQ=255)
                                                        ggTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGAtcttct          <  2:5164149/76‑1 (MQ=255)
                                                         gTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAACCTGCGCAAGAGATCTTCTg         <  1:3077048/76‑1 (MQ=255)
                                                         gTTATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTg         >  2:4967625/1‑76 (MQ=255)
                                                          ttATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt        >  1:929324/1‑76 (MQ=255)
                                                          ttATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt        <  1:2814445/76‑1 (MQ=255)
                                                          ttATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt        >  1:2315313/1‑76 (MQ=255)
                                                          ttATATGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGt        <  1:4606370/76‑1 (MQ=255)
                                                            atatGCCCGAGCAAGATCCTGCAAAACGATTGGGACCGCGGATAATAGCCTAAACTGCGCAAGAGATCTTGTGttt      >  1:4717342/1‑76 (MQ=255)
                                                            atatGCCCGATCAAGATCCGGCAAAACGATCGGGACCGCGGATCATAGCCTAAAc                           <  2:3499256/55‑1 (MQ=255)
                                                             tatGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtc     >  2:4352178/1‑76 (MQ=255)
                                                             tatGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtc     <  2:1496166/76‑1 (MQ=255)
                                                              atGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct    >  1:3040358/1‑76 (MQ=255)
                                                              atGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct    >  2:663513/1‑76 (MQ=255)
                                                              atGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct    <  1:3973707/76‑1 (MQ=255)
                                                              atGCCCGATCAAGAGCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtct    <  2:5636723/76‑1 (MQ=255)
                                                               tGCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTtctc   >  1:5540422/1‑76 (MQ=255)
                                                                gCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTca  >  2:3277179/1‑76 (MQ=255)
                                                                gCCCGATCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTca  <  1:445473/76‑1 (MQ=255)
                                                                      |                                                                      
TGTTGACGTACGTCGAAAAGACCCTGCTTGAGGGTGACGCACGAACCGCTGTCTGCGGTTATATGCCCGAGTCAAGATCCTGCAAAACGATCGGGACCGCGGATCATAGCCTAAACTGCGCAAGAGATCTTCTGTTTCTCA  >  NC_000913/3881413‑3881553

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