Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,486,279 T→C 16.7% T345A (ACA→GCA)  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,2790TC16.7% 139.9 / 21.8 60T345A (ACA→GCA) nuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base T (29/21);  new base C (6/4);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GCGCAGGTAGTAGTAGAGGCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGA  >  minE/1486211‑1486348
                                                                    |                                                                     
gcgcAGGTAGTAGTAGAGGCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCt                                                                     >  1:2129947/1‑71 (MQ=255)
            gtagAGGCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACaccgaccg                                                         >  1:2914931/1‑71 (MQ=255)
            gtagAGGCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGCCTGGACaccgacc                                                          <  1:1228133/70‑1 (MQ=255)
            gtagAGGCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGCCTGGACaccgacc                                                          <  1:2823627/70‑1 (MQ=255)
               gagGCCGATTGCCGAACCGACAACCACGGCACCAACCAGCCACCACAAGTGTGTCTGGACACCGACCGc                                                        >  1:627803/1‑69 (MQ=255)
                agGCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAg                                                     <  1:2666978/71‑1 (MQ=255)
                  gCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCa                                                   >  1:529771/1‑71 (MQ=255)
                  gCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGCCTGGACACCGACCGCCAGCa                                                   <  1:3934344/71‑1 (MQ=255)
                       ttGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACg                                                 <  1:226868/68‑1 (MQ=255)
                        tGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAgtgtgt                                                                       >  1:1826097/1‑45 (MQ=255)
                         gCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGaa                                            >  1:1774344/1‑71 (MQ=255)
                          ccGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGa                                                                  <  1:635509/48‑1 (MQ=255)
                          ccGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGa                                                                  <  1:3543831/48‑1 (MQ=255)
                          ccGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGaa                                            <  1:3376915/70‑1 (MQ=255)
                            gAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGaa                                            <  1:3094482/68‑1 (MQ=255)
                             aaCCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTa                                        <  1:2536138/71‑1 (MQ=255)
                               ccGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCt                                                                     >  1:908416/1‑40 (MQ=255)
                               ccGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACacc                                                              >  1:563275/1‑47 (MQ=255)
                                 gACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTa                                        >  1:280294/1‑67 (MQ=255)
                                 gACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTa                                        >  1:932091/1‑67 (MQ=255)
                                 gACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTAc                                       >  1:975564/1‑68 (MQ=255)
                                    aaCCACGGCACCCACCAGCCACCACAAGTGTGCCCGGACACCGACCGCCAGCACGTAGaa                                            >  1:813647/1‑60 (MQ=255)
                                        aCGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGcc                             >  1:18005/1‑71 (MQ=255)
                                        aCGGCACCCACCAGCCACCACAAGTGTGCCTGTACACCGACCGCCAGCACGTAGAACTTACCGATAAAGcc                             >  1:724591/1‑71 (MQ=255)
                                         cGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGccc                            <  1:3934108/71‑1 (MQ=255)
                                          ggCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAg                                                     >  1:718910/1‑45 (MQ=255)
                                          ggCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAg                                                     >  1:688770/1‑45 (MQ=255)
                                          ggCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAg                                                     >  1:3379106/1‑45 (MQ=255)
                                          ggCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCa                           <  1:721180/71‑1 (MQ=255)
                                            cacccaccAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGc                         <  1:2848745/71‑1 (MQ=255)
                                               ccaccaGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTc                      <  1:3887576/71‑1 (MQ=255)
                                                 accaGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCAt                    >  1:1683474/1‑71 (MQ=255)
                                                  ccaGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATc                   >  1:3293166/1‑71 (MQ=255)
                                                  ccaGCCACAACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATa                                  >  1:4107380/1‑56 (MQ=255)
                                                   caGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCg                  >  1:355469/1‑71 (MQ=255)
                                                    aGCCACCACAAGTGTGTCTGGACACCGCCCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                 >  1:1335971/1‑71 (MQ=255)
                                                    aGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCAc                                                  >  1:604865/1‑38 (MQ=255)
                                                    aGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                 <  1:153363/71‑1 (MQ=255)
                                                    aGCCACCACAAGTGTGCCTGGACACCGACCGCCAGCACGTAGaa                                            >  1:1121583/1‑44 (MQ=255)
                                                     gCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                 >  1:3169380/1‑70 (MQ=255)
                                                     gCCACCACAAGTGTGCCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                 >  1:511478/1‑70 (MQ=255)
                                                      ccaccaCAAGTGTGTCTGGACACCGACCGCCAGCACGTAGa                                             <  1:1115956/41‑1 (MQ=255)
                                                       caccacAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGAt              <  1:2407402/71‑1 (MQ=255)
                                                        accacAAGTGTGTCTGGACACCGACCGCCAGCACGTAg                                              <  1:2508730/38‑1 (MQ=255)
                                                        accacAAGTGTGTCTGGACACCGACCGCCAGCACGTAg                                              <  1:787526/38‑1 (MQ=255)
                                                        accacAAGTGTGTCTGGACACCGACCGCCAGCACGTAg                                              <  1:1040198/38‑1 (MQ=255)
                                                        accacAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATa             >  1:3637652/1‑71 (MQ=255)
                                                        accacAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATa             >  1:746673/1‑71 (MQ=255)
                                                        accacAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATa             <  1:2156619/71‑1 (MQ=255)
                                                         ccacAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGc                         <  1:648560/58‑1 (MQ=255)
                                                         ccacAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATAc            <  1:423920/71‑1 (MQ=255)
                                                          cacaAGTGTGTCTGGACACCGACCGCCAGCACGTAGaa                                            >  1:3953375/1‑38 (MQ=255)
                                                          cacaAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATaa                                 >  1:220062/1‑49 (MQ=255)
                                                          cacaAGTGTGCCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGc                              <  1:1865452/52‑1 (MQ=255)
                                                           acaAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTa                                        >  1:3876628/1‑41 (MQ=255)
                                                           acaAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTa                                        >  1:2942743/1‑41 (MQ=255)
                                                           acaAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATa                                  >  1:364834/1‑47 (MQ=255)
                                                           acaAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATa                                  >  1:2712358/1‑47 (MQ=255)
                                                           acaAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCgg                 <  1:1064423/64‑1 (MQ=255)
                                                           acaAGTGTGCCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCg          >  1:1762521/1‑71 (MQ=255)
                                                            caAGTGTGCCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCgg         >  1:1550347/1‑71 (MQ=255)
                                                               gtgtgtCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGcc       <  1:3780797/70‑1 (MQ=255)
                                                                tgtgtCTGGACACCGACCGCCAGCACGTAGAACTTACCGATaaa                                >  1:3021954/1‑44 (MQ=255)
                                                                  tgtCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCaga    <  1:2801392/70‑1 (MQ=255)
                                                                  tgtCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCaga    <  1:2129906/70‑1 (MQ=255)
                                                                   gtCTGGACACCGACCGCCAGCACGTAGAACTTACCGATa                                  <  1:259287/39‑1 (MQ=255)
                                                                    tCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCg                  >  1:2462141/1‑54 (MQ=255)
                                                                    tCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAgaga  <  1:1358225/70‑1 (MQ=255)
                                                                    |                                                                     
GCGCAGGTAGTAGTAGAGGCCGATTGCCGAACCGACAACCACGGCACCCACCAGCCACCACAAGTGTGTCTGGACACCGACCGCCAGCACGTAGAACTTACCGATAAAGCCCAGCGTCATCGGGATACCGGCCAGAGA  >  minE/1486211‑1486348

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: