Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,264,344 C→A 12.0% intergenic (‑50/+34) 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,3440CA12.0% 136.1 / 7.6 50intergenic (‑50/+34)zinU/PP_4640P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base A (4/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACA  >  NC_002947/5264195‑5264492
                                                                                                                                                      |                                                                                                                                                    
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)
                                                                                                                                                     aCCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACa  <  2:472887/150‑1 (MQ=255)
                                                                                                                                                     aCCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACa  >  1:872026/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACA  >  NC_002947/5264195‑5264492

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: