Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,486,331 C→T 31.6% P327P (CCG→CCA nuoN ← NADH:ubiquinone oxidoreductase, membrane subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,486,3310CT31.6% 69.7 / 38.0 80P327P (CCG→CCAnuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base C (15/39);  new base T (7/18);  total (22/58)
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.54e-01

AGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGACGATGCCAG  >  minE/1486263‑1486394
                                                                    |                                                               
aGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATc                                                                 >  1:815670/1‑69 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:2981004/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:1785829/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:1789374/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:1765770/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:3425503/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:382401/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:1736760/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:390162/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:1863837/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:2015344/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:265347/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:2229970/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:829550/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:128819/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:983495/69‑1 (MQ=255)
  ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                                                               <  1:2337284/69‑1 (MQ=255)
        caAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATa                                                           >  1:2341722/1‑67 (MQ=255)
        caAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATAcc                                                         >  1:1775464/1‑69 (MQ=255)
        caAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATAcc                                                         >  1:516716/1‑69 (MQ=255)
        caAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATAcc                                                         >  1:721254/1‑69 (MQ=255)
                    gACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGCCATCgg                                                               <  1:3464/51‑1 (MQ=255)
                      cacCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:351012/68‑1 (MQ=255)
                      cacCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:315940/68‑1 (MQ=255)
                      cacCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:2984052/68‑1 (MQ=255)
                      cacCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:2940157/68‑1 (MQ=255)
                      cacCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:1509149/68‑1 (MQ=255)
                      cacCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:2638094/68‑1 (MQ=255)
                      cacCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:939873/68‑1 (MQ=255)
                       accgaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGc                                            <  1:1019412/67‑1 (MQ=255)
                         cgaccgacAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAGTGAAAGcatc                                         >  1:1811262/1‑68 (MQ=255)
                         cgaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAGTGAAAGcatc                                         >  1:1042090/1‑68 (MQ=255)
                         cgaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAGTGAAAGcatc                                         >  1:3337339/1‑68 (MQ=255)
                         cgaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAGTGAAAGcatc                                         >  1:554576/1‑68 (MQ=255)
                         cgaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAGTGAAAGcatc                                         >  1:709385/1‑68 (MQ=255)
                         cgaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAGTGAAAGcatc                                         >  1:2259088/1‑68 (MQ=255)
                         cgaccgCCAGCAAGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAGTGAAAGcatc                                         >  1:1855480/1‑68 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:2869363/1‑69 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:359436/1‑69 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:3415101/1‑69 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:1871471/1‑69 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:743941/1‑69 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:2772058/1‑69 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:2771509/1‑69 (MQ=255)
                          gaccgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGcatcat                                       >  1:1234120/1‑69 (MQ=255)
                           accgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCatcatc                                      <  1:2149008/69‑1 (MQ=255)
                           accgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCatcatc                                      <  1:1752639/69‑1 (MQ=255)
                           accgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCatcatc                                      <  1:1837477/69‑1 (MQ=255)
                           accgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCatcatc                                      <  1:911567/69‑1 (MQ=255)
                           accgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCatcatc                                      <  1:142734/69‑1 (MQ=255)
                           accgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCatcatc                                      <  1:3211541/69‑1 (MQ=255)
                           accgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCatcatc                                      <  1:2914722/69‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCGAGCGTCATTGGGATACCAGCCAgtga                                                <  1:679566/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:1631934/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:1665669/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:453661/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:1748749/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:897121/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:1043650/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:1439900/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:801216/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:777245/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:348246/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:2006058/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:2051537/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:2128440/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:2562807/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:2591723/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:277232/57‑1 (MQ=255)
                             cgCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATTGGGATACCAGCCAgtga                                                <  1:2700952/57‑1 (MQ=255)
                                   cACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGtcat                              <  1:1688249/69‑1 (MQ=255)
                                   cACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGtcat                              <  1:225859/69‑1 (MQ=255)
                                                 cGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCCTCATCACCGTCATCACTGCCGCGAgg                 <  1:3239007/68‑1 (MQ=255)
                                                 cGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAgg                 <  1:2614607/68‑1 (MQ=255)
                                                  gATAAAGCCCAGCGTCATGGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAgg                 <  1:2379509/67‑1 (MQ=255)
                                                        gCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGAc          <  1:2376528/68‑1 (MQ=255)
                                                        gCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGAc          <  1:1572077/68‑1 (MQ=255)
                                                        gCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGAc          <  1:594996/68‑1 (MQ=255)
                                                                tCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGACGATGCCAg  >  1:1447005/1‑68 (MQ=255)
                                                                tCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGACGATGCCAg  >  1:2402149/1‑68 (MQ=255)
                                                                tCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGACGATGCCAg  >  1:1920881/1‑68 (MQ=255)
                                                                    |                                                               
AGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGACAGCATCATCACCGTCATCACTGCCGCGAGGATCGGACGATGCCAG  >  minE/1486263‑1486394

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: