Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,930,575 C→A 100% F137L (TTC→TTA yecE → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,930,5750CA99.1% 401.2 / ‑3.3 114F137L (TTC→TTAyecEhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (66/47);  total (67/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 114 total aligned reads displayed.

TCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCG  >  REL606/1930527‑1930622
                                                |                                               
tCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGc                                               >  1:9302527/1‑51 (MQ=255)
 ccTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGcc                                              >  1:4324745/1‑51 (MQ=255)
  cTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCa                                             >  1:1209434/1‑51 (MQ=255)
  cTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCa                                             >  1:5163907/1‑51 (MQ=255)
  cTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCa                                             >  1:948318/1‑51 (MQ=255)
  cTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCa                                             >  1:1624878/1‑51 (MQ=255)
   tGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCaa                                            >  1:3779657/1‑51 (MQ=255)
   tGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCaa                                            >  1:10535505/1‑51 (MQ=255)
   tGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCaa                                            >  1:3548508/1‑51 (MQ=255)
    ggCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCaaa                                           <  1:3104070/51‑1 (MQ=255)
     gCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAg                                          >  1:1140822/1‑51 (MQ=255)
       gAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGga                                        >  1:8885184/1‑51 (MQ=255)
       gAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGga                                        >  1:7269467/1‑51 (MQ=255)
        aaTTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGgag                                       >  1:9098950/1‑51 (MQ=255)
         aTTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGgaga                                      >  1:6674697/1‑51 (MQ=255)
          tttAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGagaa                                     >  1:7726669/1‑51 (MQ=255)
           ttAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGagaag                                    >  1:9308367/1‑51 (MQ=255)
           ttAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGagaag                                    >  1:1530922/1‑51 (MQ=255)
           ttAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGagaag                                    >  1:10282327/1‑51 (MQ=255)
           ttAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGagaag                                    >  1:2437537/1‑51 (MQ=255)
            tAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGagaaga                                   >  1:2163154/1‑51 (MQ=255)
             aaTTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAagag                                  >  1:9274833/1‑51 (MQ=255)
             aaTTATGGCGTGGAAGCCCGCCATCCACAGTTTTTAGCCAAAGGAGAagag                                  <  1:7994386/51‑1 (MQ=255)
              aTTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAgg                                 >  1:845737/1‑51 (MQ=255)
              aTTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAgg                                 >  1:2110292/1‑51 (MQ=255)
              aTTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAgg                                 >  1:3766675/1‑51 (MQ=255)
               ttATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGa                                >  1:10416645/1‑51 (MQ=255)
               ttATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGa                                >  1:6612600/1‑51 (MQ=255)
                 aTGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAAc                              >  1:1098971/1‑51 (MQ=255)
                 aTGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAAc                              >  1:5311852/1‑51 (MQ=255)
                  tGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACa                             >  1:4654521/1‑51 (MQ=255)
                  tGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACa                             >  1:3424478/1‑51 (MQ=255)
                  tGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACa                             <  1:5303097/51‑1 (MQ=255)
                   ggCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAagag                                  <  1:840708/45‑1 (MQ=255)
                   ggCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACaa                            >  1:3139777/1‑51 (MQ=255)
                   ggCGGGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACaa                            <  1:4512715/51‑1 (MQ=255)
                    gCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACaaa                           >  1:1933207/1‑51 (MQ=255)
                     cGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAAc                          >  1:3180518/1‑51 (MQ=255)
                     cGTGGAAGCCCGCCATCCACAGTTTTTAGCCAAAGGagaa                                     <  1:6094642/40‑1 (MQ=38)
                      gTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAAc                          <  1:1057190/50‑1 (MQ=255)
                      gTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                         >  1:1649530/1‑51 (MQ=255)
                      gTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                         >  1:88126/1‑51 (MQ=255)
                      gTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                         >  1:9344296/1‑51 (MQ=255)
                        ggAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCt                       >  1:3750777/1‑51 (MQ=255)
                        ggAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCt                       <  1:5296738/51‑1 (MQ=255)
                        ggAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCt                       <  1:10411402/51‑1 (MQ=255)
                        ggAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCt                       <  1:10086200/51‑1 (MQ=255)
                         gAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCtt                      >  1:8622124/1‑51 (MQ=255)
                         gAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCtt                      >  1:10458002/1‑51 (MQ=255)
                          aaGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTa                     >  1:5646000/1‑51 (MQ=255)
                          aaGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTa                     <  1:1792728/51‑1 (MQ=255)
                          aaGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTa                     >  1:1210196/1‑51 (MQ=255)
                          aaGCCCGCCATCCACAGTTTTTAGCCAAAGgag                                       <  1:826088/33‑1 (MQ=37)
                           aGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTaa                    <  1:5039295/51‑1 (MQ=255)
                           aGTCCGCCATCCACAGATTTTAGCCAAAGGAGAAGAGGAACAAACGCTTaa                    >  1:1501979/1‑51 (MQ=255)
                             tCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATc                  >  1:7174632/1‑51 (MQ=255)
                              cccccATCCACAGTTTTTAGCCAAAGGAGAAGAGGaa                               <  1:5396703/34‑1 (MQ=37)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACa                             <  1:1256110/39‑1 (MQ=39)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAAc                          <  1:692103/42‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                         <  1:1054766/43‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGc                        <  1:17570/44‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTaa                    <  1:2162602/48‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAAt                   <  1:6887680/49‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcg                 >  1:9556811/1‑51 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcg                 <  1:10882019/51‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcg                 <  1:2268880/51‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATc                  <  1:7145365/50‑1 (MQ=255)
                              ccGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATc                  <  1:2076660/50‑1 (MQ=255)
                               cGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcgc                >  1:53958/1‑51 (MQ=255)
                                gCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcgcg               >  1:865060/1‑51 (MQ=255)
                                gCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcgcg               >  1:3578690/1‑51 (MQ=255)
                                 ccATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACa                             <  1:5783322/36‑1 (MQ=255)
                                 ccATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcgcg               <  1:7530685/50‑1 (MQ=255)
                                  cATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACaa                            <  1:7398749/36‑1 (MQ=255)
                                  cATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACa                             <  1:4091962/35‑1 (MQ=255)
                                  cATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGt             >  1:5711310/1‑51 (MQ=255)
                                   aTCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTa                     <  1:7528174/42‑1 (MQ=255)
                                    tCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGttt           <  1:4083396/51‑1 (MQ=255)
                                    tCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGttt           <  1:434666/51‑1 (MQ=255)
                                    tCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGttt           >  1:2381290/1‑51 (MQ=255)
                                    tCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGttt           >  1:2455657/1‑51 (MQ=255)
                                     ccACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTa          <  1:5057520/51‑1 (MQ=255)
                                     ccACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTa          <  1:9559860/51‑1 (MQ=255)
                                     ccACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTa          <  1:2879794/51‑1 (MQ=255)
                                     ccACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTa          >  1:1330540/1‑51 (MQ=255)
                                     ccACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTa          <  1:5740563/51‑1 (MQ=255)
                                      cacaGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTAc         <  1:1472190/51‑1 (MQ=255)
                                       acaGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACa        >  1:1948013/1‑51 (MQ=255)
                                        caGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACAt       >  1:7283157/1‑51 (MQ=255)
                                        caGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACAt       <  1:5926830/51‑1 (MQ=255)
                                        caGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACAt       <  1:7913308/51‑1 (MQ=255)
                                        caGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACAt       >  1:3629890/1‑51 (MQ=255)
                                         aGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATc      <  1:8108839/51‑1 (MQ=255)
                                            ttttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgc   >  1:6548656/1‑51 (MQ=255)
                                            ttttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgc   >  1:6584832/1‑51 (MQ=255)
                                            ttttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgc   >  1:8690607/1‑51 (MQ=255)
                                             tttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgcg  >  1:5767284/1‑51 (MQ=255)
                                             tttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgcg  >  1:726384/1‑51 (MQ=255)
                                             tttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgcg  >  1:4669892/1‑51 (MQ=255)
                                             tttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgcg  >  1:89750/1‑51 (MQ=255)
                                                |                                               
TCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCG  >  REL606/1930527‑1930622

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