Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,384,511 T→C 17.6% L350L (CTA→CTG metE ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,384,5110TC17.6% 244.1 / 42.8 102L350L (CTA→CTGmetE5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (40/44);  new base C (5/13);  total (45/57)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

CTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGGTCGA  >  minE/2384444‑2384572
                                                                   |                                                             
ctGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAgg                                                            >  1:1599981/1‑71 (MQ=255)
      gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                      >  1:2436248/1‑71 (MQ=255)
      gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                      >  1:1880534/1‑71 (MQ=255)
      gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                      >  1:3786761/1‑71 (MQ=255)
      gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                      >  1:3858875/1‑71 (MQ=255)
      gagCTGCCGTATCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAg                                                      >  1:472278/1‑71 (MQ=255)
           gCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaaa                                                 >  1:1520747/1‑71 (MQ=255)
             cGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAcc                                               <  1:291584/71‑1 (MQ=255)
              gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaaa                                                 >  1:3581046/1‑68 (MQ=255)
              gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                                >  1:259854/1‑69 (MQ=255)
              gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                                >  1:1334222/1‑69 (MQ=255)
              gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                                >  1:3576565/1‑69 (MQ=255)
              gtgtCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAc                                                >  1:2318466/1‑69 (MQ=255)
              gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                              <  1:2837048/68‑1 (MQ=255)
              gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                              <  1:354384/68‑1 (MQ=255)
              gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                              <  1:2054522/68‑1 (MQ=255)
              gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                              <  1:1143566/68‑1 (MQ=255)
              gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                              <  1:3006751/68‑1 (MQ=255)
              gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                              <  1:679723/68‑1 (MQ=255)
              gtatCACCACTGTTCAACGCATCGCGCAGCAATGCCAGTTCATGACATTTTTGCAGGGCGAAGGCAAACCa                                              <  1:2547815/68‑1 (MQ=255)
                 tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                                           <  1:3208637/71‑1 (MQ=255)
                 tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                                           <  1:247978/71‑1 (MQ=255)
                 tCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCt                                           <  1:3268882/71‑1 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                         <  1:2972727/71‑1 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                         >  1:2531894/1‑71 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                         >  1:1407526/1‑71 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                         >  1:116741/1‑71 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                         >  1:503024/1‑71 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtt                                          <  1:777133/70‑1 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtt                                          >  1:2934076/1‑70 (MQ=255)
                   accacTGTTCAGCGCATCGCGCAGCAGTCCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                         <  1:513235/71‑1 (MQ=255)
                    ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt                                        <  1:3164271/71‑1 (MQ=255)
                    ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt                                        <  1:1608458/71‑1 (MQ=255)
                    ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt                                        <  1:2835817/71‑1 (MQ=255)
                    ccacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt                                        <  1:917543/71‑1 (MQ=255)
                     cacTGTTCAGCGCATCGCGCAGCAGTGCCAGTTGATGGCATTTTTGTAGGGCGAAGGCAAACCAGCttt                                         <  1:1661228/69‑1 (MQ=255)
                       cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGGaa                                                  <  1:1769867/58‑1 (MQ=255)
                       cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                                                  <  1:3650997/58‑1 (MQ=255)
                       cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                                                  <  1:570389/58‑1 (MQ=255)
                       cTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa                                                  <  1:552048/58‑1 (MQ=255)
                         gTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCa                                              >  1:3826583/1‑60 (MQ=255)
                         gTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCAttt                                   >  1:2579579/1‑71 (MQ=255)
                          ttCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGc                                                    >  1:3507569/1‑53 (MQ=255)
                          ttCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTc                                  >  1:3027212/1‑71 (MQ=255)
                            cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                                >  1:1091904/1‑71 (MQ=255)
                            cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                                >  1:1541899/1‑71 (MQ=255)
                            cAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                                >  1:2418600/1‑71 (MQ=255)
                             aGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                                <  1:1219284/70‑1 (MQ=255)
                             aGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTg                                <  1:3886369/70‑1 (MQ=255)
                                  aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:4128196/71‑1 (MQ=255)
                                  aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:1516029/71‑1 (MQ=255)
                                  aTCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:2269720/71‑1 (MQ=255)
                                   tCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                                    <  1:4097144/60‑1 (MQ=255)
                                   tCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                                    <  1:58292/60‑1 (MQ=255)
                                     gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGa                        >  1:1824521/1‑70 (MQ=255)
                                     gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAc                       >  1:1871080/1‑71 (MQ=255)
                                     gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAc                       >  1:3129345/1‑71 (MQ=255)
                                     gcgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAc                       >  1:2452422/1‑71 (MQ=255)
                                      cgcAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAcg                      <  1:3668847/71‑1 (MQ=255)
                                       gcagcaGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGAcg                      <  1:1446628/70‑1 (MQ=255)
                                         agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:595498/64‑1 (MQ=255)
                                         agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:3724605/64‑1 (MQ=255)
                                         agcagTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:2472017/64‑1 (MQ=255)
                                            agTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTc                >  1:2332054/1‑71 (MQ=255)
                                             gTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:3293700/60‑1 (MQ=255)
                                                ccAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCAc             >  1:3628290/1‑70 (MQ=255)
                                                   gTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                                    <  1:1079116/44‑1 (MQ=255)
                                                   gTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATc                            <  1:3588507/52‑1 (MQ=255)
                                                   gTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTc         <  1:3546531/71‑1 (MQ=255)
                                                    ttCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:2919233/53‑1 (MQ=255)
                                                    ttCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:2733722/53‑1 (MQ=255)
                                                    ttCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:3513786/53‑1 (MQ=255)
                                                    ttCATGACATTTTTGCAGGGCGAAGGGAAACCAGCTTTTCACtt                                   <  1:313001/44‑1 (MQ=37)
                                                    ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                                   <  1:826648/44‑1 (MQ=38)
                                                    ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                                   <  1:1579154/44‑1 (MQ=38)
                                                    ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                                   <  1:3960536/44‑1 (MQ=38)
                                                    ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                                   <  1:4089148/44‑1 (MQ=38)
                                                    ttCATGACATTTTTGCAGGGCGAAGGCAAACCAGCTTTTCACtt                                   <  1:482984/44‑1 (MQ=38)
                                                      cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACt                                    <  1:915259/41‑1 (MQ=255)
                                                      cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg      >  1:3328887/1‑71 (MQ=255)
                                                      cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg      >  1:577107/1‑71 (MQ=255)
                                                      cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg      >  1:1275993/1‑71 (MQ=255)
                                                      cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg      >  1:3260301/1‑71 (MQ=255)
                                                      cATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg      >  1:3582965/1‑71 (MQ=255)
                                                      cATGGCATTTTTGTAGGGCGAAGGCAAACAAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAgg      >  1:34398/1‑71 (MQ=255)
                                                          gCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAg       >  1:1334390/1‑66 (MQ=255)
                                                          gCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAg       >  1:739878/1‑66 (MQ=255)
                                                          gCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAg       >  1:3576373/1‑66 (MQ=255)
                                                              ttttGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTc                <  1:2025098/53‑1 (MQ=255)
                                                               tttGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGGTCGa  >  1:3482368/1‑66 (MQ=255)
                                                               tttGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGGTCGa  >  1:2382712/1‑66 (MQ=255)
                                                               tttGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGGTCGa  >  1:738791/1‑66 (MQ=255)
                                                               tttGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGGTCGa  >  1:1635108/1‑66 (MQ=255)
                                                                  gTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:1308815/39‑1 (MQ=255)
                                                                  gTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:4028440/39‑1 (MQ=255)
                                                                  gTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:1600948/39‑1 (MQ=255)
                                                                  gTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:595694/39‑1 (MQ=255)
                                                                  gTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:1270719/39‑1 (MQ=255)
                                                                  gTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:3539545/39‑1 (MQ=255)
                                                                  gTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCaa                          <  1:1937842/39‑1 (MQ=255)
                                                                   |                                                             
CTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTTCACTTCTGCATCAAGACGCGTTTCCACGCTCAGGTCGA  >  minE/2384444‑2384572

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: