Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,225,612 A→T L7Q (CTG→CAG)  yohC ← Yip1 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,225,6120AT98.5% 228.4 / ‑3.6 65L7Q (CTG→CAG) yohCYip1 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (25/39);  total (25/40)
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.97e-01

TCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCC  >  NC_000913/2225517‑2225705
                                                                                               |                                                                                             
tCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTcc                                                                                           <  2:118069/100‑1 (MQ=255)
  gccgccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTcccc                                                                                         >  1:686185/1‑100 (MQ=255)
  gccgccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTcccc                                                                                         >  1:760341/1‑100 (MQ=255)
         tCAGCAAAACGTTGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACTTGCATTTCACGATCGGGATGGTAACACTGAACCCAg                                                                                       >  2:808702/1‑95 (MQ=255)
         tCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATg                                                                                  >  1:894862/1‑100 (MQ=255)
           aGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGc                                                                                <  2:1316343/100‑1 (MQ=255)
           aGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGc                                                                                <  1:39088/100‑1 (MQ=255)
            gCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCt                                                                               <  2:54915/100‑1 (MQ=255)
            gCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCt                                                                               >  2:617344/1‑100 (MQ=255)
             cAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTc                                                                              <  2:760337/100‑1 (MQ=255)
              aaaaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCa                                                                             <  1:254976/100‑1 (MQ=255)
                 aCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATAc                                                                          >  2:1049258/1‑100 (MQ=255)
                       tgggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGcccc                                                                    >  1:712061/1‑100 (MQ=255)
                        gggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGccccc                                                                   >  1:995800/1‑100 (MQ=255)
                            gtAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCttgt                                                               >  1:1280698/1‑100 (MQ=255)
                              aaTGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCttgttg                                                             <  2:256039/100‑1 (MQ=255)
                              aaTGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCttgttg                                                             <  2:977230/100‑1 (MQ=255)
                               aTGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTtgttgt                                                            <  2:933496/100‑1 (MQ=255)
                                tGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTtgtg                                                           <  2:1442913/100‑1 (MQ=255)
                                tGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTtgtg                                                           <  1:327754/100‑1 (MQ=255)
                                      gaAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCa                                                     <  2:645352/100‑1 (MQ=255)
                                        aaTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATg                                                   <  2:122199/100‑1 (MQ=255)
                                        aaTCGTTTCGTTTTCACGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATg                                                   <  1:522544/100‑1 (MQ=255)
                                            gtttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAAcc                                               >  1:1097328/1‑100 (MQ=255)
                                                cgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGAc                                           <  1:534869/100‑1 (MQ=255)
                                                cgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGAc                                           <  1:854470/100‑1 (MQ=255)
                                                 gtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACa                                          <  1:250338/100‑1 (MQ=255)
                                                  ttttCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAg                                         <  2:925285/100‑1 (MQ=255)
                                                   tttCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGt                                        >  1:416970/1‑100 (MQ=255)
                                                    ttCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCGTGTTGGGAGTTCATGAACCTGACAGt                                        >  1:231562/1‑99 (MQ=255)
                                                    tccgcgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAg                                         >  2:1071890/3‑98 (MQ=255)
                                                     tCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTc                                      <  1:351136/100‑1 (MQ=255)
                                                     tCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTc                                      >  1:1570739/1‑100 (MQ=255)
                                                      cgcgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCa                                     >  2:1492482/1‑100 (MQ=255)
                                                        cgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTTTGCGTTCATGAACCCGACAGTTCTAg                                   >  2:467824/1‑100 (MQ=255)
                                                        cgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                                   <  1:617347/100‑1 (MQ=255)
                                                        cgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                                   <  2:553052/100‑1 (MQ=255)
                                                        cgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                                   <  2:1370112/100‑1 (MQ=255)
                                                        cgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                                   <  2:1370100/100‑1 (MQ=255)
                                                         gATTAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCTCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                                   <  1:183926/99‑1 (MQ=255)
                                                           ttAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtat                                <  2:1394643/100‑1 (MQ=255)
                                                           ttAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtat                                <  2:1230858/100‑1 (MQ=255)
                                                            tAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtata                               <  1:663660/100‑1 (MQ=255)
                                                            tAATCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtata                               <  1:615143/100‑1 (MQ=255)
                                                              aTCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAg                             <  2:959631/100‑1 (MQ=255)
                                                              aTCACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAg                             >  1:851369/1‑100 (MQ=255)
                                                                cACCTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTc                           <  1:729330/100‑1 (MQ=255)
                                                                  ccTGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAg                         >  2:367436/1‑100 (MQ=255)
                                                                    tGCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCt                       <  1:1071897/100‑1 (MQ=255)
                                                                     gCATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCtt                      >  2:1190491/1‑100 (MQ=255)
                                                                      cATTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtg                     <  2:1365665/100‑1 (MQ=255)
                                                                       aTTTCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgt                    >  1:1520718/1‑100 (MQ=255)
                                                                         ttCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGa                  >  1:401055/1‑100 (MQ=255)
                                                                          tCACGATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGat                 <  1:1190498/100‑1 (MQ=255)
                                                                              gATCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattatt             >  1:1151241/1‑100 (MQ=255)
                                                                               aTCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAttt            <  2:378324/100‑1 (MQ=255)
                                                                                tCGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAtttt           >  2:1310850/1‑100 (MQ=255)
                                                                                 cGGGATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAttttt          <  2:1044994/100‑1 (MQ=255)
                                                                                  gggATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAtttttt         >  1:1283854/1‑100 (MQ=255)
                                                                                  gggATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGCATACGTCCGCTTGTGCTTAtttt           >  2:1481280/1‑98 (MQ=255)
                                                                                   ggATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTg        <  1:1419165/100‑1 (MQ=255)
                                                                                   ggATGGGAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTg        <  1:1478433/100‑1 (MQ=255)
                                                                                   ggATGGAAAAACTGTCTCCAGACATGGCTCAGACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAtttttt         <  2:619768/99‑1 (MQ=255)
                                                                                         gAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGcc  <  2:1087162/100‑1 (MQ=255)
                                                                                         gAAAACTGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGcc  <  1:474729/100‑1 (MQ=255)
                                                                                               |                                                                                             
TCGCCGCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTTCGCC  >  NC_000913/2225517‑2225705

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