Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,381,721 G→C 11.3% D156H (GAC→CAC)  dinB → DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,381,7210GC11.3% 152.0 / 11.4 53D156H (GAC→CAC) dinBDNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base C (2/4);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

TCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGAAAC  >  NC_002947/1381579‑1381857
                                                                                                                                              |                                                                                                                                        
tCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGgcg                                                                                                                                   <  2:46063/150‑1 (MQ=255)
  tGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTCACTGgcgct                                                                                                                                 <  2:382011/150‑2 (MQ=255)
  tGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTCACTGgcgat                                                                                                                                 <  2:593190/150‑3 (MQ=255)
        tCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCa                                                                                                                           >  1:879963/1‑150 (MQ=255)
           gCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATg                                                                                                                        >  2:1001889/1‑150 (MQ=255)
                    aGTGGTATTCGGGCAGCGCCACTCGCATTTCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCg                                                                                                               <  2:1150526/150‑1 (MQ=255)
                    aGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGAGATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCg                                                                                                               >  1:1150526/1‑150 (MQ=255)
                          aTTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGNCAAGATNGCCAGTGACTGGNGCAAACCCAATGGGCTGTTCGTGATTa                                                                                                         <  1:1058207/150‑1 (MQ=255)
                           ttCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATggg                                                                                                                      <  2:285166/136‑1 (MQ=255)
                           ttCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATggg                                                                                                                      >  1:285166/1‑136 (MQ=255)
                           ttCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTAc                                                                                                        <  1:289935/150‑1 (MQ=255)
                           ttCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTCACTGgcgatc                                                                                                                                <  2:72009/126‑4 (MQ=255)
                           ttCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTCACTGgcgatc                                                                                                                                >  1:72009/1‑123 (MQ=255)
                             cGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTAccc                                                                                                      <  2:554499/150‑1 (MQ=255)
                                          gcgcATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCaaa                                                                                                  >  1:1121330/1‑141 (MQ=255)
                                          gcgcATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCaaa                                                                                                  <  2:1121330/141‑1 (MQ=255)
                                          gcgcATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGa                                                                                         <  2:1332790/150‑1 (MQ=255)
                                             cATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAAc                                                                                      <  1:115082/150‑1 (MQ=255)
                                             cATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAAc                                                                                      <  1:1140421/150‑1 (MQ=255)
                                                    gAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAAcgttcgtt                                                                               >  1:667265/1‑150 (MQ=255)
                                                               ccGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGcc                                                                    <  2:879963/150‑1 (MQ=255)
                                                                cGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGTCTGTTCGTGATTACCCCAAATGAAGTGGAAAcgtt                                                                                   >  2:320775/1‑134 (MQ=255)
                                                                cGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGTCTGTTCGTGATTACCCCAAATGAAGTGGAAAcgtt                                                                                   <  1:320775/134‑1 (MQ=255)
                                                                  gCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGt                                                                 <  1:1112342/150‑1 (MQ=255)
                                                                  gCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGt                                                                 <  2:1146670/150‑1 (MQ=255)
                                                                   cGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACGCAATGGGCTGTTCCTGATTGCCTCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTg                                                                <  2:410173/150‑1 (MQ=255)
                                                                    gtgtTGCCCGCCCCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTgg                                                               <  1:1001889/150‑1 (MQ=255)
                                                                          cccGCACCCTGCATATCACAGTGTGGGCTGGCGTGACGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTGGTGATTATCCAAAATGCAGTGGAAACGTTCGTTGAAGCCTTGCCGGTGGCCAGGt                                                         >  1:1148406/1‑150 (MQ=255)
                                                                           ccGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAAc                                                                                      <  2:857502/120‑1 (MQ=255)
                                                                           ccGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAAc                                                                                      >  1:857502/1‑120 (MQ=255)
                                                                               aCCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCAt                                                    >  2:381885/1‑150 (MQ=255)
                                                                                   tGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATgggg                                                >  2:1204735/1‑150 (MQ=255)
                                                                                    gCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGAATGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATggggg                                               >  1:1098337/1‑149 (MQ=255)
                                                                                            ccGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAgg                                       <  1:1354049/150‑1 (MQ=255)
                                                                                            ccGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAgg                                       <  1:1204735/150‑1 (MQ=255)
                                                                                            ccGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATAGGGTGGGTAAgg                                       <  2:707179/150‑1 (MQ=255)
                                                                                                   ggCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGc                                >  1:1263977/1‑150 (MQ=255)
                                                                                                            ggcgCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACt                       <  1:689772/150‑1 (MQ=255)
                                                                                                             gcgcCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTg                      >  2:295515/1‑150 (MQ=255)
                                                                                                               gcCTAAAAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCACAGCGCTACAAAGAATTAACCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGAc                    >  2:1058783/1‑150 (MQ=255)
                                                                                                                     cAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAAc                                                                                                                              >  2:540391/1‑38 (MQ=255)
                                                                                                                     cAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAAc                                                                                                                              <  1:540391/38‑1 (MQ=255)
                                                                                                                      aaGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTg             >  1:6136/1‑150 (MQ=255)
                                                                                                                       aGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTgg            <  1:1139019/150‑1 (MQ=255)
                                                                                                                          ttCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCa         >  1:546619/1‑150 (MQ=255)
                                                                                                                           tCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCAt        <  2:1028095/150‑1 (MQ=255)
                                                                                                                             tGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAAcgttc                                                                                  >  1:1157103/1‑74 (MQ=255)
                                                                                                                             tGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAAcgttc                                                                                  <  2:1157103/74‑1 (MQ=255)
                                                                                                                               gCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGaa    >  2:1176258/1‑150 (MQ=255)
                                                                                                                               gCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTCCAGCCTTGCCGGTGGCCAGGTTGCATGGGGGAAGTAAGGTGACGGGGGGCAAACTGACGCGGCTGGGCATCGaa    >  1:804640/1‑150 (MQ=255)
                                                                                                                                 cAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGAAAc  >  2:393997/1‑150 (MQ=255)
                                                                                                                                  gtttgcGCCAGTCACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGAc                                   >  1:1233810/6‑116 (MQ=255)
                                                                                                                                  gtttgcGCCAGTCACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGAc                                   <  2:1233810/111‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                        
TCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCCAAATGAAGTGGAAACGTTCGTTGCAGCCTTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGAAAC  >  NC_002947/1381579‑1381857

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