Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 520,993 T→G V193G (GTC→GGC)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913520,9930TG94.2% 171.3 / ‑3.2 52V193G (GTC→GGC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (24/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACG  >  NC_000913/520897‑521088
                                                                                                |                                                                                               
cGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAc                                                                                              >  1:1011941/1‑100 (MQ=255)
 ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                             >  1:1228804/1‑100 (MQ=255)
 ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                             <  1:573004/100‑1 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                        >  1:783619/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                        >  1:679680/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                        >  2:879341/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                        >  2:925746/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                        >  1:1406895/1‑100 (MQ=255)
        ggTTTTAACCCCGTCCAGAGGGCTCCGCGTCTGATGATGAAGCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGtt                                                                                      <  2:1517804/100‑1 (MQ=255)
           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGAGGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                   <  2:478283/100‑1 (MQ=255)
              aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                >  2:475688/1‑100 (MQ=255)
              aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                >  2:227873/1‑100 (MQ=255)
               aCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAtt                                                                               <  1:925751/100‑1 (MQ=255)
               aCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAtt                                                                               >  1:646900/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                              >  2:106185/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                              >  2:967570/1‑100 (MQ=255)
                 cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                             <  2:143923/100‑1 (MQ=255)
                  ccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcgg                                                                            >  1:1579653/1‑100 (MQ=255)
                      ccAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCggcggc                                                                        >  1:864304/1‑100 (MQ=255)
                         gATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAAc                                                                     <  1:789312/100‑1 (MQ=255)
                               tCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAAc                                                               <  2:1579643/100‑1 (MQ=255)
                                 cgcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCa                                                             >  2:210775/1‑100 (MQ=255)
                                    gTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCt                                                          <  1:566524/100‑1 (MQ=255)
                                        gatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACGAGGCGTTATAAATTCGGCGGCAACGAGCACTAGCACGa                                                       >  2:845019/1‑99 (MQ=255)
                                        gatgatgaATCTGGCGGATGGCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAACTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGAc                                                      >  2:1500131/1‑100 (MQ=255)
                                          tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                    <  1:234808/100‑1 (MQ=255)
                                            atgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCt                                                  <  2:1322238/100‑1 (MQ=255)
                                            atgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCt                                                  <  2:639962/100‑1 (MQ=255)
                                            atgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCt                                                  >  1:645821/1‑100 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                 <  2:16771/100‑1 (MQ=255)
                                                         aTGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgtt                                     <  1:794145/100‑1 (MQ=255)
                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                    <  2:1555529/100‑1 (MQ=255)
                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                    <  2:379140/100‑1 (MQ=255)
                                                           gTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgt                                   >  2:958390/1‑100 (MQ=255)
                                                            tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                  <  2:747538/100‑1 (MQ=255)
                                                               aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcc                               <  2:657088/100‑1 (MQ=255)
                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                            >  1:599878/1‑100 (MQ=255)
                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                            >  1:367160/1‑100 (MQ=255)
                                                                   aaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                            >  1:470987/1‑99 (MQ=255)
                                                                    aCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGc                          <  1:1084363/100‑1 (MQ=255)
                                                                      cGGAGCCGTGCAGCCGGGGAGTCGGGTCTCCTTGTGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAACTGGTTGTTACCTCAGCtt                        <  2:1308695/100‑1 (MQ=255)
                                                                      cGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                         <  2:695518/99‑1 (MQ=255)
                                                                          gCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                    <  1:434217/100‑1 (MQ=255)
                                                                           ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                   <  2:853059/100‑1 (MQ=255)
                                                                             gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCg                 >  1:1140285/1‑100 (MQ=255)
                                                                               gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa               <  1:475689/100‑1 (MQ=255)
                                                                               gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCCCAGCTTAAACCCgca               >  2:252715/1‑98 (MQ=255)
                                                                                  gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa            >  1:1332207/1‑100 (MQ=255)
                                                                                   ccGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa            >  2:781978/1‑99 (MQ=255)
                                                                                         aGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGt     <  1:340737/100‑1 (MQ=255)
                                                                                          gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa    <  1:227873/100‑1 (MQ=255)
                                                                                            cGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACg  <  2:470986/100‑1 (MQ=255)
                                                                                                |                                                                                               
CGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACG  >  NC_000913/520897‑521088

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: