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% 132.1 / NA 41I21N (ATC→AAC) rpsR30S ribosomal protein S18
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base A (16/23);  minor base C (0/1);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCA  >  REL606/4403982‑4404126
                                                                         |                                                                       
gagaCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACg                                                                       <  2:643169/76‑1 (MQ=255)
 agaCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGa                                                                      <  1:668362/76‑1 (MQ=255)
  gaCTAGCCATATGGTACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGAc                                                                     <  2:811255/76‑1 (MQ=255)
  gaCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGAc                                                                     <  2:14414/76‑1 (MQ=255)
     tAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtat                                                                  >  2:380295/1‑76 (MQ=255)
      aGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACtata                                                                 <  1:1127130/76‑1 (MQ=255)
        ccATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATaaa                                                               <  2:354060/76‑1 (MQ=255)
         cATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAg                                                              <  1:52281/76‑1 (MQ=255)
         cATATGGCACGNTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAg                                                              >  1:166154/1‑76 (MQ=255)
          atatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGa                                                             >  1:455784/1‑76 (MQ=255)
           tatGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGat                                                            >  2:585664/1‑76 (MQ=255)
            atGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGata                                                           >  1:748482/1‑76 (MQ=255)
             tGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGatat                                                          <  1:32691/76‑1 (MQ=255)
               gCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCg                                                        <  1:459543/76‑1 (MQ=255)
                 aCGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCt                                                      <  1:244673/76‑1 (MQ=255)
                  cGTTATTTCCGTCGTCGCAAGTTCTGGCGTTTCNCCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTa                                                     <  1:962529/76‑1 (MQ=255)
                   gTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGACATCGCTAc                                                    >  2:1119669/1‑76 (MQ=255)
                         tCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaa                                              <  1:755212/76‑1 (MQ=255)
                           cgtcgtCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGaaaa                                            <  2:114795/76‑1 (MQ=255)
                              cgtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACt                                         >  1:697120/1‑76 (MQ=255)
                              cgtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACt                                         >  1:453128/1‑76 (MQ=255)
                              cgtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACt                                         >  1:1168420/1‑76 (MQ=255)
                               gtcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTa                                        >  2:808606/1‑76 (MQ=255)
                                tcgCAAGTTCTGCCGTTTCACCGCGGAAGGCGTNCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTAc                                       <  1:968919/76‑1 (MQ=255)
                                  gCAATTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACAt                                     <  2:576814/76‑1 (MQ=255)
                                   cAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAANATATCGCTACGCTGAAAAACTACTTc                                    <  2:1032332/76‑1 (MQ=255)
                                   cAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATc                                    >  2:79844/1‑76 (MQ=255)
                                       ttCTGCCGTTTCATCGCGGAACGCGTTCACGAGACCGACTACAAAGATATCGCTACGCTGAAAAACTGTATCACCg                                <  2:1068897/76‑1 (MQ=255)
                                         cTGCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaa                              <  1:835101/76‑1 (MQ=255)
                                          tGCCGTTTCACCGCGGAAGGCGCTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaaa                             <  1:861183/76‑1 (MQ=255)
                                           gCCGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAg                            >  2:424268/1‑76 (MQ=255)
                                             cGTTTCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg                          <  2:453128/76‑1 (MQ=255)
                                                ttCACCGCGGAAGGCGTTCAAGAGAACGACTATAAAGCTATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTa                       <  1:793186/76‑1 (MQ=255)
                                                     cgcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTNAGAtt                  >  2:959009/1‑76 (MQ=255)
                                                     cgcgGAAGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGAtt                  <  1:380295/76‑1 (MQ=255)
                                                          aaGGCGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTccc             >  2:907988/1‑76 (MQ=255)
                                                              cGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGc         <  1:566666/76‑1 (MQ=255)
                                                              cGTTCAAGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGc         >  1:401297/1‑76 (MQ=255)
                                                               gTTCAAGAGAACGACTATAAAGATATCGCTATGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGcc        >  2:362780/1‑76 (MQ=255)
                                                                   aaGAGAACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTAt    >  2:1144336/1‑76 (MQ=255)
                                                                     gagaACGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCa  <  2:615726/76‑1 (MQ=255)
                                                                         |                                                                       
GAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCGGTAAGATTGTCCCAAGCCGTATCA  >  REL606/4403982‑4404126

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