Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 218,030 A→G 8.7% Y6C (TAC→TGC)  yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported218,0300AG 8.7% 125.5 / 3.3 46Y6C (TAC→TGC) yafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/17);  new base G (2/2);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAC  >  Exported/217960‑218099
                                                                      |                                                                     
aTCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                                                       >  1:383285/1‑71 (MQ=255)
aTCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                                                       >  1:3480905/1‑71 (MQ=255)
      aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGTACGTGGTAATGAATGAAGACGACTACaaaa                                                                  >  1:2242711/1‑70 (MQ=255)
      aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                                  >  1:1292249/1‑70 (MQ=255)
      aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                                  <  1:3385195/70‑1 (MQ=255)
      aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa                                                                  <  1:2723146/70‑1 (MQ=255)
         cACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt                                                                 >  1:2429495/1‑68 (MQ=255)
          aCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATaaga                                                             <  1:1578328/71‑1 (MQ=255)
           cTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGCCGACTACAAAATaag                                                              >  1:349178/1‑69 (MQ=255)
            tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                                                       >  1:544938/1‑59 (MQ=255)
            tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                                                       >  1:2053138/1‑59 (MQ=255)
              gCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAagaaga                                                          <  1:3267037/70‑1 (MQ=255)
              gCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                                                         >  1:3335552/1‑71 (MQ=255)
               cACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                                                         >  1:1323149/1‑70 (MQ=255)
               cACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                                                         >  1:2942886/1‑70 (MQ=255)
               cACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAgg                                                        <  1:25306/71‑1 (MQ=255)
                   tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTa                                                                       >  1:2129699/1‑52 (MQ=255)
                       gaATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaa                                                      >  1:2931290/1‑65 (MQ=255)
                       gaATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgca                                                 <  1:302336/70‑1 (MQ=255)
                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag                                             <  1:1231078/71‑1 (MQ=255)
                             gCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                          <  1:880816/71‑1 (MQ=255)
                              cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACg                                                   >  1:2291139/1‑61 (MQ=255)
                                    aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    >  1:533800/1‑70 (MQ=255)
                                    aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                    >  1:530657/1‑70 (MQ=255)
                                        gTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTc                                   <  1:1043071/67‑1 (MQ=255)
                                        gTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAgg                                >  1:2221498/1‑70 (MQ=255)
                                         tCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                              >  1:2323309/1‑71 (MQ=255)
                                         tCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAatta                                        <  1:2190656/61‑1 (MQ=255)
                                          cTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcagca                                              >  1:2149157/1‑54 (MQ=255)
                                           tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                            >  1:1771386/1‑71 (MQ=255)
                                           tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                            >  1:2858918/1‑71 (MQ=255)
                                           tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                            >  1:447723/1‑71 (MQ=255)
                                            ggACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa                                                       <  1:2026663/43‑1 (MQ=255)
                                                    tAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag                                                <  1:114624/42‑1 (MQ=255)
                                                    tAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAtt                   >  1:361684/1‑71 (MQ=255)
                                                    aaatgaatgAAGACGACTACAAAATAAGAAGAGGAAAc                                                    >  1:1901886/2‑38 (MQ=255)
                                                    aaatgaatgAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                           <  1:3374404/46‑1 (MQ=255)
                                                     aatgaatgAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                           >  1:2683627/1‑62 (MQ=255)
                                                     aatgaatgAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTg                  >  1:304912/1‑71 (MQ=255)
                                                      atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTg                  >  1:418109/1‑70 (MQ=255)
                                                      atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                 >  1:870195/1‑71 (MQ=255)
                                                              aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                    <  1:2903562/60‑1 (MQ=255)
                                                               gacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTAttt                                     <  1:1918273/42‑1 (MQ=255)
                                                               gacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTAttt                                     <  1:2059739/42‑1 (MQ=255)
                                                                 cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAatta                                        <  1:809062/37‑1 (MQ=37)
                                                                  gacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATAttt      >  1:1931122/1‑70 (MQ=255)
                                                                  gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt     >  1:2770476/1‑71 (MQ=255)
                                                                  gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt     >  1:2775026/1‑71 (MQ=255)
                                                                   acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTa            <  1:1174980/63‑1 (MQ=255)
                                                                    cTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTg    <  1:99112/70‑1 (MQ=255)
                                                                     tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAat          >  1:1919138/1‑63 (MQ=255)
                                                                      aCAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc  >  1:564082/1‑70 (MQ=255)
                                                                      aCAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc  >  1:988763/1‑70 (MQ=255)
                                                                      |                                                                     
ATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAC  >  Exported/217960‑218099

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: