Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,308,182 A→T L145Q (CTG→CAG)  yojI ← microcin J25 efflux ABC transporter permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,308,1820AT97.8% 156.8 / ‑3.6 46L145Q (CTG→CAG) yojImicrocin J25 efflux ABC transporter permease/ATPase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (27/18);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GTCGCCATATGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACAT  >  NC_000913/2308086‑2308275
                                                                                                |                                                                                             
gTCGCCATATGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCa                                                                                            >  2:812710/1‑100 (MQ=255)
   gCCATATGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGGCCGACTGCAt                                                                                           >  2:1063438/1‑98 (MQ=255)
   gCCATATGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATcc                                                                                         >  1:755013/1‑100 (MQ=255)
   gCCATATGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATcc                                                                                         <  2:1236774/100‑1 (MQ=255)
    ccATATGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCca                                                                                        >  1:964824/1‑100 (MQ=255)
         tGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGAt                                                                                   >  1:257504/1‑100 (MQ=255)
          gTTTGTACACGCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCAGCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGAt                                                                                   >  2:1275382/1‑99 (MQ=255)
            ttGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATAcg                                                                                <  2:305381/100‑1 (MQ=255)
                acacCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATAcgccgc                                                                            <  1:29275/100‑1 (MQ=255)
                acacCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATAcgccgc                                                                            <  1:977974/100‑1 (MQ=255)
                 cacCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCt                                                                           >  1:474699/1‑100 (MQ=255)
                   cccGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTg                                                                          >  2:1015667/1‑99 (MQ=255)
                     cgcgCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAAc                                                                       >  1:1006116/1‑100 (MQ=255)
                      gcgcCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAAcc                                                                      <  1:617427/100‑1 (MQ=255)
                            ccaGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGt                                                                <  2:780304/100‑1 (MQ=255)
                                cacaAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTgaga                                                            <  1:533221/100‑1 (MQ=255)
                                  caAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAAt                                                          >  1:1009316/1‑100 (MQ=255)
                                  caAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAAt                                                          <  1:360146/100‑1 (MQ=255)
                                     aaCCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGAt                                                       >  1:597976/1‑100 (MQ=255)
                                     aaCCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGAt                                                       <  2:755008/100‑1 (MQ=255)
                                       ccGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATcc                                                     >  2:158275/1‑100 (MQ=255)
                                         gCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATcccc                                                   <  2:1285223/100‑1 (MQ=255)
                                          ccccAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCt                                                  <  1:1287913/100‑1 (MQ=255)
                                             cAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCa                                               <  1:369140/100‑1 (MQ=255)
                                                     gATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTcc                                       <  1:439759/100‑1 (MQ=255)
                                                     gATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTcc                                       >  2:14789/1‑100 (MQ=255)
                                                     gATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTcc                                       >  2:64235/1‑100 (MQ=255)
                                                     gATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTcc                                       >  1:212740/1‑100 (MQ=255)
                                                      aTCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCg                                      <  1:1021859/100‑1 (MQ=255)
                                                        cGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGc                                    <  2:395178/100‑1 (MQ=255)
                                                         gCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCa                                   <  1:30984/100‑1 (MQ=255)
                                                         gCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCa                                   <  2:96709/100‑1 (MQ=255)
                                                           cATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGTACCGAGAGTGAGAATCATCCCCCGCACCACTTCCGGc                                    >  2:736172/1‑97 (MQ=255)
                                                           cATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCa                                   >  2:1095455/1‑98 (MQ=255)
                                                            aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGAc                                >  1:512563/1‑100 (MQ=255)
                                                            aTCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGAc                                >  2:169023/1‑100 (MQ=255)
                                                             tCCAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACg                               >  1:796386/1‑100 (MQ=255)
                                                               cAGATAGCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGca                             >  2:326396/1‑100 (MQ=255)
                                                                     gCCGTTACCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAg                       >  2:1295955/1‑100 (MQ=255)
                                                                           aCCAGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATgg                 <  2:10402/100‑1 (MQ=255)
                                                                              aGCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGa              >  1:1127995/1‑100 (MQ=255)
                                                                               gCAACATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGat             >  2:1097558/1‑100 (MQ=255)
                                                                                  aCATTTTGCCCGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATAtt          >  2:819508/1‑100 (MQ=255)
                                                                                          cccGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACCAAAGCAATGGTGATATTGCGCACAt  >  1:1198114/1‑100 (MQ=255)
                                                                                          cccGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACAt  >  2:666930/1‑100 (MQ=255)
                                                                                          cccGACTGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACAt  >  1:322946/1‑100 (MQ=255)
                                                                                                |                                                                                             
GTCGCCATATGTTTGTACACCCGCGCCACCAGCACAAAACCGCCCCAGATGGTGATCGCCATCCAGATAGCCGTTACCAGCAACATTTTGCCCGACAGCATCCACAGATACGCCGCTGAACCGATAGTGAGAATGATCCCCTGCACCAGTTCCGGCAGACGCACAAAAGCAATGGTGATATTGCGCACAT  >  NC_000913/2308086‑2308275

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