Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,881,837 G→T T162N (ACC→AAC)  dnaN ← DNA polymerase III, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,881,8370GT95.5% 221.3 / ‑0.0 67T162N (ACC→AAC) dnaNDNA polymerase III, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (26/38);  total (27/40)
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

ATCACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGGTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTGCCTGCG  >  NC_000913/3881741‑3881930
                                                                                                |                                                                                             
aTCACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGtttt                                                                                            <  2:469707/100‑1 (MQ=255)
aTCACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGtttt                                                                                            <  1:654488/100‑1 (MQ=255)
aTCACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGtttt                                                                                            >  1:1112417/1‑100 (MQ=255)
 tCACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTc                                                                                           <  2:1544424/100‑1 (MQ=255)
 tCACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTc                                                                                           <  2:1383801/100‑1 (MQ=255)
  cACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGttttga                                                                                          >  2:735954/1‑98 (MQ=255)
  cACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCa                                                                                          >  1:1422900/1‑100 (MQ=255)
  cACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCa                                                                                          >  2:1419545/1‑100 (MQ=255)
  cACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCa                                                                                          >  2:862889/1‑100 (MQ=255)
   aCCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCaa                                                                                         <  2:31298/100‑1 (MQ=255)
    ccGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCTGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCan                                                                                         <  2:1197267/99‑2 (MQ=255)
     cGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAAc                                                                                       <  1:1176179/100‑1 (MQ=255)
        aTGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGc                                                                                    >  2:725202/1‑100 (MQ=255)
           gCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATa                                                                                 >  1:1503897/1‑100 (MQ=255)
           gCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATa                                                                                 >  2:782394/1‑100 (MQ=255)
           gCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATa                                                                                 >  1:690731/1‑100 (MQ=255)
              tGGCAAAGCTTGACCAAGTGGCATTGAACAGACCGCCAAACGGTGGCCGTCGGTTGCAACGGTGCGCAGTTCTTGACCTTCGTTTTCAAACAGCATACCa                                                                              >  2:1454078/1‑100 (MQ=255)
               ggCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCAt                                                                             >  1:840731/1‑100 (MQ=255)
               ggCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCAt                                                                             >  2:1349503/1‑100 (MQ=255)
                  aaaGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTa                                                                          <  2:753084/100‑1 (MQ=255)
                    aGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAg                                                                        <  1:728119/100‑1 (MQ=255)
                    aGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAg                                                                        <  1:246864/100‑1 (MQ=255)
                    aGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAg                                                                        <  2:789293/100‑1 (MQ=255)
                     gATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGt                                                                       <  1:20038/100‑1 (MQ=255)
                      aTTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTa                                                                      <  1:862896/100‑1 (MQ=255)
                       ttGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTaa                                                                     >  2:71814/1‑100 (MQ=255)
                       ttGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTaa                                                                     >  2:1413144/1‑100 (MQ=255)
                       ttGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTaa                                                                     <  2:582811/100‑1 (MQ=255)
                         gACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATa                                                                   >  2:1229404/1‑100 (MQ=255)
                         gACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATa                                                                   <  1:1262615/100‑1 (MQ=255)
                         ctccAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATa                                                                   <  1:1250347/98‑1 (MQ=255)
                          aCCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAg                                                                  <  2:496145/100‑1 (MQ=255)
                           ccAATTGGCATTGAACAGACCGCCCGACTGTGGCCGTCCGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAg                                                                  <  2:497074/99‑1 (MQ=255)
                           ccAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATaag                                                                 >  2:1250339/1‑98 (MQ=255)
                           ccAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGc                                                                 >  1:105704/1‑100 (MQ=255)
                            cAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCg                                                                <  2:1218343/100‑1 (MQ=255)
                            cAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCg                                                                <  2:1422890/100‑1 (MQ=255)
                              aTTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGaa                                                              <  1:716409/100‑1 (MQ=255)
                                 ggCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGGTTTCAAACAGCATACCATTTAAGTAATAGCGAACGt                                                           >  1:1261822/1‑100 (MQ=255)
                                    aTTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCt                                                        <  1:1014878/100‑1 (MQ=255)
                                    aTTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCt                                                        <  1:1411529/100‑1 (MQ=255)
                                     ttGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCt                                                        <  1:1438482/99‑1 (MQ=255)
                                      tGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGa                                                      >  1:1002433/1‑100 (MQ=255)
                                      tGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGa                                                      >  2:421954/1‑100 (MQ=255)
                                       gAACAGGCCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGAt                                                     <  1:317220/100‑1 (MQ=255)
                                       gAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGAt                                                     >  2:1051203/1‑100 (MQ=255)
                                         aCAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATgc                                                   >  2:1481499/1‑100 (MQ=255)
                                            gACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGcc                                                <  1:196070/100‑1 (MQ=255)
                                              ccGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCAt                                              >  2:1013870/1‑100 (MQ=255)
                                                         tGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCTTAGAAAACTggg                                   <  1:1030776/100‑1 (MQ=255)
                                                         tGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTggg                                   <  2:1384609/100‑1 (MQ=255)
                                                         tGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAACTAATAGCGAACGTCCTGATGCGCCATAGAAAACTggg                                   <  2:841096/100‑1 (MQ=255)
                                                          ggCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGt                                  <  2:558156/100‑1 (MQ=255)
                                                             cGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGc                               >  2:770490/1‑100 (MQ=255)
                                                               tCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCtt                             <  1:1103661/100‑1 (MQ=255)
                                                               tCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCtt                             <  1:707714/100‑1 (MQ=255)
                                                                cGGTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTc                            <  2:105704/100‑1 (MQ=255)
                                                                 ggTTGCCACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCa                           <  1:554994/100‑1 (MQ=255)
                                                                      ccACGGTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGAGCCATAGAAAACTGGGTCGCTTCAATc                        >  2:944011/1‑98 (MQ=255)
                                                                          ggTGCGCAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGc                  <  1:1127622/100‑1 (MQ=255)
                                                                             gggcAGTTCTTCACCTTCGGTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTc               <  1:680877/98‑1 (MQ=255)
                                                                             gcgcAGTTCTTCACCTTCGGTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTc               <  2:1261813/100‑1 (MQ=255)
                                                                               gcAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCAt             <  1:71814/100‑1 (MQ=255)
                                                                                cAGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATc            <  2:161918/100‑1 (MQ=255)
                                                                                 aGTTCTTCACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCg           <  1:942277/100‑1 (MQ=255)
                                                                                   ttcttcACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGtt         <  2:418814/100‑1 (MQ=255)
                                                                                     cttcACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTGc       >  1:287586/1‑100 (MQ=255)
                                                                                      ttcACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTg        >  2:1442304/1‑98 (MQ=255)
                                                                                        cACCTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTGCCTg    >  1:754893/1‑100 (MQ=255)
                                                                                          ccTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTGCCTGCg  >  1:463319/1‑100 (MQ=255)
                                                                                          ccTTCGTTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTGCCTGCg  <  2:809065/100‑1 (MQ=255)
                                                                                                |                                                                                             
ATCACCGAATGGCTTGGCAAAGATTGACCAATTGGCATTGAACAGACCGCCAGACGGTGGCCGTCGGTTGCCACGGTGCGCAGTTCTTCACCTTCGGTTTCAAACAGCATACCATTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTGCCTGCG  >  NC_000913/3881741‑3881930

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