Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,045,956 C→A 29.9% L16F (TTG→TTT yghJ ← predicted inner membrane lipoprotein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,045,9560C→A29.9% 11.3 ‑2.7 ‑0.13048yghJpredicted inner membrane lipoprotein
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (8/13):  total (8/22)
Fisher's exact test strand distribution p-value = 6.65e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.73e-01

CCACCGCCATCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTT  >  REL606/3045921‑3045988
                                   |                                
ccACCGCCATCACAACCGGCTAACAGGGTTGCGCTc                                  <  1:6240348/36‑1 (MQ=255)
 cACCGCCATCACAACCGGCTAACAGGGTTGCGCTCa                                 >  1:782667/1‑36 (MQ=255)
  cccGCCATCACAACCGGCTAACAGGGTTGCGCTCaa                                <  1:7751995/35‑1 (MQ=255)
  aCCGCCATCACAACCGGCTAACAGGGTTGCGCTCaa                                <  1:214407/36‑1 (MQ=255)
   ccGCCATCACAGCCGGCTAACAGAGTTGCGCTCaaa                               <  1:279349/36‑1 (MQ=38)
   ccGCCATCACAACCGGCTAACAGTGTTGCGCTCaaa                               <  1:3618977/36‑1 (MQ=255)
   ccGCCATCACAACCGGCTAACAGGGTTGCGCTCaaa                               <  1:6587471/36‑1 (MQ=255)
    cGCCATCACAACCGGCTAACAGGGTTGCGCTCaaac                              <  1:6448156/36‑2 (MQ=255)
      ccATCACAACCGGCTAACAGGGTTGCGCTCAAAATa                            <  1:266243/36‑1 (MQ=255)
         tCACAACCGGCTAACAGGGTTGCGCTCAAAATAGcc                         <  1:5390493/36‑1 (MQ=255)
         tCACAACCGGCTAACAGGGTTGCGCTCAAAATAGcc                         >  1:3655042/1‑36 (MQ=255)
         tCACAACCGGCTAACAGGGTTGCGCTCAAAATAGcc                         <  1:764267/36‑1 (MQ=255)
         tCACAACCGGCAAACAGGGTTGCGCTCAAAATAGcc                         <  1:3866821/36‑1 (MQ=255)
            aaaCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCt                      <  1:2925647/35‑1 (MQ=255)
              aCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTaa                    <  1:7463922/36‑1 (MQ=255)
              aCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTaa                    >  1:3982595/1‑36 (MQ=255)
              aCCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTaa                    >  1:1511353/1‑36 (MQ=255)
              aCCGGATAACAGGGTTGCGCTAAAAATAGCCGCTaa                    <  1:1644932/36‑1 (MQ=37)
               ccGGCTAACAGGTTTGCGCTAAAAATAGCCGCTaaa                   <  1:862269/36‑1 (MQ=37)
               ccGGCTAACAGGGTTGCGCTCAAAATAGCCGCTaaa                   <  1:858968/36‑1 (MQ=255)
                 ggatAACAGGGTTGCGCTAAAAATAGCCGCTAAAAg                 <  1:4132944/33‑1 (MQ=37)
                   cTAACAGGGTTGCGCTAAAAATAGCCGCTAAAAGCg               <  1:1254683/36‑1 (MQ=255)
                    tAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGa              >  1:6819245/1‑36 (MQ=255)
                     aaaaGGGTTGCCCTAAAAATAGCCGCTAAAAGCGAt             <  1:362515/33‑1 (MQ=25)
                      aaaGGGTTGCGCTAAAAATAGCCGCTAAAAGCGAtt            <  1:7861976/34‑1 (MQ=37)
                        aGGGTGGCGCTAAAAATAGCCGCTAAAAGCGATTTc          <  1:3040733/36‑1 (MQ=37)
                         gggTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCt         >  1:3042496/1‑36 (MQ=255)
                             tGCGCAAAAAATAGCCGCTAAAAGCGATTTCTtata     <  1:1554402/36‑1 (MQ=37)
                               cgcTCAAAATAGCCGCTAAAAGCGATTTCTTATAtt   >  1:3326483/1‑36 (MQ=255)
                               cgcTAAAAATAGCCGCTAAAAGCGATTTCTTATAtt   <  1:7030195/36‑1 (MQ=255)
                                gcTCAAAATAGCCGCTAAAAGCGATTTCTTATAttt  >  1:3371244/1‑36 (MQ=255)
                                   |                                
CCACCGCCATCACAACCGGCTAACAGGGTTGCGCTCAAAATAGCCGCTAAAAGCGATTTCTTATATTT  >  REL606/3045921‑3045988

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG