Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 741,873 T→C 12.2% Y124Y (TAT→TAC ycfL → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE741,8730TC12.2% 102.2 / 6.5 41Y124Y (TAT→TACycfLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/21);  new base C (3/2);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TTCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTC  >  minE/741807‑741943
                                                                  |                                                                      
ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGTCGCACAAAGTCAGACTTTATCTATATTTGt                                                                    <  1:1438884/71‑1 (MQ=255)
ttCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGt                                                                    <  1:190783/71‑1 (MQ=255)
  cggcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGt                                                                    >  1:177853/1‑69 (MQ=255)
  cggcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGt                                                                    >  1:2670828/1‑69 (MQ=255)
  cggcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTa                                                                   <  1:2700431/70‑1 (MQ=255)
    gcggTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAgg                                                                <  1:3311783/71‑1 (MQ=255)
              tGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAggg                                                               >  1:3189512/1‑62 (MQ=255)
              tGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATACTTGTAAGGGGTTAATCtt                                                      >  1:348066/1‑71 (MQ=255)
               gTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAggg                                                               <  1:996971/61‑1 (MQ=255)
               gTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCtt                                                      >  1:25091/1‑70 (MQ=255)
                tACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTtga                                                    >  1:770832/1‑71 (MQ=255)
                  cGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTtgat                                                   <  1:3828999/70‑1 (MQ=255)
                  cGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTtgat                                                   <  1:1255771/70‑1 (MQ=255)
                      aGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTaa                                                                  >  1:2766518/1‑51 (MQ=255)
                        cgcCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACaaaa                                            <  1:1216523/71‑1 (MQ=255)
                          ccAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAAt                                          <  1:936020/71‑1 (MQ=255)
                           cAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATACTTGt                                                                    <  1:2464113/44‑1 (MQ=255)
                             aTTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTaa                                                                  >  1:2321391/1‑44 (MQ=255)
                             aTTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAg                                       <  1:62784/71‑1 (MQ=255)
                                   tGGGGGCGCACAAAGTCAGACTTTATCTATACTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTa                                 >  1:3719145/1‑71 (MQ=255)
                                      gggCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGt                                                             <  1:2807340/40‑1 (MQ=255)
                                        gcgcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCt                            <  1:3327637/71‑1 (MQ=255)
                                          gcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGaa                                                          >  1:2098241/1‑39 (MQ=255)
                                          gcACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGaa                                                          >  1:251041/1‑39 (MQ=255)
                                            acaAAGTCAGACTTTATCTATATTTGTAAGGGGTGAAt                                                         <  1:1753250/38‑1 (MQ=255)
                                             caAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGAtt                       <  1:2671302/71‑1 (MQ=255)
                                              aaaGTCAGACTTTATCTATATTTGTAAGGGGTTAAt                                                         <  1:1061602/36‑1 (MQ=255)
                                              aaaGTCAGACTTTATCTATACTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTa                      >  1:3710508/1‑71 (MQ=255)
                                               aaGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTAc                     >  1:2120750/1‑71 (MQ=255)
                                               aaGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTAc                     >  1:1901098/1‑71 (MQ=255)
                                                aGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTAcc                    <  1:3370963/71‑1 (MQ=255)
                                                  tCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTc                                     <  1:219827/52‑1 (MQ=255)
                                                  tCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTg                          <  1:767438/63‑1 (MQ=255)
                                                    aGACTTTATCTATACTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCg                                    <  1:2667697/51‑1 (MQ=255)
                                                      aCTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTa                      >  1:3213847/1‑63 (MQ=255)
                                                        tttATCTATATTTGTAAGGGGTGAATCTTGATGGCAAAAATGAGTCGCTACGCCTTGATTACCGCGCTgg             >  1:655555/1‑70 (MQ=255)
                                                        tttATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTgg             >  1:2383108/1‑70 (MQ=255)
                                                          tATCTATATTTGTAAGGGGTGAATCTTGATGACaa                                              >  1:2149536/1‑35 (MQ=255)
                                                          tATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGc            <  1:3426848/69‑1 (MQ=255)
                                                          tATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGc            <  1:1850994/69‑1 (MQ=255)
                                                            tCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGAt         <  1:2086760/70‑1 (MQ=255)
                                                                tatTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGt       <  1:2056640/68‑1 (MQ=255)
                                                                  tttGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTtctc  >  1:2956081/1‑71 (MQ=255)
                                                                  |                                                                      
TTCGGCGGTAACGCTGTACGGCAGCGCCAATTTTCTGGGGGCGCACAAAGTCAGACTTTATCTATATTTGTAAGGGGTGAATCTTGATGACAAAAATGAGTCGCTACGCCTTGATTACCGCGCTGGCGATGTTTCTC  >  minE/741807‑741943

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