Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 980,180 A→G 60.3% D365G (GAC→GGC)  rpsA → 30S ribosomal protein S1

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606980,1800A→G60.3% 55.9 0.0 ‑0.2531094rpsA30S ribosomal protein S1
Reads supporting (aligned to +/- strand):  new base (18/14):  ref base (11/10):  total (29/24)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.70e-01

GTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACT  >  REL606/980145‑980215
                                   |                                   
gTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGa                                     >  1:5087102/1‑36 (MQ=255)
  ggCAGCAGTTCGCGGAAACCCACAACAAGGGCGAcc                                   >  1:5748648/1‑36 (MQ=255)
  ggCAGCAGTTCGCGGAAACCCACAACAAGGGCGAcc                                   >  1:5826712/1‑36 (MQ=255)
   gcagcaGTTCGCGGAAACCCACAACAAGGGCGGCCg                                  <  1:4130517/36‑1 (MQ=255)
   gcagcaGTTCGCGGAAACCCACAACAAGGGCGGCCg                                  <  1:5059063/36‑1 (MQ=255)
   gcagcaGTTCGCGGAAACCCACAACAAGGGCGAGCg                                  >  1:4678622/1‑36 (MQ=255)
     agcagTTCGCGGAAACCCACAACAAGGGCGGCCgtg                                >  1:7020381/1‑36 (MQ=255)
      gcagTTCGCGGAAACCCACAACAAGGGCGACCgtgt                               <  1:595479/36‑1 (MQ=255)
       cagTTCGCGGAAACCCACAACAAGGGCGACCGTGtt                              <  1:2535172/36‑1 (MQ=255)
       cagTTCGCGGAAACCCACAACAAGGGCGACCGTGtt                              <  1:319842/36‑1 (MQ=255)
         gTTCGCGGAAACCCACAACAAGTGCGGCCGTGTTGa                            >  1:7352656/1‑36 (MQ=37)
         gTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGa                            <  1:6348788/36‑1 (MQ=255)
          ttCGCGGAAACCCACAACAAGGGCGGCCGTGTTGaa                           >  1:1267470/1‑36 (MQ=255)
          ttCGCGGAAACCCACAACAAGGGCGGCCGTGTTGaa                           <  1:5701923/36‑1 (MQ=255)
           tCGCGGAAACCCACAACAAGGGCGGCCGTGTTGAAg                          >  1:2517229/1‑36 (MQ=255)
           tCGCGGAAACCCACAACAAGGGCGGCCGTGTTGAAg                          <  1:721061/36‑1 (MQ=255)
               ggAAACCCACAACAAGGGCGGTCGTGTTGAAGGTaa                      <  1:799523/36‑1 (MQ=37)
               ggAAACCCACAACAAGGGCGGCCGTGTTGAAGGTaa                      >  1:112366/1‑36 (MQ=255)
               ggAAACCCACAACAAGGGCGGCCGTGTTGAAGGTaa                      >  1:1319080/1‑36 (MQ=255)
               ggAAACCCACAACAAGGGCGACCGTGTTGAAGGTaa                      >  1:3446028/1‑36 (MQ=255)
               ggAAACCCACAACAAGCGCGACCGTGTTGAAGGTaa                      >  1:6004390/1‑36 (MQ=255)
                gAAACCCACAACCAGGGCGGACGTGTTGAAGGTaaa                     >  1:3901292/1‑36 (MQ=25)
                gAAACCCACAACAAGTGCGGTCGTGTTGAAGGTaaa                     <  1:6911984/36‑1 (MQ=25)
                gAAACCCACAACAAGGGCGGCCGTGTTGAAGGTaaa                     >  1:41290/1‑36 (MQ=255)
                 aaaCCCACAACAAGGGCGACCGTGTTGAAGGTaaaa                    <  1:7148543/36‑1 (MQ=255)
                  aaCCCACAACAAGGGCGGCCGTGTTGAAGGTAAAAt                   <  1:5989734/36‑1 (MQ=255)
                   aCCCACAACAAGGGCGGCCGTGTTGAAGGTAAAATc                  >  1:6219883/1‑36 (MQ=255)
                   aCCCACAACAAGGGCGGCCGTGTTGAAGGTAAAATc                  <  1:5681730/36‑1 (MQ=255)
                   aCCCACAACAAGGGCGACCGTGTTGAAGGTAAAATc                  >  1:5569054/1‑36 (MQ=255)
                    cccACAACAAGGGCGGCCGTGTTGAAGGTAAAATCa                 >  1:1501215/1‑36 (MQ=255)
                    cccACAACAAGGGCGACCGTGTTGAAGGTAAAATCa                 <  1:6021214/36‑1 (MQ=255)
                    cccACAACAAGGGCGACCGTGTTGAAGGTAAAATCa                 <  1:7676580/36‑1 (MQ=255)
                     ccACAACAAGGGCGGCCGTGTTGAAGGTAAAATCaa                <  1:3080709/36‑1 (MQ=255)
                     ccACAACAAGGGCGGCCGTGTTGAAGGTAAAATCaa                >  1:1330632/1‑36 (MQ=255)
                      cacaACAAGGGCGGCCGTGTTGAAGGTAAAATCAAg               <  1:4643220/36‑1 (MQ=255)
                      cacaACAAGGGCGACCGTGTTGAAGGTAAAATCAAg               <  1:7459167/36‑1 (MQ=255)
                       acaacaAGGGCGGCCGTGTTGAAGGTAAAATCAAGt              >  1:5170519/1‑36 (MQ=255)
                        caacaaGGGCGGCCGTGTTGAAGGTAAAATCAAGTc             >  1:3807374/1‑36 (MQ=255)
                        caacaaGGGCGGCCGTGTTGAAGGTAAAATCAAGTc             <  1:1848761/36‑1 (MQ=255)
                        caacaaGGGCGACCGTGTTGAAGGTAAAATCAAGTc             >  1:6118661/1‑36 (MQ=255)
                        aaacaaGGGCGACCGTGTTGAAGGTAAAATCAAGTc             <  1:3357886/35‑1 (MQ=255)
                         aacaaGGGCGGCCGTGTTGAAGGTAAAATCAAGTCt            >  1:4584028/1‑36 (MQ=255)
                         aacaaGGGCGGCCGTGTTGAAGGTAAAATCAAGTCt            >  1:380114/1‑36 (MQ=255)
                         aacaaGGGCGACCGTGTTGAAGGTAAAATCAAGTCt            >  1:4160870/1‑36 (MQ=255)
                          acaaGGGCGGCCGTGTTGAAGGTAAAATCAAGTCTa           <  1:7216927/36‑1 (MQ=255)
                            aaGGGCGGCCGTGTTGAAGGTAAAATCAAGTCTATc         >  1:4173429/1‑36 (MQ=255)
                            aaGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATc         >  1:1382847/1‑36 (MQ=255)
                             aGGGCGGCCGTGTTGAAGGTAAAATCAAGTCTATCa        >  1:2019730/1‑36 (MQ=255)
                              gggCGGCCGTGTTGAAGGTAAAATCAAGTCTATCAc       <  1:4613980/36‑1 (MQ=255)
                              gggCGACCGTGTTGAAGGTAAAATCAAGTCTATCAc       >  1:7131077/1‑36 (MQ=255)
                              gggCGACCGTGTTGAAGGTAAAATCAAGTCTATCAc       <  1:4850411/36‑1 (MQ=255)
                               ggCGGCCGTGTTGAAGGTAAAATCAAGTCTATCACt      <  1:4490612/36‑1 (MQ=255)
                                gCGGCCGTGTTGAAGGTAAAATCAAGTCTATCACTg     >  1:7884312/1‑36 (MQ=255)
                                  gACCGTGTTGAAGGGAAAATCAAGTCTATCACTGAc   >  1:4945689/1‑36 (MQ=255)
                                   aCCGTGTTGAAGGTAAAATCAAGTCTATCACTGACt  >  1:904945/1‑36 (MQ=255)
                                   |                                   
GTGGCAGCAGTTCGCGGAAACCCACAACAAGGGCGACCGTGTTGAAGGTAAAATCAAGTCTATCACTGACT  >  REL606/980145‑980215

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG