Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,867,615 G→A T429I (ACC→ATC)  ygcU ← predicted FAD containing dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,867,6150GA98.9% 267.7 / ‑4.1 88T429I (ACC→ATC) ygcUpredicted FAD containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (41/46);  total (41/47)
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.93e-01

GTTTGCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTA  >  minE/1867547‑1867677
                                                                    |                                                              
gTTTGCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGAtt                                                              >  1:2885378/1‑71 (MQ=255)
 tttGCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGAtt                                                              >  1:2410904/1‑70 (MQ=255)
  ttGCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttc                                                            >  1:1030593/1‑71 (MQ=255)
  ttGCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttc                                                            >  1:161670/1‑71 (MQ=255)
   tGCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttc                                                            <  1:499784/70‑1 (MQ=255)
    gCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttct                                                           <  1:1480263/70‑1 (MQ=255)
     cTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttcttc                                                         <  1:1747161/71‑1 (MQ=255)
     cTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttctt                                                          <  1:1344324/70‑1 (MQ=255)
     cTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttctt                                                          <  1:1333709/70‑1 (MQ=255)
       ccAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATttc                                                            >  1:1628306/1‑66 (MQ=255)
        cAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTcaca                                                      <  1:960063/71‑1 (MQ=255)
        cAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTca                                                        >  1:2027525/1‑69 (MQ=255)
        cAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTca                                                        <  1:2891474/69‑1 (MQ=255)
        cAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTca                                                        >  1:854797/1‑69 (MQ=255)
        cAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTca                                                        >  1:843887/1‑69 (MQ=255)
        cAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTca                                                        >  1:2470520/1‑69 (MQ=255)
          gTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACAGCCGAGGCGAATGATTTCTTca                                                        >  1:271666/1‑67 (MQ=255)
           tGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAgat                                                   <  1:1394784/71‑1 (MQ=255)
             aaCGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAgatga                                                 <  1:1624226/71‑1 (MQ=255)
              aCGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAgatgat                                                >  1:3047559/1‑71 (MQ=255)
                gcgATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAgatg                                                  >  1:2437252/1‑67 (MQ=255)
                   aTGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgt                                           >  1:1670628/1‑71 (MQ=255)
                    tGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAga                                                    >  1:244135/1‑61 (MQ=255)
                     gTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgtt                                          <  1:721122/70‑1 (MQ=255)
                      tttACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgttg                                         <  1:1086561/70‑1 (MQ=255)
                      tttACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGa                                        <  1:1579292/71‑1 (MQ=255)
                      tttACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACACATGATCTTGTTGa                                        <  1:327407/71‑1 (MQ=255)
                      tttACCGATACCGTGGTGGTGCACCATCGAACCACCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGa                                        <  1:2072917/71‑1 (MQ=255)
                       ttACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgttg                                         <  1:2898252/69‑1 (MQ=255)
                       ttACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAg                                       >  1:192209/1‑71 (MQ=255)
                        tACCGATACCGTGGTGGTGCACCATCGAACCGCCGATGCGAATGATTTCTTCACAGATGATCttgtt                                          >  1:1175606/1‑67 (MQ=255)
                        tACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgtt                                          >  1:2552469/1‑67 (MQ=255)
                        tACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAg                                       <  1:2120143/70‑1 (MQ=255)
                        tACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAg                                       <  1:2582905/70‑1 (MQ=255)
                        tACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTAACAGATGATCttgttg                                         <  1:2213430/68‑1 (MQ=255)
                             aTACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAtt                                 >  1:1520887/1‑71 (MQ=255)
                             aTACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAtt                                 <  1:1824527/71‑1 (MQ=255)
                               aCCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATtgt                               >  1:1803627/1‑71 (MQ=255)
                               aCCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATtgt                               >  1:2565530/1‑71 (MQ=255)
                                ccGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATtgtg                              >  1:512879/1‑71 (MQ=255)
                                   tggtggtgCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGa                                        >  1:2312521/1‑58 (MQ=255)
                                     gtggtgCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTAc                          <  1:4879/70‑1 (MQ=255)
                                         tgCACCATCGAACCGCCGAGGCGAATGATTACTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTc                     >  1:2587855/1‑71 (MQ=255)
                                          gCACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATtg                                <  1:2762438/59‑1 (MQ=255)
                                           cACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCg                    <  1:493761/70‑1 (MQ=255)
                                           cACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCg                    <  1:908010/70‑1 (MQ=255)
                                           cACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCg                    <  1:1520673/70‑1 (MQ=255)
                                           cACCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGa                   <  1:362584/71‑1 (MQ=255)
                                            aCCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATtgtg                              >  1:2297949/1‑59 (MQ=255)
                                            aCCATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAt                  <  1:2957830/71‑1 (MQ=255)
                                             ccATCGAACCGCCGAGGCGAATGATTTCTTCACAGATTATCTTGTTGAGCGGa                                   <  1:2795076/53‑1 (MQ=255)
                                             ccATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGa                                   <  1:1439406/53‑1 (MQ=255)
                                              cATCGAACCGCCGAGGCGAATGATTTCTTCACAgatg                                                  >  1:2001174/1‑37 (MQ=255)
                                              cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCt                                              <  1:2141096/41‑1 (MQ=255)
                                              cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCgg                                    <  1:2235257/51‑1 (MQ=255)
                                              cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTAc                          <  1:2165901/61‑1 (MQ=255)
                                              cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAttt                >  1:2882396/1‑71 (MQ=255)
                                              cATCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAttt                >  1:1387738/1‑71 (MQ=255)
                                                tCGAACCGCCGAGGCGAATGATTTCTTCACAgatg                                                  <  1:2767814/35‑1 (MQ=255)
                                                tCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgttg                                         <  1:1358825/44‑1 (MQ=255)
                                                tCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAt                                  >  1:2558152/1‑51 (MQ=255)
                                                tCGAACCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTcc              <  1:76840/71‑1 (MQ=255)
                                                  gAACTGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCctc            <  1:625346/71‑1 (MQ=255)
                                                  gAACCGCCGAGGCGAATGATTTCTTCACAGATGATCttgt                                           >  1:2211538/1‑40 (MQ=255)
                                                   aaCCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAt                                  >  1:221418/1‑48 (MQ=255)
                                                   aaCCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAt                                  >  1:744029/1‑48 (MQ=255)
                                                   aaCCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCctct           >  1:1545209/1‑71 (MQ=255)
                                                    aCCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCtt          <  1:926603/71‑1 (MQ=255)
                                                    aCCGCCGAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCtt          <  1:2563037/71‑1 (MQ=255)
                                                     ccgccgAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTc         >  1:1361568/1‑71 (MQ=255)
                                                      ctccgAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTg                       <  1:2557250/54‑1 (MQ=255)
                                                       gccgAGTCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCg       >  1:393318/1‑71 (MQ=255)
                                                       gccgAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTTGATTTCCTCTTCCg       <  1:1465633/71‑1 (MQ=255)
                                                        ccgAGGCGAATGATTTCTTCACAGATTATCTTGTTGAGCGGATTGTGGTa                           <  1:652634/50‑1 (MQ=255)
                                                        ccgAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCgg      >  1:547498/1‑71 (MQ=255)
                                                        ccgAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCgg      >  1:2716134/1‑71 (MQ=255)
                                                         cgAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCgg      <  1:1996333/70‑1 (MQ=255)
                                                         cgAGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGc     >  1:2994882/1‑71 (MQ=255)
                                                          gAGGCGAATGATTTCTTCACAGATGATCTCGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCt    <  1:2025295/71‑1 (MQ=255)
                                                           aGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGa                                   <  1:1574831/39‑1 (MQ=255)
                                                           aGGCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCt    >  1:2095834/1‑70 (MQ=255)
                                                             gCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAtt                                 >  1:1217671/1‑39 (MQ=255)
                                                             gCGAATGATTTCTTCACAGATGATCTTGTTGAGCGGAtt                                 >  1:2890581/1‑39 (MQ=255)
                                                              cGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTa                           <  1:1417979/44‑1 (MQ=255)
                                                              cGAATGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTa  >  1:2889499/1‑69 (MQ=255)
                                                                aaTGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGAt                  <  1:772494/51‑1 (MQ=255)
                                                                 aTGATTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCtt   <  1:1696104/65‑1 (MQ=255)
                                                                  tGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGa                   <  1:977322/48‑1 (MQ=255)
                                                                    |                                                              
GTTTGCTCCAGTGAACGCGATGTTTACCGATACCGTGGTGGTGCACCATCGAACCGCCGAGGCGAATGGTTTCTTCACAGATGATCTTGTTGAGCGGATTGTGGTACTTGTCGATTTCCTCTTCCGGCTTA  >  minE/1867547‑1867677

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: