Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 704,276 T→C 9.1% L123L (CTA→CTG ymcA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE704,2760TC 9.1% 153.8 / 6.8 55L123L (CTA→CTGymcAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/22);  new base C (3/2);  total (31/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 = 9.43e-01

GATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGC  >  minE/704215‑704346
                                                             |                                                                      
gATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTaa                                                               <  1:3910614/71‑1 (MQ=255)
gATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTaa                                                               >  1:102284/1‑71 (MQ=255)
 aTACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAAc                                                              >  1:1144146/1‑71 (MQ=255)
       aCCGTCGAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACcttctt                                                        >  1:532811/1‑71 (MQ=255)
       aCCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCACACGCCACTTGCGGTAGCCAGTa                                                                >  1:3084523/1‑63 (MQ=255)
       aCCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACcttc                                                          <  1:3864751/69‑1 (MQ=255)
        ccGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGc                                                                     <  1:4122212/57‑1 (MQ=255)
          gTCAAACAGACCCGTACCAGCAATATAACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTcc                                                      >  1:3779271/1‑70 (MQ=255)
             aaaCAGCCCCGTACCGCCGATATCCCGCGCACCAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCcac                                                   >  1:3502649/1‑70 (MQ=255)
             aaaCAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca                                                    >  1:1997378/1‑69 (MQ=255)
             aaaCAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca                                                    >  1:2456672/1‑69 (MQ=255)
              aaCAGCCCCGTACCGCCGATATCCCGCGCACCAACCGCCACTTGCGGCAGCCAGt                                                                 <  1:1429573/55‑1 (MQ=255)
               aCAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTaa                                                               >  1:622001/1‑56 (MQ=255)
                cAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGc                                               >  1:1491974/1‑71 (MQ=255)
                        tACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGt                                                                 <  1:1116272/45‑1 (MQ=255)
                        tACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTATCCTTCTTCCCACAGGCGAAGTTTa                                       <  1:2250353/71‑1 (MQ=255)
                          ccAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaa                                      >  1:347774/1‑70 (MQ=255)
                          ccAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaa                                      >  1:2822752/1‑70 (MQ=255)
                          ccAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaa                                      >  1:2687091/1‑70 (MQ=255)
                            aGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAg                                           <  1:3394286/63‑1 (MQ=255)
                             gCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATc                                  >  1:3844021/1‑71 (MQ=255)
                               aaTATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCg                                 <  1:1360885/70‑1 (MQ=255)
                                atatCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGaa                               >  1:558766/1‑71 (MQ=255)
                                 tatCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAAt                              <  1:923859/71‑1 (MQ=255)
                                   tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACctt                                                           >  1:441709/1‑40 (MQ=255)
                                   tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca                                                    <  1:3641498/47‑1 (MQ=255)
                                   tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCg                                              <  1:3827924/53‑1 (MQ=255)
                                   tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGt                                          >  1:822092/1‑57 (MQ=255)
                                   tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaaaa                                    <  1:2337680/63‑1 (MQ=255)
                                   tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATg                             >  1:2516100/1‑70 (MQ=255)
                                     aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAAcc                                                             >  1:448107/1‑36 (MQ=255)
                                     aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAtt                          >  1:3129330/1‑71 (MQ=255)
                                     aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAtt                          >  1:3979769/1‑71 (MQ=255)
                                      cGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAttt                         <  1:1233849/71‑1 (MQ=255)
                                       gTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTg                        >  1:1693256/1‑71 (MQ=255)
                                       gTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTg                        >  1:491916/1‑71 (MQ=255)
                                          tACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca                                                    >  1:3164502/1‑40 (MQ=255)
                                          tACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTc                      <  1:650526/70‑1 (MQ=255)
                                          tACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCt                     >  1:4029543/1‑71 (MQ=255)
                                            ccAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCg                                              <  1:3318279/44‑1 (MQ=255)
                                             cAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGttt                                        <  1:2605070/49‑1 (MQ=16)
                                             cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGa                                <  1:2675394/57‑1 (MQ=255)
                                             cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg                  <  1:266453/71‑1 (MQ=255)
                                              aaaCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCt                     <  1:3525145/67‑1 (MQ=255)
                                               accGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGaa                               >  1:611367/3‑56 (MQ=255)
                                               accGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGa                                >  1:1838457/3‑55 (MQ=255)
                                               aaCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTa                >  1:666771/1‑71 (MQ=255)
                                                aCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt               <  1:574692/71‑1 (MQ=255)
                                                aCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt               <  1:751886/71‑1 (MQ=255)
                                                    cACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt               >  1:995098/1‑67 (MQ=255)
                                                     aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGt                 <  1:1677144/64‑1 (MQ=255)
                                                     aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg          <  1:194918/71‑1 (MQ=255)
                                                     aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg          <  1:3033067/71‑1 (MQ=255)
                                                      cTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa         >  1:4055974/1‑71 (MQ=255)
                                                         gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTa                >  1:1718357/1‑61 (MQ=255)
                                                             tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcgccgc  >  1:2069970/1‑71 (MQ=255)
                                                             |                                                                      
GATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGC  >  minE/704215‑704346

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: