Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,591,415 G→T K203N (AAG→AAT yfeS → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,591,4150GT98.2% 158.3 / ‑3.6 55K203N (AAG→AATyfeSconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (23/31);  total (23/32)
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.99e-01

TACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAG  >  minE/1591356‑1591481
                                                           |                                                                  
tACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCaa                                                         >  1:3042096/1‑71 (MQ=255)
tACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCaa                                                         >  1:2472876/1‑71 (MQ=255)
  cGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAAt                                                        >  1:1588715/1‑70 (MQ=255)
       cTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                                    >  1:1079883/1‑69 (MQ=255)
        tGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAg                                                 >  1:1174246/1‑71 (MQ=255)
           aTTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAgc                                                <  1:1601716/69‑1 (MQ=255)
               ggTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAgc                                                >  1:106730/1‑65 (MQ=255)
                 tCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtcatc                                        <  1:2330334/71‑1 (MQ=255)
                 tCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtcat                                         >  1:1100386/1‑70 (MQ=255)
                   aGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                                  >  1:102929/1‑43 (MQ=255)
                   aGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtcatca                                       >  1:2089805/1‑70 (MQ=255)
                      ttAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCtt                                   <  1:2850477/71‑1 (MQ=255)
                       tAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTa                                  >  1:2270493/1‑71 (MQ=255)
                       aaaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGcat                                   >  1:2933962/2‑68 (MQ=255)
                        aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTa                                  <  1:646976/70‑1 (MQ=255)
                        aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTa                                  <  1:2546270/70‑1 (MQ=255)
                        aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTa                                  <  1:2779084/70‑1 (MQ=255)
                           gTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                                  >  1:2865133/1‑35 (MQ=255)
                              aTGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCt                                    >  1:2261093/1‑62 (MQ=255)
                              aTGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGc                              <  1:1186060/68‑1 (MQ=255)
                              aTGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaa                           <  1:2585721/71‑1 (MQ=255)
                              aTGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaa                            >  1:142248/1‑70 (MQ=255)
                                 ggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTa                                                  <  1:2715510/45‑1 (MQ=255)
                                 ggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTa                                                  <  1:2828138/45‑1 (MQ=255)
                                 ggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTa                                                  <  1:114502/45‑1 (MQ=255)
                                 ggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAt                        <  1:2121851/71‑1 (MQ=255)
                                  gTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTa                                                  <  1:734439/44‑1 (MQ=255)
                                  gTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAtt                       <  1:2704798/71‑1 (MQ=255)
                                   tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCtt                                   <  1:1767549/58‑1 (MQ=255)
                                    aaaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttct                     >  1:1930439/1‑71 (MQ=255)
                                     aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttctt                    >  1:902088/1‑71 (MQ=255)
                                     aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttct                     <  1:608253/70‑1 (MQ=255)
                                     aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttct                     <  1:611769/70‑1 (MQ=255)
                                      aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc                   <  1:1826004/71‑1 (MQ=255)
                                          tcCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc                   <  1:98702/67‑1 (MQ=255)
                                            ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAAt                <  1:545244/68‑1 (MQ=255)
                                            ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAAt                <  1:929364/68‑1 (MQ=255)
                                            ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAAt                <  1:2353562/68‑1 (MQ=255)
                                            ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAAt                <  1:1123477/68‑1 (MQ=255)
                                            ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATggg             <  1:2822895/71‑1 (MQ=255)
                                            ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATggg             >  1:2570542/1‑71 (MQ=255)
                                            ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATggg             >  1:1561292/1‑71 (MQ=255)
                                             cATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtcat                                         >  1:2743022/1‑42 (MQ=255)
                                             cATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATgggg            >  1:2459427/1‑71 (MQ=255)
                                             cATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATggg             >  1:1123557/1‑70 (MQ=255)
                                              aTAAACTTAAAAATATGGGACTCAATGCACTAGCGCGt                                            <  1:3047187/38‑1 (MQ=255)
                                               tAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtc                                           <  1:1045840/38‑1 (MQ=255)
                                               tAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTa                                  <  1:1904426/47‑1 (MQ=255)
                                                ataCTTAAAAAGATGGGACTCAATGCACTAGCGCGTCTTAAGCTTa                                  <  1:25735/44‑1 (MQ=255)
                                                aaaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAtt                       >  1:2969358/1‑57 (MQ=255)
                                                aaaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAAt                >  1:2482766/1‑64 (MQ=255)
                                                 aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCt                                    >  1:1101410/1‑43 (MQ=255)
                                                   cTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACgg      <  1:194531/71‑1 (MQ=255)
                                                      aaaaaTATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTc    <  1:1856938/70‑1 (MQ=255)
                                                       aaaaTATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGt     <  1:1829063/68‑1 (MQ=255)
                                                        aaaTATGGGACTCAATGCACTAGCGCGTCATCAGc                                     >  1:3023474/1‑35 (MQ=255)
                                                        aaaTATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAg  <  1:2082720/70‑1 (MQ=255)
                                                           |                                                                  
TACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAG  >  minE/1591356‑1591481

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