Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,303,332 T→G V145G (GTG→GGG)  oppB → oligopeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,303,3320TG94.0% 170.0 / 2.7 50V145G (GTG→GGG) oppBoligopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (21/26);  total (22/28)
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.13e-01

TAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGTGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAATTCATGATATTG  >  NC_000913/1303237‑1303426
                                                                                               |                                                                                              
tAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGtt                                                                                            <  2:1154894/100‑1 (MQ=255)
tAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGtt                                                                                            <  1:736127/100‑1 (MQ=255)
 aGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTg                                                                                           <  1:622317/100‑1 (MQ=255)
   tGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCa                                                                                         <  1:1062093/100‑1 (MQ=255)
     cTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCAcc                                                                                       <  1:702638/100‑1 (MQ=255)
      tGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGGTTTGGGGTTGCAcc                                                                                       >  2:16098/1‑99 (MQ=255)
        gCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGTGGTTGCACCatt                                                                                    <  2:240682/100‑1 (MQ=255)
         cGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCatta                                                                                   >  2:826980/1‑100 (MQ=255)
             aTTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAg                                                                               >  2:442687/1‑100 (MQ=255)
               ttccGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCatta                                                                                   <  2:1303737/92‑1 (MQ=255)
                 ccGAATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCAt                                                                           <  1:114700/100‑1 (MQ=255)
                  cGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATg                                                                          <  2:189733/100‑1 (MQ=255)
                     aTTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGata                                                                       >  1:612138/1‑100 (MQ=255)
                     aTTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGata                                                                       <  1:1225070/100‑1 (MQ=255)
                     aTTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGata                                                                       <  1:1226551/100‑1 (MQ=255)
                      ttAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGatat                                                                      >  2:204748/1‑100 (MQ=255)
                         aaacaaaacaCCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTg                                                                   >  1:200593/1‑100 (MQ=255)
                         aaacaaaacaCCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTg                                                                   >  1:461296/1‑100 (MQ=255)
                           acaaaacaCCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCg                                                                 <  1:929312/100‑1 (MQ=255)
                            caaaacaCCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGTGGTTGCACCATTATTAGTCATGATATTTGCGa                                                                <  1:747858/100‑1 (MQ=255)
                             aaaacaCCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGAt                                                               <  2:1303260/100‑1 (MQ=255)
                                  acCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCAttt                                                          <  1:174847/100‑1 (MQ=255)
                                   ccAAATGGGATTATACCGTGATGGTGCTGGCAATGACCGGGGTCGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCAtttt                                                         <  2:831929/100‑1 (MQ=255)
                                    cAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTg                                                        <  2:401939/100‑1 (MQ=255)
                                      aaTGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCa                                                      <  2:890718/100‑1 (MQ=255)
                                      aaTGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCa                                                      <  1:628427/100‑1 (MQ=255)
                                        tGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCAtt                                                    >  1:983007/1‑100 (MQ=255)
                                          ggACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTgg                                                  >  2:248684/1‑100 (MQ=255)
                                          ggACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTgg                                                  <  1:495614/100‑1 (MQ=255)
                                               ataCCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGcc                                             <  2:1110013/100‑1 (MQ=255)
                                                taCCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCg                                            >  2:337251/1‑100 (MQ=255)
                                                  ccGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCggg                                          >  2:704544/1‑100 (MQ=255)
                                                      gATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGt                                      >  2:15989/1‑100 (MQ=255)
                                                           ggCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTg                                 >  1:338271/1‑100 (MQ=255)
                                                            gCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTgg                                <  1:339768/100‑1 (MQ=255)
                                                             cTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGa                               <  1:1018225/100‑1 (MQ=255)
                                                              tGGCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGaa                              <  1:946974/100‑1 (MQ=255)
                                                               ggCAATGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCCATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAAt                             >  1:1242424/1‑100 (MQ=255)
                                                                gCAATGACCGGGGTTGTTATCCCCAGTTTTGTGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAAtg                            >  1:104015/1‑100 (MQ=255)
                                                                  aaTGACCGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAAtggt                          >  2:338367/1‑100 (MQ=255)
                                                                       ccGGGGTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGgc                     >  2:941209/1‑100 (MQ=255)
                                                                          gggTTGTTATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCt                  <  1:338370/100‑1 (MQ=255)
                                                                              tgttATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTaaa              >  2:754284/1‑100 (MQ=255)
                                                                              tgttATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTaaa              <  2:842053/100‑1 (MQ=255)
                                                                               gttATCCCCAGTTTTGGGGTTGCACCCTTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAAt             >  1:793908/1‑100 (MQ=255)
                                                                               gttATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAAt             >  1:1280826/1‑100 (MQ=255)
                                                                                 tATCCCCAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAATTc           >  2:6568/1‑100 (MQ=255)
                                                                                    ccccAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAATTCATg        >  2:1114760/1‑100 (MQ=255)
                                                                                     cccAGTTTTGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAATTCATGa       <  1:630041/100‑1 (MQ=255)
                                                                                          ttttGGGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAATTCATGATATTg  <  1:6568/100‑1 (MQ=255)
                                                                                               |                                                                                              
TAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTTATCCCCAGTTTTGTGGTTGCACCATTATTAGTCATGATATTTGCGATCATTTTGCATTGGCTGCCGGGCGGTGGCTGGAATGGTGGGGCGCTTAAATTCATGATATTG  >  NC_000913/1303237‑1303426

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