Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,509,926 C→T 51.1% P364P (CCG→CCA yicH ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,509,9260CT51.1% ‑4.5 / 81.2 88P364P (CCG→CCAyicHconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/29);  new base T (15/30);  total (29/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGTTTT  >  minE/2509868‑2509979
                                                          |                                                     
gTTGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGta                                              <  1:3386592/68‑1 (MQ=255)
gTTGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGta                                              <  1:462253/68‑1 (MQ=255)
 ttGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGtata                                            <  1:1707016/69‑1 (MQ=255)
 ttGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGtata                                            <  1:905848/69‑1 (MQ=255)
 ttGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGtata                                            <  1:2161375/69‑1 (MQ=255)
 ttGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGtata                                            <  1:1711040/69‑1 (MQ=255)
 ttGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGtata                                            <  1:1241828/69‑1 (MQ=255)
 ttGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGtata                                            <  1:460387/69‑1 (MQ=255)
 ttGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGtata                                            <  1:1958771/69‑1 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCAGAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:3084733/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:287777/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:3044067/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:3095057/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:2765476/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:2111874/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:2431010/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:2315441/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:2284408/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:1519562/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:132504/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:742069/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:835406/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:873557/1‑68 (MQ=255)
           tttAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGtcc                                   >  1:94191/1‑68 (MQ=255)
               accaaccGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCa                               <  1:1290971/68‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:1082017/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:3239417/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:3085564/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:571287/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:3032384/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:639610/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:2924485/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:2539483/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:763844/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:1455071/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:1477621/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:206780/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGcc                         <  1:1127129/61‑1 (MQ=255)
                            cGTTTCCATCCACAACTGTTGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAATTGc                          <  1:2602243/60‑1 (MQ=255)
                             gTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCAtcatca                <  1:2789965/69‑1 (MQ=255)
                                tccatccaCAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAgg             >  1:1212361/1‑69 (MQ=255)
                                tccatccaCAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAgg             >  1:1129709/1‑69 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCCAGGTatat                                           <  1:3122499/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTttat                                           <  1:1011083/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTatat                                           <  1:2515564/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTatat                                           <  1:600991/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTatat                                           <  1:211811/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTatat                                           <  1:202747/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTatat                                           <  1:857909/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTatat                                           <  1:543741/35‑1 (MQ=255)
                                    tccaCAACTGTTGCCAGTTTTTCGGCAAGGTatat                                           <  1:2326083/35‑1 (MQ=255)
                                     ccaCAACTGTTGTCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCaa        >  1:2705480/1‑69 (MQ=255)
                                     ccaCAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAATGATCaa        >  1:498093/1‑69 (MQ=255)
                                     ccaCAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCaa        >  1:1421643/1‑69 (MQ=255)
                                     ccaCAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCaa        >  1:2690506/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:755192/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:550299/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:53286/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:2125555/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:2558819/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:1821019/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:802144/1‑69 (MQ=255)
                                      cacaACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAAc       >  1:2006287/1‑69 (MQ=255)
                                           cTGTTGCCAGTTTTTCGGCAATGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:1105808/69‑1 (MQ=255)
                                           cTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:2405867/69‑1 (MQ=255)
                                           cTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:2258927/69‑1 (MQ=255)
                                           cTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:120531/69‑1 (MQ=255)
                                           cTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:692860/69‑1 (MQ=255)
                                           cTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:1848658/69‑1 (MQ=255)
                                           cTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:609806/69‑1 (MQ=255)
                                            tGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:162141/68‑1 (MQ=255)
                                            tGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGtttt  <  1:2505471/68‑1 (MQ=255)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:701157/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:93860/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:1644997/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:488363/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:1515185/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:362494/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:1651096/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:1726314/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:3040178/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:1928340/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:2922387/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:2665125/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:2075215/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:2427239/40‑1 (MQ=37)
                                              ttGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:2210811/40‑1 (MQ=37)
                                               tGCCAGTTTTTTGGCAAGGTATATTCCAGTCCGGCAAtt                            <  1:908100/39‑1 (MQ=37)
                                                          |                                                     
GTTGCAGGCTGTTTAACCAACCGGGTGTCGTTTCCATCCACAACTGTTGCCAGTTTTTCGGCAAGGTATATTCCAGCCCGGCAATTGCCGCATCATCAAGGATCAACGTTTT  >  minE/2509868‑2509979

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