Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,380,909 C→A 14.7% F198L (TTC→TTA ompG → outer membrane porin

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,380,9090C→A14.7% 2.0 ‑2.2 ‑0.027594ompGouter membrane porin
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (12/11):  total (12/15)
Fisher's exact test strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.87e-01

GCCTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTC  >  REL606/1380874‑1380944
                                   |                                   
gCCTTGCGAGTGAACTATTATCTCGAGCGCGGCTTc                                     >  1:5287815/1‑36 (MQ=255)
 ccTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCa                                    >  1:5874935/1‑36 (MQ=255)
  cTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCaa                                   <  1:1380968/36‑1 (MQ=255)
  cTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCaa                                   <  1:3738880/36‑1 (MQ=255)
     gCGAGTGAACTATTATCTCGAGCGCGGCTTCAatat                                >  1:5421782/1‑36 (MQ=255)
      cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATg                               <  1:5851210/36‑1 (MQ=255)
      cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATg                               <  1:3850412/36‑1 (MQ=255)
       gAGTGAACTATTATCTCGAGCGCGGCTTCAATATgg                              <  1:5540923/36‑1 (MQ=255)
       gAGTGAACTATTATCTCGAGCGCGGCTTCAATATgg                              >  1:1935733/1‑36 (MQ=255)
        aGTGAACTATTATCTCGAGCGCGGCTTCAATCTGga                             >  1:2708406/1‑36 (MQ=255)
         gTGAACTATTATCTCGAGCGCGGCTTCAATATGgac                            <  1:2071815/36‑1 (MQ=255)
          tGAACTATTATCTCGAGCGCGGCTTCAATATGgacg                           <  1:4918757/36‑1 (MQ=255)
           gAACTATTATCTCGAGCGCGGCTTCAATATGgacga                          <  1:6169989/36‑1 (MQ=255)
             aCTATTATCTCGAGCGCGGCTTCAATATGGACGACa                        <  1:5298926/36‑1 (MQ=255)
               tattatCTCGAGCGCGGCTTCAATATGGACGACAGc                      <  1:5436288/36‑1 (MQ=255)
                attatCTCGAGCGCGGCTTCAATATGGACGACAGcc                     >  1:3648200/1‑36 (MQ=255)
                  tatCTCGAGCGCGGCTTCAATATGGACGACAGCCGc                   <  1:1718834/36‑1 (MQ=255)
                  tatCTCGAGCGCGGCTTCAATATGGACGACAGCCGc                   >  1:1170128/1‑36 (MQ=255)
                    tctcGAGCGCGGCTTCAATATGGACGACAGCCGCaa                 >  1:3311216/1‑36 (MQ=255)
                     ctcGAGCGCGGCTTCAATATGGACGACAGCCGCAAt                >  1:4021804/1‑36 (MQ=255)
                        gAGCGCGGCTTCAATATGGACGACAGCCGCAATAAc             >  1:2191471/1‑36 (MQ=255)
                        gAGCGCGGCTTCAATATGGACGACAGCCGCAATAAc             >  1:6005704/1‑36 (MQ=255)
                        gAGCGCGGCTTCAATATGGAAGACAGCCGCAATaag             >  1:4462954/1‑35 (MQ=255)
                        aaGCGCGGCTTAAATATGGACAACAGCCGCAATAAc             <  1:1203030/35‑1 (MQ=37)
                          gggcgGCTTAAATATGAACGACAGCCGCAATAACgg           <  1:5876834/34‑1 (MQ=37)
                           cgcgGCTTAAATATGGACGACAGCCGCAATAACGGt          <  1:345294/36‑1 (MQ=255)
                             cgGCTTCAATATGGACGACAGCCGCAATAACGGTGa        >  1:1230581/1‑36 (MQ=255)
                             cgGATTAAATAGGGACGCCAGCCGCAATAACGGTGa        <  1:3843543/36‑1 (MQ=19)
                                   cAATATGGACGAAAGCCGCAATAACGGTGAGttttt  >  1:4567906/1‑35 (MQ=255)
                                   |                                   
GCCTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTC  >  REL606/1380874‑1380944

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG