Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,484,894 T→C 8.9% M316V (ATG→GTG)  menF ← isochorismate synthase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,484,8940TC 8.9% 186.3 / 9.2 67M316V (ATG→GTG) menFisochorismate synthase 2
Reads supporting (aligned to +/- strand):  ref base T (24/37);  new base C (4/2);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGA  >  Exported/1484824‑1484961
                                                                      |                                                                   
aaCTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAt                                                                     <  1:672867/71‑1 (MQ=255)
aaCTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAt                                                                     <  1:1895515/71‑1 (MQ=255)
 aCTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAtt                                                                    >  1:897544/1‑71 (MQ=255)
 aCTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAt                                                                     <  1:3581774/70‑1 (MQ=255)
  cTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAttt                                                                   <  1:808499/71‑1 (MQ=255)
     aTGTAAACAGATCACATCATCCTCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAtttt                                                                  <  1:1462030/69‑1 (MQ=255)
     aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAtttt                                                                  <  1:2757203/69‑1 (MQ=255)
     aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAtttt                                                                  <  1:1448201/69‑1 (MQ=255)
     aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                >  1:764038/1‑71 (MQ=255)
     aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                >  1:733399/1‑71 (MQ=255)
     aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                >  1:2979491/1‑71 (MQ=255)
     aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTa                                                                 <  1:3665677/70‑1 (MQ=255)
     aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGAGGCGAAGATGCTGCAtttt                                                                  <  1:3080508/69‑1 (MQ=255)
      tGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTacg                                                               <  1:785080/71‑1 (MQ=255)
       gTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTacg                                                               >  1:1679906/1‑70 (MQ=255)
       gTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTacg                                                               >  1:263157/1‑70 (MQ=255)
          aaCAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAttt                                                                   <  1:3051751/63‑1 (MQ=255)
          aaCAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTacgcaga                                                           >  1:1226270/1‑71 (MQ=255)
            cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTacgcagacg                                                         <  1:695789/71‑1 (MQ=255)
            cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                >  1:3469751/1‑64 (MQ=255)
            cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                >  1:821182/1‑64 (MQ=255)
            cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                >  1:2745422/1‑64 (MQ=255)
            cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                >  1:2499814/1‑64 (MQ=255)
                tCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGt                                                     <  1:1803752/71‑1 (MQ=255)
                tCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgcag                                                      <  1:1919014/70‑1 (MQ=255)
                 cacaTCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGCCGCAGTa                                                    >  1:458836/1‑71 (MQ=255)
                    atcatcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCt                                                 <  1:1143469/71‑1 (MQ=255)
                    atcatcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCt                                                 <  1:1308160/71‑1 (MQ=255)
                        tcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTgcgg                                             <  1:2505282/71‑1 (MQ=255)
                        tcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTgcgg                                             <  1:1871909/71‑1 (MQ=255)
                        tcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTgcgg                                             <  1:3046337/71‑1 (MQ=255)
                                gttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTaaa                                      <  1:1640190/70‑1 (MQ=255)
                                 ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTaaaa                                     <  1:2440728/70‑1 (MQ=255)
                                 ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTaaaa                                     <  1:2496865/70‑1 (MQ=255)
                                 ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTaaaa                                     <  1:40118/70‑1 (MQ=255)
                                 ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAAc                                    <  1:1619856/71‑1 (MQ=255)
                                 ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAAc                                    <  1:623734/71‑1 (MQ=255)
                                 ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAAc                                    <  1:3535883/71‑1 (MQ=255)
                                 ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAAc                                    <  1:2843561/71‑1 (MQ=255)
                                  tgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACa                                   >  1:3136979/1‑71 (MQ=255)
                                  tgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACa                                   >  1:524159/1‑71 (MQ=255)
                                   gAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTa                                                                 >  1:3490981/1‑40 (MQ=255)
                                   gAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTgc                                               <  1:1875674/58‑1 (MQ=255)
                                   gAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgcag                                                      <  1:2867992/51‑1 (MQ=255)
                                    aGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACAt                                  <  1:210053/70‑1 (MQ=255)
                                    aGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACAt                                  <  1:2818457/70‑1 (MQ=255)
                                      tGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCa                               <  1:1649358/71‑1 (MQ=255)
                                       gAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTgcgg                                             >  1:436384/1‑56 (MQ=255)
                                       gAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTgcgg                                             >  1:1742661/1‑56 (MQ=255)
                                       gAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAg                              >  1:2099186/1‑71 (MQ=255)
                                        aaGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTaa                                       <  1:2072095/61‑1 (MQ=255)
                                        aaGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAg                              <  1:2673456/70‑1 (MQ=255)
                                        aaGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAg                              <  1:103333/70‑1 (MQ=255)
                                           tCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTc                          <  1:3362574/71‑1 (MQ=255)
                                            ccAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCt                         >  1:3633208/1‑71 (MQ=255)
                                             cAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCt                         <  1:3565508/70‑1 (MQ=255)
                                             cAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCt                         <  1:814115/70‑1 (MQ=255)
                                                aTACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTggg                      >  1:1780123/1‑70 (MQ=255)
                                                aTACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTggg                      >  1:1311442/1‑70 (MQ=255)
                                                aTACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGt                     <  1:3439938/71‑1 (MQ=255)
                                                  aCAGCGGCGAAGCTGCTGCATTTTACGAAGCCGCAGTACCTgc                                               >  1:943043/1‑43 (MQ=37)
                                                  aCAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTggg                      >  1:1274529/1‑68 (MQ=255)
                                                    aGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCg                 >  1:1390266/1‑71 (MQ=255)
                                                        gcgAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGcc              <  1:2570148/70‑1 (MQ=255)
                                                        gcgAAGATGCTGCATTTTACGCAGACGCAGTACCAgcgg                                             <  1:998742/39‑1 (MQ=255)
                                                          gAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTg            <  1:2802601/70‑1 (MQ=255)
                                                          gAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGt           <  1:658150/71‑1 (MQ=255)
                                                          gAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGt           <  1:2646405/71‑1 (MQ=255)
                                                            aGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTaa         >  1:1249313/1‑71 (MQ=255)
                                                            aGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaa                                       >  1:711608/1‑41 (MQ=255)
                                                            aGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaa                                       >  1:715942/1‑41 (MQ=255)
                                                            aGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTa                                        >  1:76419/1‑40 (MQ=255)
                                                               tgctgcATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCg      <  1:2454933/71‑1 (MQ=255)
                                                               tgctgcATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCg      <  1:1726280/71‑1 (MQ=255)
                                                                 ctgcATTTTACGCAGACGCAGTACCTGCGGCGGTaaa                                      >  1:2280043/1‑37 (MQ=255)
                                                                 ctgcATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGtt    <  1:2683127/71‑1 (MQ=255)
                                                                  tgcATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTg            <  1:227048/62‑1 (MQ=255)
                                                                      ttttACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGa  >  1:196974/1‑68 (MQ=255)
                                                                      |                                                                   
AACTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGA  >  Exported/1484824‑1484961

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: ‑