Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,386,365 T→A F169Y (TTC→TAC)  yfbL → putative M28A family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,3650TA98.0% 169.2 / ‑3.5 51F169Y (TTC→TAC) yfbLputative M28A family peptidase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (26/24);  total (27/24)
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.58e-01

GTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAC  >  NC_000913/2386271‑2386459
                                                                                              |                                                                                              
gtgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTc                                                                                           <  2:976029/100‑1 (MQ=255)
gtgGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTc                                                                                           >  1:1283917/1‑100 (MQ=255)
  ggCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGt                                                                                         >  2:1072814/1‑100 (MQ=255)
   gCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTa                                                                                        <  1:1310937/100‑1 (MQ=255)
    cAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGATCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAg                                                                                       <  2:174974/100‑1 (MQ=255)
    cAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACAGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAg                                                                                       <  2:175675/100‑1 (MQ=255)
              cGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATg                                                                             >  2:442288/1‑100 (MQ=255)
               gAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATgg                                                                            <  2:1283912/100‑1 (MQ=255)
                aaCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATggg                                                                           <  2:1481855/100‑1 (MQ=255)
                 aCTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATgggg                                                                          <  2:1249494/100‑1 (MQ=255)
                  cTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGa                                                                         >  1:1287069/1‑100 (MQ=255)
                    ggCACGTTTGTTACATCAGCACGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAgc                                                                       >  2:291905/1‑100 (MQ=255)
                     gCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAgcg                                                                      <  2:1212148/100‑1 (MQ=255)
                        cGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCgg                                                                   >  1:714621/1‑100 (MQ=255)
                         gTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGt                                                                  <  1:645228/100‑1 (MQ=255)
                          tttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtg                                                                 <  2:134001/100‑1 (MQ=255)
                          tttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtg                                                                 <  2:737559/100‑1 (MQ=255)
                           ttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgc                                                                >  2:692243/1‑100 (MQ=255)
                           ttGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgc                                                                <  2:801306/100‑1 (MQ=255)
                            tGTTACATCAGCAAGGACCGAAAACAGGCGTGCAGCTGGTGGCCTATGCGTCGGAAGAGCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgca                                                               <  2:778932/100‑1 (MQ=255)
                             gTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcat                                                              >  2:477018/1‑100 (MQ=255)
                               tACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGtgcatgc                                                            >  1:783105/1‑100 (MQ=255)
                                  aTCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGc                                                         <  1:1245248/100‑1 (MQ=255)
                                        aaGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCt                                                   >  2:1350981/1‑100 (MQ=255)
                                          gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGTGCGGTGCATGCAGCTTCGCTTg                                                 >  2:513063/1‑100 (MQ=255)
                                          gTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTg                                                 <  2:1473836/100‑1 (MQ=255)
                                            aCCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAg                                               >  1:200665/1‑100 (MQ=255)
                                              cGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCg                                             <  2:1027290/100‑1 (MQ=255)
                                               gAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGt                                            >  1:637517/1‑100 (MQ=255)
                                               gAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGAGGATCGCGGTGCATGCAGCTTCGCTTGAGCGt                                            >  2:1024834/1‑100 (MQ=255)
                                                      ggCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGa                                     >  2:1165784/1‑100 (MQ=255)
                                                        cGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGaaa                                   <  2:776185/100‑1 (MQ=255)
                                                            cAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTaa                               >  1:395502/1‑100 (MQ=255)
                                                                 ggTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                          <  1:1376651/100‑1 (MQ=255)
                                                                 ggTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                          >  1:380172/1‑100 (MQ=255)
                                                                 ggTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATa                          <  2:251036/100‑1 (MQ=255)
                                                                  gTCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGGGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAg                         <  2:783101/100‑1 (MQ=255)
                                                                   tCGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGc                        >  2:1355850/1‑100 (MQ=255)
                                                                    cGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCa                       >  2:703752/1‑100 (MQ=255)
                                                                    cGCCTATGCGTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCa                       >  2:1194267/1‑100 (MQ=255)
                                                                     gCCTATGCTTCGGAAGAACCGCCCTATTTTCTTAGCGATGAAATGTGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCAc                      <  2:1208711/100‑1 (MQ=255)
                                                                             ggcGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatg              >  1:1060231/3‑100 (MQ=255)
                                                                             gTCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatg              <  2:379896/100‑1 (MQ=255)
                                                                              tCGGAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAgatga             >  1:376181/1‑100 (MQ=255)
                                                                                ggAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGAtt           >  1:1370897/1‑100 (MQ=255)
                                                                                ggAAGAACCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGAtt           <  1:721507/100‑1 (MQ=255)
                                                                                     aaCCGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTa      >  2:1260230/1‑100 (MQ=255)
                                                                                      aCCGCCCTACTTTCGTAGCGATGAAATGGGGAGTGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTAt     <  1:513064/100‑1 (MQ=255)
                                                                                       ccGCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTAtt    >  2:1237341/1‑100 (MQ=255)
                                                                                         gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc  >  2:223789/1‑100 (MQ=255)
                                                                                         gCCCTACTTTCGTAGCGATGAAATGGGGAGCGCGGTACATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAc  >  2:1442850/1‑100 (MQ=255)
                                                                                              |                                                                                              
GTGGCAGGATTACTCGAACTGGCACGTTTGTTACATCAGCAAGTACCGAAAACAGGCGTGCAGCTGGTCGCCTATGCGTCGGAAGAACCGCCCTTCTTTCGTAGCGATGAAATGGGGAGCGCGGTGCATGCAGCTTCGCTTGAGCGTCCAGTGAAATTAATGATAGCACTGGAGATGATTGGCTATTAC  >  NC_000913/2386271‑2386459

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: