Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,446,544 T→A 16.1% L689L (CTT→CTA ydbD → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,446,5440T→A16.1% 6.1 ‑3.7 0.0342067ydbDhypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (15/13):  total (15/19)
Fisher's exact test strand distribution p-value = 2.39e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGC  >  REL606/1446509‑1446576
                                   |                                
aaaaGCAATCTCACATTTTTCAAGTTCGGCGCACtt                                  <  1:7340087/36‑1 (MQ=255)
 aaaGCAATCTCACATTTTTCAAGTTCGGCGCACTTa                                 >  1:22991/1‑36 (MQ=255)
 aaaGCAATCTCACATTTTTCAAGTTCGGCGCACTTa                                 <  1:2804747/36‑1 (MQ=255)
 aaaGCAATCTCACATTTTTCAAGTTCGGCGCACTTa                                 <  1:2978276/36‑1 (MQ=255)
  aaGCAATCTCACATTTTTCAAGTTCGGCGCACTTaa                                <  1:1234431/36‑1 (MQ=255)
  aaGCAATCTCACATTTTTCAAGTTCGGCGCACTTaa                                <  1:858327/36‑1 (MQ=255)
  aaGCAATCTCACATTTTTCAAGTTCGGCGCACTTaa                                >  1:7275608/1‑36 (MQ=255)
  aaGCAATCTCACATTTTTCAAGTTCGGCGCACTTaa                                <  1:4165785/36‑1 (MQ=255)
   aGCAATCTCACATTTTTCAAGTTCGGCGCACTTaac                               <  1:1849563/36‑1 (MQ=255)
    gCAATCTCACATTTTTCAAGTTCGGCGCACTTaaca                              <  1:7142951/36‑1 (MQ=255)
      aaTCTCACATTTTTCAAGTTCGGCGCACTTaacaac                            <  1:1530683/36‑1 (MQ=255)
      aaTCTCACATTTTTCAAGTTCGGCGCACTTaacaac                            <  1:4235727/36‑1 (MQ=255)
       aTCTCACATTTTTCAAGTTCGGCGCACTTAACAAcc                           >  1:5712676/1‑36 (MQ=255)
       aTCTCACATTTTTCAAGTTCGGCGCACTTAACAAcc                           <  1:5559988/36‑1 (MQ=255)
       aTCTCACATTTTTCAAGTTCGGCGCACTTAACAAcc                           <  1:5468891/36‑1 (MQ=255)
        tctcACATTTTTCAAGTTCGGCGCACTTAACAACCg                          <  1:858387/36‑1 (MQ=255)
        tctcACATTTTTCAAGTTCGGCGCACTTAACAACCg                          <  1:2750473/36‑1 (MQ=255)
          tcACATTTTTCAAGTTCGGCGCACTTAACAACCGcc                        >  1:5912729/1‑36 (MQ=255)
          tcACATTTTTCAAGTTCGGCGCACTTAACAACCGcc                        >  1:7306546/1‑36 (MQ=255)
               tttttCAAGTTCGGCGCTCTTAACAACCGCcacaca                   >  1:1494410/1‑36 (MQ=255)
               ttttaCAAGTCCGGCCCATAAAAAAACCGCcacaca                   <  1:4843199/36‑1 (MQ=0)
                ttttCAAGTTCGGCGCACTTAACAACCGCCacacac                  >  1:7938408/1‑36 (MQ=255)
                 tttCAAGTTCGGCGCACTTAACAACCGACACACAcc                 >  1:2123531/1‑36 (MQ=255)
                  ttCAAGTTCGGCGCACTTAACAACCGCCACACACCg                >  1:1401799/1‑36 (MQ=255)
                  ttCAAGTTCGGCGCACTTAACAACCGCCACACACCg                >  1:2506481/1‑36 (MQ=255)
                   caaaGTCCGCCCCACTAAAAAACCGCCACACACCGc               <  1:4734222/34‑1 (MQ=12)
                     aaGTTCGGCGCACTTAACAACCGCCACACACCGCCt             >  1:2722676/1‑36 (MQ=255)
                     aaGTTCGGCCCACTAAAAAACCGCCACACACCGCCt             <  1:682853/36‑1 (MQ=37)
                      aGTTCGGCGCACTTAACAACCGCCACACACCGCCTg            >  1:4545787/1‑36 (MQ=255)
                           ggCGCACTTAACAACCGCCACACACCGCCTGCAGTa       >  1:3989293/1‑36 (MQ=255)
                           gccgcACTAAAAAACCGCCACACACCGCCTGCAGTa       <  1:7541386/34‑1 (MQ=37)
                             ggcACTAAAACACCGCAACACACCGCCTGCAGTACt     <  1:6977806/35‑1 (MQ=16)
                               cACTTAACAACCGCCACACACCGCCTGCAGTACTgg   >  1:67217/1‑36 (MQ=255)
                                aCTTAACAACCGCCACACACCGCCTGCAGTACTggc  >  1:4658035/1‑36 (MQ=255)
                                aCTAAAAAACCGCCACACACCGCCTGCAGTACTggc  <  1:144748/36‑1 (MQ=37)
                                   |                                
AAAAGCAATCTCACATTTTTCAAGTTCGGCGCACTTAACAACCGCCACACACCGCCTGCAGTACTGGC  >  REL606/1446509‑1446576

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