Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,264,347 Δ1 bp 11.1% intergenic (‑53/+31) 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,3470A.11.1% 141.2 / 14.8 54intergenic (‑53/+31)zinU/PP_4640P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base . (4/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

CTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACATAG  >  NC_002947/5264202‑5264495
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                          <  2:700613/64‑1 (MQ=18)
                                                                                                                              cTGTGGGAGCCTGGGTAACCGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                          >  1:700613/1‑64 (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)
                                                                                                                                          cgatctccAGGCTCCCACAGGTTTACACGACCCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTc         <  2:389162/144‑1 (MQ=255)
                                                                                                                                           gTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                               <  1:1323828/95‑1 (MQ=255)
                                                                                                                                           gTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                               >  2:1323828/1‑95 (MQ=255)
                                                                                                                                              aCCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACa     <  2:472887/150‑1 (MQ=255)
                                                                                                                                              aCCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACa     >  1:872026/1‑150 (MQ=255)
                                                                                                                                                 cAGGCTCCCACAGGTTTACGCGAACCGCAGGCTCAAAAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTGCTGACGTTCGCGGAAGAATGCCTCGTAGCTCAGCACCGGCCTGTCCGGATGCGTTTTCTACAGGTGTTCGAGAtcc  >  1:389162/1‑148 (MQ=255)
                                                                                                                                                 cAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTcc                                                                             <  2:552838/75‑1 (MQ=255)
                                                                                                                                                 cAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTcc                                                                             >  1:552838/1‑75 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCGGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACATAG  >  NC_002947/5264202‑5264495

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: