Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,570,800 C→T 5.8% I235I (ATC→ATT)  yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,570,8000CT 5.8% 227.9 / 2.9 69I235I (ATC→ATT) yibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/26);  new base T (2/2);  total (41/28)
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.16e-01

GCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGG  >  minE/2570732‑2570869
                                                                    |                                                                     
gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCg                                                                     >  1:3770683/1‑71 (MQ=255)
 cAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                                                                    <  1:607265/71‑1 (MQ=255)
  aaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                                                                   <  1:1629912/71‑1 (MQ=255)
   aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                                                                   >  1:3893099/1‑70 (MQ=255)
   aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                                                                   >  1:3217210/1‑70 (MQ=255)
   aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                                                                   >  1:1231582/1‑70 (MQ=255)
     atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                                                                >  1:3573855/1‑71 (MQ=255)
      tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                                                                >  1:746901/1‑70 (MQ=255)
      tataACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCggg                                                               >  1:90012/1‑71 (MQ=255)
       ataACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCggg                                                               >  1:4124072/1‑70 (MQ=255)
        taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGccc                                                                  <  1:24079/66‑1 (MQ=255)
        taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                                             <  1:125495/71‑1 (MQ=255)
        taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                                             >  1:145549/1‑71 (MQ=255)
        taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGt                                                              <  1:647673/70‑1 (MQ=255)
        taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGt                                                              <  1:173734/70‑1 (MQ=255)
          aCCTTGAGCAGACTGTCATTCGCGCTCCGAGCAATGGCTACGTTACTCAGGTACTGATCC‑CCCAGTTACa                                                           <  1:1564847/70‑1 (MQ=255)
          aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                                                           <  1:364936/71‑1 (MQ=255)
          aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                                                           <  1:369893/71‑1 (MQ=255)
             ttGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                                            >  1:91420/1‑67 (MQ=255)
                     aCTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCa                                                                 <  1:352150/54‑1 (MQ=255)
                     aCTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGc                                                <  1:3178123/71‑1 (MQ=255)
                         tCATTCGCTCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCtt                                             <  1:2425951/70‑1 (MQ=255)
                         tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                                                                       >  1:1157155/1‑44 (MQ=255)
                          cATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                                                                       >  1:781991/1‑43 (MQ=255)
                          cATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGc                                           >  1:1136439/1‑71 (MQ=255)
                           aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATcccc                                                                    <  1:1082063/45‑3 (MQ=255)
                           aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGccc                                                                  <  1:2686883/47‑1 (MQ=255)
                            ttCGCGCGCCGAGCAATGGCTACGTTACTCAGTTACTGATCCGCCCAGGTa                                                             <  1:2577331/51‑1 (MQ=255)
                            ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                                                                   <  1:3357757/45‑1 (MQ=255)
                            ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                 >  1:893123/1‑63 (MQ=255)
                              cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAg                                                    >  1:1584515/1‑58 (MQ=255)
                              cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:3002171/1‑59 (MQ=255)
                              cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt                                       >  1:3700257/1‑71 (MQ=255)
                              cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt                                       >  1:3876204/1‑71 (MQ=255)
                                 gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCg                                    <  1:1270308/71‑1 (MQ=255)
                                     cGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCCTCCg                                <  1:2447190/71‑1 (MQ=255)
                                       aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                                             <  1:417223/40‑1 (MQ=255)
                                       aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                                            >  1:3563087/1‑41 (MQ=255)
                                       aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGt                              <  1:1065265/71‑1 (MQ=255)
                                       aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGt                              <  1:456533/71‑1 (MQ=255)
                                         cAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:2046087/1‑48 (MQ=255)
                                         cAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:381454/1‑48 (MQ=255)
                                         cAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:1657286/1‑48 (MQ=255)
                                         cAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:3440506/1‑48 (MQ=255)
                                         cAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt                                       <  1:2007969/60‑1 (MQ=255)
                                             ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATgg                         >  1:2618532/1‑70 (MQ=255)
                                             ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                        >  1:1853236/1‑71 (MQ=255)
                                              gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                <  1:3670160/62‑1 (MQ=255)
                                              gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                       >  1:879129/1‑71 (MQ=255)
                                              gCTACATTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                       >  1:898074/1‑71 (MQ=255)
                                               cTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCt                      <  1:3171041/71‑1 (MQ=255)
                                                 aCGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                                                      >  1:418328/1‑37 (MQ=255)
                                                 aCGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                                                      >  1:3672401/1‑37 (MQ=255)
                                                  cGTTACTCATGTACTGATCCGCCCAGTTACATACGCAGCTGCCTTGCCGCt                                       <  1:2603132/51‑1 (MQ=255)
                                                  cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                    <  1:3110670/70‑1 (MQ=255)
                                                     tACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                >  1:1107052/1‑55 (MQ=255)
                                                     tACTCAAGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa                            >  1:3583/1‑59 (MQ=255)
                                                      aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc               >  1:182383/1‑71 (MQ=255)
                                                      aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATccc               >  1:2161884/1‑71 (MQ=255)
                                                           ggTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCAt                  >  1:1712282/1‑63 (MQ=255)
                                                            gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGc          >  1:2296211/1‑70 (MQ=255)
                                                             tACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                    >  1:1665525/1‑59 (MQ=255)
                                                              aCTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATgg                         >  1:1208620/1‑53 (MQ=255)
                                                              aCTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaa        >  1:619582/1‑70 (MQ=255)
                                                              aCTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaa        >  1:2762689/1‑70 (MQ=255)
                                                               cTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                                          <  1:1082422/35‑1 (MQ=255)
                                                               cTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                                          <  1:2132146/35‑1 (MQ=255)
                                                                 gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa     >  1:3632640/1‑70 (MQ=255)
                                                                 gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc    >  1:1082184/1‑71 (MQ=255)
                                                                 gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAAc    >  1:931954/1‑71 (MQ=255)
                                                                  aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACg   <  1:2514797/71‑1 (MQ=255)
                                                                  aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACg   >  1:3659507/1‑71 (MQ=255)
                                                                    ccGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACgg  >  1:2772083/1‑70 (MQ=255)
                                                                    |                                                                     
GCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGG  >  minE/2570732‑2570869

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: ‑