Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,045,069 G→T 46.7% T312N (ACC→AAC)  yghJ ← predicted inner membrane lipoprotein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,045,0690G→T46.7% 75.6 ‑0.2 ‑0.256935yghJpredicted inner membrane lipoprotein
Reads supporting (aligned to +/- strand):  new base (14/11):  ref base (15/14):  total (31/25)
Fisher's exact test strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.24e-01

CCAGTTCAAAGGTATCGATACCAAAGGAGATGGTTTCGCCCCAGCTAAAGGAAAATTCACCATTTTCC  >  REL606/3045037‑3045104
                                |                                   
ccAGTTCAAAGGTATCGATACCAAAGGAGATGGttt                                  <  1:4913273/36‑1 (MQ=255)
 cAGTTCAAAGGTATCGATACCAAAGGAGATGGTTTc                                 <  1:3366589/36‑1 (MQ=255)
  aGTTCAAAGGTATCGATACCAAAGGAGATGTTTTCg                                <  1:2997166/36‑1 (MQ=255)
   gTTCAAAGGTATCGATACCAAAGGAGATGTTTTCGc                               <  1:4850598/36‑1 (MQ=255)
   gTTCAAAGGTATCGATACCAAAGGAGATGGTTTCGc                               <  1:7055199/36‑1 (MQ=255)
   gTTCAAAGGTATCGATACCAAAGGAGATGGTTTCGc                               >  1:3514723/1‑36 (MQ=255)
    ttCAAAGGTATCGATACCACAGGAGATGGTTTctgg                              >  1:2512036/1‑33 (MQ=37)
    ttCAAAGGTATCGATACCAAAGGAGATGTTTTTGcc                              >  1:6913860/1‑36 (MQ=37)
     tCAAAGGTATCGATACCAAAGGAGATGGTTTCGccc                             <  1:2004965/36‑1 (MQ=255)
      cAAAGGTATCGATACCAAAGGAGATGTTTTCGcccc                            <  1:5830675/36‑1 (MQ=255)
      cAAAGGTATCGATACCAAAGGAGATGGTTTCGcccc                            >  1:432456/1‑36 (MQ=255)
       aaaGGTATCGATACCAAAGGAGATGTTTTCGCCCCa                           >  1:154060/1‑36 (MQ=255)
       aaaGGTATCGATACCAAAGGAGATGTTTTCGCCCCa                           <  1:7143015/36‑1 (MQ=255)
       aaaGGTATCGATACCAAAGGAGATGGTTTCGCCCCa                           <  1:1957360/36‑1 (MQ=255)
       aaaGGTATCGATACCAAAGGAGATGGTTTCGCCCCa                           >  1:7881462/1‑36 (MQ=255)
        aaGGTATCGATACCAAAGGAGATGTTTTCGCCCCAg                          >  1:5823956/1‑36 (MQ=255)
        aaGGTATCGATACCAAAGGAGATGGTTTCGCCCCAg                          >  1:5892772/1‑36 (MQ=255)
         aGGTATCGATACCAAAGGAGATGTTTTCGCCCCAGc                         >  1:7558780/1‑36 (MQ=255)
         aGGTATCGATACCAAAGGAGATGGTTTCGCCCCAGc                         <  1:336813/36‑1 (MQ=255)
         aGGTATCGATACCAAAGGAGATGGTTTCGCCCCAGc                         <  1:4599470/36‑1 (MQ=255)
          ggTATCGATCCCAAAGGAGATGGTTTCGCCCCAGCt                        <  1:2949733/36‑1 (MQ=255)
          ggTATCGATACCAAAGGAGATGGTTTCGCCCCAGCt                        >  1:6824456/1‑36 (MQ=255)
           gTATCGATACCAAAGGAGATGGTTTCGCCCCAGCTa                       >  1:2963881/1‑36 (MQ=255)
            tATCGATACCAAAGGAGATGTTTTCGCCCCAGCTaa                      <  1:2821334/36‑1 (MQ=255)
             aTCGATACCAAAGGAGATGTTTTTGCCCCAGCTaaa                     >  1:3531877/1‑36 (MQ=37)
             aTCGATACCAAAGGAGATGTTTTCGCCCCAGCTaaa                     >  1:4579195/1‑36 (MQ=255)
             aTCGATACCAAAGGAGATGGTTTCGCCCCAGCTaaa                     >  1:4251632/1‑36 (MQ=255)
             aTCGATACCAAAGGAGATGGTTTCGCCCCAGCTaaa                     <  1:4221845/36‑1 (MQ=255)
             aTCGATACCAAAGGAGATGGTTTCGCCCCAGCTaaa                     <  1:4533299/36‑1 (MQ=255)
             aTCGATACCAAAGGAGATGGTTTCGCCCCAGCTaaa                     >  1:5561446/1‑36 (MQ=255)
             aTCGATACAAAAGGAGATGTTTTCGCCCCAGCTaaa                     <  1:3159666/36‑1 (MQ=37)
              tCGATACCAAAGGAGATGTTTTCGCCCCAGCTAAAg                    <  1:5368395/36‑1 (MQ=255)
                gATACCAAAGGAGATGGTTTCGCCCCAGCTAAAGGa                  >  1:6960725/1‑36 (MQ=255)
                gATACCAAAGGAGATGGTTTCGCCCCAGCTAAAGGa                  >  1:5420904/1‑36 (MQ=255)
                 aTACCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaa                 >  1:6130641/1‑36 (MQ=255)
                 aTACCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaa                 >  1:1352055/1‑36 (MQ=255)
                  tACCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaaa                <  1:7313240/36‑1 (MQ=255)
                  tACCAAAGGAGATGTTTTCGCCCCAGCTAAAGGaaa                >  1:1087765/1‑36 (MQ=255)
                  tACCAAAGGAGATGGTTTCGCCCCAGCTAAAGGaaa                <  1:5622201/36‑1 (MQ=255)
                   aCCAAAGGAGATGCTTTCGCCACAGCTAAAGGaaaa               >  1:460120/1‑36 (MQ=37)
                   aCCAAAGGAGATGATTTCGCCCCAGCTAAAGGaaaa               >  1:5884483/1‑36 (MQ=255)
                    ccAAAGGAGATGGTTTCGCCCCAGCTAAAGGAAAAt              >  1:1414248/1‑36 (MQ=255)
                    ccAAAGGAGATGGTTTCGCCCCAGCTAAAGGAAAAt              <  1:365115/36‑1 (MQ=255)
                     cAAAGGAGATGTTTTCGCCCCAGCTAAAGGAAAAtt             >  1:4129377/1‑36 (MQ=255)
                       aaGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTCa           >  1:7891893/1‑36 (MQ=255)
                       aaGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTCa           >  1:7018695/1‑36 (MQ=255)
                       aaGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTCa           >  1:4662239/1‑36 (MQ=255)
                       aaGGAGATGGTTTCGCCCCAGCTAAAGGAAAATTCa           >  1:6634368/1‑36 (MQ=255)
                        aGGAGATGTTTTCTCCCCAGCTAAAGGAAAATACAc          >  1:185300/1‑36 (MQ=25)
                        aGGAGATGTTTTCGCCCCAGCTAAAGGAAAATTCAc          <  1:6724225/36‑1 (MQ=255)
                        aGGAGATGGTTTCGCCCCAGCTAAAGGAAAATTCAc          <  1:6883251/36‑1 (MQ=255)
                          gagaTGTTTTCGCCCCAGCTAAAGGAAAATTCACCa        <  1:3956510/36‑1 (MQ=255)
                           agaTGTTTTCGCCCCAGCTAAAGGAAAATTCACCAt       <  1:6071353/36‑1 (MQ=255)
                           agaTGGTTTCGCCCCAGCTAAAGGAAAATTCACCAt       >  1:2302397/1‑36 (MQ=255)
                            gaTGGTTTCGCCCCAGCTAAAGGAAAATTCACCAtt      >  1:6171162/1‑36 (MQ=255)
                            gaTGGTTTCGCCCCAGCTAAAGGAAAATTCACCAtt      <  1:4658622/36‑1 (MQ=255)
                               ggTTTCGCCCCAGCTAAAGGAAAATTCACCATTTTc   <  1:237747/36‑1 (MQ=255)
                               ggTTTCGCCCCAGCTAAAGGAAAATTCACCATTTTc   <  1:4421466/36‑1 (MQ=255)
                                gTTTCGCCCCAGCTAAAGGAAAATTCACCATTTTcc  <  1:3539430/36‑1 (MQ=255)
                                |                                   
CCAGTTCAAAGGTATCGATACCAAAGGAGATGGTTTCGCCCCAGCTAAAGGAAAATTCACCATTTTCC  >  REL606/3045037‑3045104

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG