Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,264,341 T→G 13.6% intergenic (‑47/+37) zinU ← / ← PP_4640 P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,264,3410TG13.6% 114.0 / 11.7 50intergenic (‑47/+37)zinU/PP_4640P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/22);  new base G (5/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

CGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGT  >  NC_002947/5264195‑5264486
                                                                                                                                                   |                                                                                                                                                 
cGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTAc                                                                                                                                                <  1:491205/150‑1 (MQ=255)
  gCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTAccc                                                                                                                                              >  2:1119417/1‑150 (MQ=255)
       cTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCt                                                                                                                                         <  2:663066/150‑1 (MQ=255)
        ttCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTc                                                                                                                                        >  2:439461/1‑150 (MQ=255)
              ccATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                  <  1:765167/150‑1 (MQ=255)
              ccATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGAGGCCACAg                                                                                                                                  >  1:60290/1‑150 (MQ=255)
               cATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCGTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAgg                                                                                                                                 >  2:1174534/1‑150 (MQ=255)
                            gAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCg                                                                                                                    <  1:439461/150‑1 (MQ=255)
                               cATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAg                                                                                                                 <  2:2863/150‑1 (MQ=255)
                                aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAAGACTCCTGAATCTGAAAATGACTCTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAgg                                                                                                                <  2:640365/150‑1 (MQ=255)
                                   tACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTc                                                                                                             >  2:1267096/1‑150 (MQ=255)
                                           cAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCAc                                                                                                     >  1:140796/1‑150 (MQ=255)
                                              ggtggtCTTGCCGGCACCGAGGAAGCCTGTAAATCTACACACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTg                                                                                                  <  2:1181912/150‑1 (MQ=255)
                                                   tCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTGCTGAATCTGCAAGTGACTGTGGGAGCCGGTTACCCAGGCCCCCACAGGTTTACACGAACCGCAGGCCGAGCCGCACCTGCGccc                                                                                             >  2:48433/1‑150 (MQ=255)
                                                     ttGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCgg                                                                                           >  2:618881/1‑150 (MQ=255)
                                                          ggCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTg                                                                                      >  2:796842/1‑150 (MQ=255)
                                                             aCCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCa                                                                                   <  1:1267096/150‑1 (MQ=255)
                                                                     aGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACg                                                                           <  2:1330853/150‑1 (MQ=255)
                                                                      gCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCCACGAGGTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGt                                                                          >  1:375038/1‑150 (MQ=255)
                                                                            gAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTcc                                                                    <  2:688630/150‑1 (MQ=255)
                                                                               gACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCcaca                                                                                                                                   >  2:964784/1‑84 (MQ=25)
                                                                               gACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACCGGCTCCcaca                                                                                                                                   <  1:964784/84‑1 (MQ=25)
                                                                                    ttACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTc                                                            <  1:796842/150‑1 (MQ=255)
                                                                                    ttACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGAACTCGCGTCCGTAACGTGCTTCCTGACGttt                                                            >  1:714704/1‑149 (MQ=255)
                                                                                       cGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGttagga                                                         <  1:1119417/150‑5 (MQ=255)
                                                                                        gggAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCgg                                                        >  2:1335337/1‑150 (MQ=255)
                                                                                                gCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGTGTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTTCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAAttc                                                >  2:200530/1‑148 (MQ=255)
                                                                                                       cacagaaaaTACGCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                         <  1:48433/145‑1 (MQ=255)
                                                                                                       cACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                         <  1:823991/150‑1 (MQ=255)
                                                                                                       cACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                         <  1:1174534/150‑1 (MQ=255)
                                                                                                         cGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAg                                                                                                        <  2:610177/85‑1 (MQ=255)
                                                                                                         cGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAg                                                                                                        >  1:610177/1‑85 (MQ=255)
                                                                                                            aaaaTACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                  >  1:1380525/1‑56 (MQ=255)
                                                                                                            aaaaTACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAg                                                                                                                                  <  2:1380525/56‑1 (MQ=255)
                                                                                                               aTACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCAc                                 <  2:60290/150‑1 (MQ=255)
                                                                                                                tactacGGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCAcc                                >  1:1058984/1‑150 (MQ=255)
                                                                                                                 gcTCCTGAATCGGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCg                               <  2:454302/149‑1 (MQ=255)
                                                                                                                       gAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGt                         >  2:869762/1‑150 (MQ=255)
                                                                                                                          tCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCCCATCACCGCCTTGTCCg                      >  2:227101/1‑150 (MQ=255)
                                                                                                                             gAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGAt                   <  1:618881/150‑1 (MQ=255)
                                                                                                                              tgtaaaccTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATg                  >  2:842008/8‑150 (MQ=38)
                                                                                                                                     cTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                 >  1:700613/1‑64 (MQ=18)
                                                                                                                                     cTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                 <  2:700613/64‑1 (MQ=18)
                                                                                                                                     cTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTc           >  1:821490/1‑150 (MQ=255)
                                                                                                                                      tgtgGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt          >  1:96869/1‑150 (MQ=255)
                                                                                                                                          ggAGCCGGTTACCCAGGCTCCCACAGGCATGTACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt      >  1:37787/1‑150 (MQ=255)
                                                                                                                                            aGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAtgt    <  2:140796/150‑1 (MQ=255)
                                                                                                                                              ccGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATgtgt  >  1:472887/1‑150 (MQ=255)
                                                                                                                                                  gTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                      >  2:1323828/1‑95 (MQ=255)
                                                                                                                                                  gTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                      <  1:1323828/95‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGT  >  NC_002947/5264195‑5264486

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