Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,340,622:1 +A 10.0% coding (249/1587 nt) ybiT ← ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,340,6221.A10.0% 173.0 / 15.7 50coding (249/1587 nt)ybiTABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base . (22/23);  new base A (3/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

TCGGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGC  >  NC_002947/3340476‑3340768
                                                                                                                                                   |                                                                                                                                                  
tCGGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGCTACGGTAACGCTAAGCTTTGACTTACCACGGCTGTCCGTGGTTCATGTTCACCGTGTCg                                                                                                                                                 >  2:615270/1‑150 (MQ=255)
    cggTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCa                                                                                                                                             <  1:932389/150‑1 (MQ=255)
        gtAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGt                                                                                                                                         <  2:1070907/150‑1 (MQ=255)
         tAGCCGTCCATTTCGGTGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTg                                                                                                                                        <  2:279123/150‑1 (MQ=255)
              gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACTGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGatc                                                                                                                                     >  2:581670/1‑146 (MQ=255)
              gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACTGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGatc                                                                                                                                     <  1:581670/148‑3 (MQ=255)
              gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTc                                                                                                                                   >  2:409070/1‑150 (MQ=255)
               tCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGatcatg                                                                                                                                  <  1:162867/150‑6 (MQ=255)
                ccATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCt                                                                                                                                 <  1:20491/150‑1 (MQ=255)
                     tCGGCGAATTCTTTTTCCCCCTCGGCGCCGGCCATGCCGTCTTCTTCGGTCATTCCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCGGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAg                                                                                                                            <  1:535679/150‑1 (MQ=255)
                       ggCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGt                                                                                                                                         >  1:989075/1‑135 (MQ=255)
                       ggCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGt                                                                                                                                         <  2:989075/135‑1 (MQ=255)
                                ggTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTc                                                                                                                 <  1:1170291/150‑1 (MQ=255)
                                     ccAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGgc                                                                                                            <  1:668085/150‑1 (MQ=255)
                                         cTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGa                                                                                                                                                <  1:1125932/110‑1 (MQ=255)
                                         cTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGa                                                                                                                                                >  2:1125932/1‑110 (MQ=255)
                                            ggCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCtt                                                                                                     >  2:883876/1‑150 (MQ=255)
                                                    ccATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGc                                                                                             >  1:1158133/1‑150 (MQ=255)
                                                       tGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCggg                                                                                          >  1:687174/1‑150 (MQ=255)
                                                                         tttCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGc                                                                                                                          >  1:476134/1‑100 (MQ=255)
                                                                         tttCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGc                                                                                                                          <  2:476134/100‑1 (MQ=255)
                                                                             cGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTg                                                                    >  1:950769/1‑150 (MQ=255)
                                                                               gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                  <  1:524462/150‑1 (MQ=255)
                                                                                  gCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACt                                                                                                                     <  2:899083/96‑1 (MQ=255)
                                                                                  gCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACt                                                                                                                     >  1:899083/1‑96 (MQ=255)
                                                                                     aGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCgg                                                            >  1:880891/1‑150 (MQ=255)
                                                                                          aTACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGc                                                       >  2:404197/1‑150 (MQ=255)
                                                                                               gTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAg                                                  <  1:409070/150‑1 (MQ=255)
                                                                                                 cACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGt                                                <  2:687174/150‑1 (MQ=255)
                                                                                                          cTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGa                                       <  2:1095447/150‑1 (MQ=255)
                                                                                                          cTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCCTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGa                                       <  1:615270/150‑1 (MQ=255)
                                                                                                           tttGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAg                                      <  2:1158133/150‑1 (MQ=255)
                                                                                                           tttGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAg                                      >  2:414500/1‑150 (MQ=255)
                                                                                                               aCTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                              >  2:1202534/1‑106 (MQ=255)
                                                                                                               aCTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                              <  1:1202534/106‑1 (MQ=255)
                                                                                                               aCTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGAtt                                  >  2:376147/1‑150 (MQ=255)
                                                                                                                 tttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGAtttt                                >  2:387838/1‑150 (MQ=255)
                                                                                                                  ttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTc                               >  1:547208/1‑150 (MQ=255)
                                                                                                                        aGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAg                         >  2:1055961/1‑150 (MQ=255)
                                                                                                                          cTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTTACCGCCGAGGATTTTCATGAAGGt                       <  2:950769/150‑1 (MQ=255)
                                                                                                                                gTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACtt                 <  1:376147/150‑1 (MQ=255)
                                                                                                                                   gCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTACGCTAACTGGTCCTGACGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCCGACGGCTCCAGGTCACCCCCGAGGATTTTCCTGAACGTCGACTTGcc              >  1:151954/1‑150 (MQ=255)
                                                                                                                                      ggaattCACCGTGATCGACACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCGGCa           >  1:123050/6‑150 (MQ=255)
                                                                                                                                      cATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAAc                                                                                                                      >  2:977197/1‑43 (MQ=255)
                                                                                                                                      cATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAAc                                                                                                                      <  1:977197/43‑1 (MQ=255)
                                                                                                                                        aattCACCGTGATCGACACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAAc         >  2:770135/4‑150 (MQ=255)
                                                                                                                                           tCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCgg                                                                                   <  2:1030533/73‑1 (MQ=255)
                                                                                                                                           tCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCgg                                                                                   >  1:1030533/1‑73 (MQ=255)
                                                                                                                                             aCCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTg    <  2:880891/150‑1 (MQ=255)
                                                                                                                                               cGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGgc  >  1:476716/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGC  >  NC_002947/3340476‑3340768

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: