Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 264,098 G→T L50L (CTC→CTA ykfH ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913264,0980GT97.7% 139.9 / NA 43L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (22/20);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT  >  NC_000913/264036‑264162
                                                              |                                                                
gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                             >  1:1115576/1‑68 (MQ=255)
  gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTc                                                           >  2:722735/1‑68 (MQ=255)
     gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                        >  2:328598/1‑68 (MQ=255)
       tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCg                                                      >  2:918037/1‑68 (MQ=255)
          tcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGata                                                   >  2:710792/1‑68 (MQ=255)
           caggcagtTTTTGTGCCCATTGCGGTCCTGGCGGCGGTCGTCGTATACCGCTAGGACTCCGTCGatat                                                  <  2:245/68‑1 (MQ=255)
           caggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGatat                                                  <  2:743540/68‑1 (MQ=255)
           caggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACGCCGGCGatat                                                  <  2:1278720/68‑1 (MQ=255)
            aggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                 <  2:661265/68‑1 (MQ=255)
               caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTTGGACTCCGGCGATATCCgg                                              >  1:187646/1‑68 (MQ=255)
               caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                              <  1:1112014/68‑1 (MQ=255)
               caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                              >  2:416083/1‑68 (MQ=255)
               caggAAAAGGGCACATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                              >  1:190189/1‑68 (MQ=25)
                aggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGa                                             <  2:224241/68‑1 (MQ=255)
                    ataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                         <  1:1161654/68‑1 (MQ=255)
                    ataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                         >  2:904089/1‑68 (MQ=255)
                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                     <  2:454999/68‑1 (MQ=255)
                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                     >  1:647446/1‑68 (MQ=255)
                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                     >  1:455001/1‑68 (MQ=255)
                         gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  2:848873/1‑68 (MQ=255)
                         gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  2:1303912/1‑68 (MQ=255)
                         gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    <  2:684289/68‑1 (MQ=255)
                              ttGCGGTCCTGGCGGCTTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                               <  1:1308513/68‑1 (MQ=255)
                              ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGACCGTCCCTGCCGgt                               >  2:1113430/1‑68 (MQ=255)
                              ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                               <  2:82472/68‑1 (MQ=255)
                               tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTTGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                              <  2:187646/68‑1 (MQ=255)
                               tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                              <  2:1189166/68‑1 (MQ=255)
                               tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                              <  1:1052440/68‑1 (MQ=255)
                                gCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                             >  1:894621/1‑68 (MQ=255)
                                gCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                             <  2:366532/68‑1 (MQ=255)
                                     ccTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                        >  2:1214993/1‑68 (MQ=255)
                                        ggcggcGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAg                     <  2:190189/68‑1 (MQ=255)
                                          cggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAggc                   >  2:89693/1‑66 (MQ=255)
                                          cggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                   >  1:1116503/1‑68 (MQ=255)
                                             cgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt                <  2:1228752/68‑1 (MQ=255)
                                             cgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt                <  1:424465/68‑1 (MQ=255)
                                                  ggCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGt           <  1:416085/68‑1 (MQ=255)
                                                   gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg          <  1:1064443/68‑1 (MQ=255)
                                                     gTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGcc        >  1:416954/1‑68 (MQ=255)
                                                        taCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACt     >  2:955405/1‑68 (MQ=255)
                                                        taCCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACt     <  1:1303913/68‑1 (MQ=255)
                                                         aCCGCTGGGACTCCGGCGGAATCCGGAACGTCCCTGCCGGTGTAGTGACAGCCGCTGCCGTGCCACTg    >  2:1328341/1‑68 (MQ=17)
                                                         aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg    >  2:1086008/1‑68 (MQ=255)
                                                           cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  2:647444/68‑1 (MQ=255)
                                                           cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  2:1170385/68‑1 (MQ=255)
                                                           cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  1:604751/68‑1 (MQ=255)
                                                           cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  1:1278709/68‑1 (MQ=255)
                                                              |                                                                
GTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT  >  NC_000913/264036‑264162

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