Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 218,049 C→T 9.1% N12N (AAC→AAT yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported218,0490CT 9.1% 165.1 / 6.3 55N12N (AAC→AATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/17);  new base T (2/3);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCA  >  Exported/217983‑218116
                                                                  |                                                                   
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:2858918/1‑71 (MQ=255)
                    tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                             >  1:447723/1‑71 (MQ=255)
                    tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                             >  1:1771386/1‑71 (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:2059739/42‑1 (MQ=255)
                                        gacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTAttt                                                      <  1:1918273/42‑1 (MQ=255)
                                          cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAatta                                                         <  1:809062/37‑1 (MQ=37)
                                           gacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATAttt                       >  1:1931122/1‑70 (MQ=255)
                                           gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                      >  1:2775026/1‑71 (MQ=255)
                                           gacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                      >  1:2770476/1‑71 (MQ=255)
                                            acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTa                             <  1:1174980/63‑1 (MQ=255)
                                             ctgcAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGa                    >  1:3232581/4‑71 (MQ=255)
                                             cTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTg                     <  1:99112/70‑1 (MQ=255)
                                              tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAat                           >  1:1919138/1‑63 (MQ=255)
                                               aCAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                   >  1:988763/1‑70 (MQ=255)
                                               aCAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                   >  1:564082/1‑70 (MQ=255)
                                                  aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCa                                                   >  1:2034124/1‑35 (MQ=255)
                                                  aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCa                                                   >  1:3544056/1‑35 (MQ=255)
                                                   aaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt               >  1:2258372/1‑70 (MQ=255)
                                                   aaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt               >  1:1186757/1‑70 (MQ=255)
                                                   aaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATg              <  1:2824569/71‑1 (MQ=255)
                                                     tAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgag            <  1:3443837/71‑1 (MQ=255)
                                                       agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                            <  1:1975509/53‑1 (MQ=255)
                                                         aagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg        <  1:277144/71‑1 (MQ=255)
                                                          agagGAAACGCAGCAGAATTATTTTCAGGGATACGTCACATTGCTATTAATATTTTGACGAATGAGAAgg        <  1:3413778/70‑1 (MQ=255)
                                                             ggAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAtt    >  1:2187636/1‑71 (MQ=255)
                                                             ggAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt     >  1:2769298/1‑70 (MQ=255)
                                                              gAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAtt    >  1:1591808/1‑70 (MQ=255)
                                                               aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                   >  1:3228906/1‑38 (MQ=255)
                                                               aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                   >  1:3127976/1‑38 (MQ=255)
                                                               aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt     <  1:688660/68‑1 (MQ=255)
                                                               aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCa  <  1:803445/71‑1 (MQ=255)
                                                                aaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGaa          >  1:3737457/1‑62 (MQ=255)
                                                                aaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc   >  1:345098/1‑69 (MQ=255)
                                                                aaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc   >  1:2253703/1‑69 (MQ=255)
                                                                aaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCa  >  1:1010457/1‑70 (MQ=255)
                                                                 aCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAtt    >  1:1661006/1‑67 (MQ=255)
                                                                  |                                                                   
GAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCA  >  Exported/217983‑218116

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: