Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,005,755 G→C 7.1% intergenic (‑161/‑363) pntA ← / → ydgH pyridine nucleotide transhydrogenase, alpha subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,005,7550GC 7.1% 258.4 / 7.8 84intergenic (‑161/‑363)pntA/ydgHpyridine nucleotide transhydrogenase, alpha subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (47/31);  new base C (2/4);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TAAAATAAATTAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAG  >  minE/1005685‑1005818
                                                                      |                                                               
tAAAATAATTAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 >  1:327077/1‑70 (MQ=255)
 aaaaTAAATTAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:2953343/70‑1 (MQ=255)
 aaaaTAAATTAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:772792/70‑1 (MQ=255)
          tAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGt                                                                <  1:1448057/62‑1 (MQ=255)
              aaCAAACGCCCTATAACGTACTGAAAATTATGCCTGTGAGCTAGCGCCAAAAATCAGTATTTCGGCGTg                                                     <  1:1767510/69‑1 (MQ=255)
               aCAAACGCCCTATAACTTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAAc                                                  <  1:2668686/71‑1 (MQ=255)
                cAAACGCCCTATAACGTACTGAAAATTATGTCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACt                                                 >  1:1052079/1‑71 (MQ=255)
                cAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACt                                                 >  1:2276066/1‑71 (MQ=255)
                  aaCGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTa                                                               >  1:798552/1‑55 (MQ=255)
                   aCGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGtt                                                     <  1:2519115/64‑2 (MQ=255)
                    cGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCg                                               >  1:471676/1‑69 (MQ=255)
                    cGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGAGCCAAAAATCAGTATTTCGGCGTGAACTCGc                                              >  1:1308692/1‑70 (MQ=255)
                       ccTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTa                                                               >  1:1252668/1‑50 (MQ=255)
                       ccTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCaaaa                                          <  1:2851099/71‑1 (MQ=255)
                       ccTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCaaaa                                          <  1:254056/71‑1 (MQ=255)
                       ccTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCaaaa                                          <  1:899598/71‑1 (MQ=255)
                       ccTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCaaaa                                          <  1:1132493/71‑1 (MQ=255)
                         tataACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:2600662/46‑1 (MQ=255)
                         tataACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:3225206/46‑1 (MQ=255)
                         tataACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:58217/46‑1 (MQ=255)
                          ataACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAata                                        >  1:1616043/1‑70 (MQ=255)
                          ataACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAAGCAGTATTTCGGCGTGAACTCGCAAAatat                                       >  1:2958394/1‑71 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCaaaa                                          <  1:676598/67‑1 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAatat                                       >  1:663737/1‑70 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAatat                                       >  1:918048/1‑70 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAatat                                       >  1:1324093/1‑70 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAata                                        >  1:842081/1‑69 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATAtt                                      <  1:2225457/71‑1 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCACTAt                                                              <  1:2128834/47‑1 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCACTAt                                                              <  1:482110/47‑1 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCACTAt                                                              <  1:2125061/47‑1 (MQ=255)
                           taACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCACTAt                                                              <  1:2898071/47‑1 (MQ=255)
                           taACGTACTGAAAATAATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAata                                        >  1:2150000/1‑69 (MQ=255)
                            aaCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTa                                                               >  1:3152261/1‑45 (MQ=255)
                            aaCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTc                                                >  1:1240581/1‑60 (MQ=255)
                            aaCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTc                                                >  1:165816/1‑60 (MQ=255)
                            aaCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTc                                                >  1:723943/1‑60 (MQ=255)
                            aaCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTc                                                >  1:2205340/1‑60 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCaaaa                                          >  1:128138/1‑65 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTaa                                    >  1:1920525/1‑71 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTaa                                    >  1:422457/1‑71 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTaa                                    >  1:93476/1‑71 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTaa                                    >  1:2028042/1‑71 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTaa                                    >  1:2914229/1‑71 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCACTAtca                                                            >  1:119636/1‑45 (MQ=255)
                             aCGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCACTAtca                                                            >  1:415643/1‑45 (MQ=255)
                              cGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAAc                                                  <  1:1385220/56‑1 (MQ=255)
                              cGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAAc                                   <  1:900850/71‑1 (MQ=255)
                              cGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAAc                                   >  1:2655556/1‑71 (MQ=255)
                                 aCTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:1115252/38‑1 (MQ=255)
                                 aCTGAAAATTATGCCTATGATCTAGCGCCAAAAATCAg                                                                 <  1:2618322/38‑1 (MQ=255)
                                  cTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGt                                                                <  1:2598103/38‑1 (MQ=255)
                                  cTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCg                                               <  1:3229112/55‑1 (MQ=255)
                                   tGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:1417437/36‑1 (MQ=255)
                                   tGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:2487071/36‑1 (MQ=255)
                                   tGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAg                                                                 <  1:3120905/36‑1 (MQ=255)
                                   tGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTg                                                     <  1:1118219/48‑1 (MQ=255)
                                   tGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAAc                                   <  1:306723/66‑1 (MQ=255)
                                    gAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTa                                                               >  1:3296434/1‑37 (MQ=255)
                                    gAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTAt                                                              >  1:3283218/1‑38 (MQ=255)
                                    gAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGaa                                                   >  1:3018393/1‑49 (MQ=255)
                                    gAAAATTATGCCTGTCATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGAtt                               >  1:42623/1‑69 (MQ=255)
                                       aaTTATGCCTGTGATCTAGCGCCAAAAATCAGTAtt                                                             >  1:1218017/1‑36 (MQ=255)
                                       aaTTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCg                                               >  1:3349092/1‑50 (MQ=255)
                                       aaTTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCg                                               >  1:1640068/1‑50 (MQ=255)
                                       aaTTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCg                                               >  1:1661099/1‑50 (MQ=255)
                                       aaTTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGcc                          >  1:2602782/1‑71 (MQ=255)
                                         ttATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCg                                               >  1:553323/1‑48 (MQ=255)
                                         ttATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGt                        <  1:2877703/71‑1 (MQ=255)
                                         ttATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGt                        <  1:2936196/71‑1 (MQ=255)
                                         ttATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAACATTAACGATTCAGCCGt                        <  1:793839/71‑1 (MQ=255)
                                           aTGCCTGTGATCTAGCGCCAAAAATCAGTGTTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGa                      >  1:1481146/1‑71 (MQ=255)
                                                tgtgATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTg                 >  1:2986122/1‑71 (MQ=255)
                                                  tgATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGt                  >  1:2777130/1‑68 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATaa           >  1:2957142/1‑68 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATa            <  1:2532129/67‑1 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATa            <  1:1460574/67‑1 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATa            <  1:1396213/67‑1 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATa            <  1:722975/67‑1 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATa            <  1:172268/67‑1 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATa            <  1:50158/67‑1 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATa            <  1:2579504/67‑1 (MQ=255)
                                                         gcgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAaca        >  1:948219/1‑71 (MQ=255)
                                                          cgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAacac       <  1:2956776/71‑1 (MQ=255)
                                                          cgcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAacac       <  1:409079/71‑1 (MQ=255)
                                                           tccAAAAATCAGTATTTCGGCGTTAACTCGCAAAATATTAACGATTCAGCCGTg                       <  1:1960605/53‑1 (MQ=255)
                                                           gcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAAc                                   >  1:239971/1‑42 (MQ=255)
                                                           gcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAAc                                   >  1:297719/1‑42 (MQ=255)
                                                           gcCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAAc                                   >  1:2886271/1‑42 (MQ=255)
                                                              aaaaaTCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCt     >  1:2528633/1‑69 (MQ=255)
                                                              aaaaaTCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCt     >  1:1764813/1‑69 (MQ=255)
                                                              aaaaaTCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCt     >  1:1801827/1‑69 (MQ=255)
                                                              aaaaaTCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCt     >  1:2566339/1‑69 (MQ=255)
                                                               aaaaTCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAg  <  1:766664/71‑1 (MQ=255)
                                                               aaaaTCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAg  <  1:9605/71‑1 (MQ=255)
                                                                      |                                                               
TAAAATAAATTAGTAACAAACGCCCTATAACGTACTGAAAATTATGCCTGTGATCTAGCGCCAAAAATCAGTATTTCGGCGTGAACTCGCAAAATATTAACGATTCAGCCGTGATAGTGGGATAAACACCTTAG  >  minE/1005685‑1005818

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: