Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,404,055 T→A I21N (ATC→AAC)  rpsR → 30S ribosomal protein S18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,404,0550TA97.5% 150.2 / ‑5.8 41I21N (ATC→AAC) rpsR30S ribosomal protein S18
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base A (17/22);  minor base G (1/0);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATC  >  REL606/4403983‑4404125
                                                                        |                                                                      
agaCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGa                                                                     >  1:784390/1‑76 (MQ=255)
  aCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACt                                                                   <  1:663140/76‑1 (MQ=255)
    tAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtat                                                                 >  1:593450/1‑76 (MQ=255)
     aGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtata                                                                >  2:669954/1‑76 (MQ=255)
      gCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCNGTGTCACCGCGGAAGGCGTTCAAGAGAACGACTATaa                                                               <  1:1233256/76‑1 (MQ=255)
        cATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAg                                                             <  1:1200489/76‑1 (MQ=255)
        cATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAg                                                             <  2:358433/76‑1 (MQ=255)
              gCACGTTATTTCCGTCGTCCCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCg                                                       >  1:884443/1‑76 (MQ=255)
                      tttCCTTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTAGGCTg                                               <  2:149191/76‑1 (MQ=255)
                      tttCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTg                                               >  1:972085/1‑76 (MQ=255)
                          cgtcgtCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaaaa                                           <  1:120086/76‑1 (MQ=255)
                              gtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTa                                       >  2:525025/1‑76 (MQ=255)
                              gtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTa                                       >  2:37831/1‑76 (MQ=255)
                                 gCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACAt                                    <  1:239950/76‑1 (MQ=255)
                                  cAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATc                                   <  1:745018/76‑1 (MQ=255)
                                  cAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATc                                   <  1:673715/76‑1 (MQ=255)
                                   aaGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCa                                  <  1:659149/76‑1 (MQ=255)
                                    aGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAGAACTACATCAc                                 <  1:35035/76‑1 (MQ=255)
                                       tCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGa                              >  1:969668/1‑76 (MQ=255)
                                         tGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaaa                            <  1:364173/76‑1 (MQ=255)
                                          gCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATAGCGCTACGCTGAAAAACTACATCACCGAAAg                           >  1:27210/1‑76 (MQ=255)
                                            cGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg                         >  1:357274/1‑76 (MQ=255)
                                              tttCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGt                       <  1:543935/76‑1 (MQ=255)
                                               ttCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCNGTa                      <  1:1142204/76‑1 (MQ=255)
                                               ttCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTa                      >  1:592268/1‑76 (MQ=255)
                                               ttCACCGCGGAAGGCGTTAAAGGGAACNACNACNAACATATNGCGACGATGAAAAACTACATCACCGAAAGCGGTa                      <  1:1101773/76‑1 (MQ=255)
                                                tCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTaa                     >  1:115345/1‑76 (MQ=255)
                                                tCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTaa                     >  2:31728/1‑76 (MQ=255)
                                                 cACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAg                    <  1:261654/76‑1 (MQ=255)
                                                 cACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAg                    >  2:419224/1‑76 (MQ=255)
                                                 cACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAg                    >  1:388792/1‑76 (MQ=255)
                                                    cgcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGAtt                 <  2:507097/76‑1 (MQ=255)
                                                     gcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTg                <  1:699241/76‑1 (MQ=255)
                                                     gcgGAAGGCGTCCCAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGAttc                <  2:16150/76‑2 (MQ=255)
                                                      cgGAAGGCGTTCAAGAGAGCGNCTNTAAAGATTTCGCTACGCTGAANAACTACANCNCCGAAAGCGGTAAGATTGt               >  2:182360/1‑76 (MQ=255)
                                                        gAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTcc             >  2:489048/1‑76 (MQ=255)
                                                             cGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGc        <  1:837808/76‑1 (MQ=255)
                                                               ttCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCg      <  1:231624/76‑1 (MQ=255)
                                                               ttCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCg      <  2:511197/76‑1 (MQ=255)
                                                                tCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGt     >  1:833141/1‑76 (MQ=255)
                                                                  aaGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTAt   <  1:514330/76‑1 (MQ=255)
                                                                  aaGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTAt   >  1:525475/1‑76 (MQ=255)
                                                                   agagAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATc  <  2:520111/76‑1 (MQ=255)
                                                                   agagAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATc  <  2:563630/76‑1 (MQ=255)
                                                                        |                                                                      
AGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATC  >  REL606/4403983‑4404125

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: