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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,404,0300TA98.6% 273.8 / ‑3.1 71F13I (TTC→ATC) rpsR30S ribosomal protein S18
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (36/34);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

GAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCG  >  REL606/4403957‑4404102
                                                                         |                                                                        
gAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATc                                                                        <  1:39883/76‑1 (MQ=255)
gAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATc                                                                        >  1:603880/1‑76 (MQ=255)
 aGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCa                                                                       >  2:473238/1‑76 (MQ=255)
  gCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCAc                                                                      <  1:491897/76‑1 (MQ=255)
   cAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCAcc                                                                     >  2:581200/1‑76 (MQ=255)
   cAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCAcc                                                                     >  1:734763/1‑76 (MQ=255)
    aGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACcg                                                                    >  1:1231420/1‑76 (MQ=255)
     gATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACcgc                                                                   <  2:897981/76‑1 (MQ=255)
       ttgaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCgg                                                                 >  2:1118333/1‑76 (MQ=255)
       ttgaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCgg                                                                 >  1:1481834/1‑76 (MQ=255)
       ttgaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCgg                                                                 >  1:110071/1‑76 (MQ=255)
       ttgaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCgg                                                                 >  1:992651/1‑76 (MQ=255)
       ttgaattgaTAGATGCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCgg                                                                 <  1:138196/76‑1 (MQ=255)
        tgaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGa                                                                >  1:753380/1‑76 (MQ=255)
        tgaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGa                                                                >  1:373358/1‑76 (MQ=255)
        agaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGa                                                                >  1:824880/2‑76 (MQ=255)
         gaattgaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGaa                                                               >  1:672366/1‑76 (MQ=255)
              gaTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGt                                                          <  1:93635/76‑1 (MQ=255)
               aTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGtt                                                         >  2:1443480/1‑76 (MQ=255)
               aTAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGtt                                                         <  1:390903/76‑1 (MQ=255)
                  gATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCaa                                                      <  1:56016/76‑1 (MQ=255)
                  gATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCaa                                                      >  2:721565/1‑76 (MQ=255)
                  gATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCaa                                                      <  2:94617/76‑1 (MQ=255)
                   aTTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAag                                                     <  2:1069825/76‑1 (MQ=255)
                   aTTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAag                                                     >  1:948311/1‑76 (MQ=255)
                   aTTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAag                                                     >  1:258196/1‑76 (MQ=255)
                   aTTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAag                                                     <  1:1282521/76‑1 (MQ=255)
                   aTTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAag                                                     <  2:894258/76‑1 (MQ=255)
                      cTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAgaga                                                  <  2:1208829/76‑1 (MQ=255)
                         gagaCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCg                                               <  1:284785/76‑1 (MQ=255)
                          agaCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGa                                              <  1:369745/76‑1 (MQ=255)
                            aCTAGCCATATGGCACGTTATTTCCGTCGTTGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACt                                            >  1:912656/1‑76 (MQ=255)
                            aCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACt                                            <  2:436052/76‑1 (MQ=255)
                            aCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACt                                            <  1:608684/76‑1 (MQ=255)
                            aCTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACt                                            <  2:1132839/76‑1 (MQ=255)
                             cTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACta                                           >  2:330867/1‑76 (MQ=255)
                             cTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACta                                           <  1:1097300/76‑1 (MQ=255)
                              tAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACtat                                          >  2:785336/1‑76 (MQ=255)
                              tAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACtat                                          >  2:114569/1‑76 (MQ=255)
                               aGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACtata                                         >  1:682833/1‑76 (MQ=255)
                               aGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACtata                                         >  1:1237564/1‑76 (MQ=255)
                                gCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATaa                                        <  2:1427227/76‑1 (MQ=255)
                                 ccATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATaaa                                       <  2:1051327/76‑1 (MQ=255)
                                  cATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAg                                      >  1:1445644/1‑76 (MQ=255)
                                  cATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAg                                      >  1:340728/1‑76 (MQ=255)
                                     atGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGata                                   <  1:1082029/76‑1 (MQ=255)
                                      tGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGttat                                  <  2:982411/76‑1 (MQ=255)
                                        gCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCg                                >  1:90491/1‑76 (MQ=255)
                                          aCGTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCt                              <  1:623090/76‑1 (MQ=255)
                                           ctttATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTa                             >  1:1083278/3‑76 (MQ=255)
                                            gTTATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTAc                            <  2:505587/76‑1 (MQ=255)
                                             ttATTTCCGTCGTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACg                           <  2:570778/76‑1 (MQ=255)
                                                   ccGTCTTCGCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGaaa                     <  2:887768/76‑1 (MQ=255)
                                                   ccGTCGTCGCAAGTTCTGCCGTATCTCCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGaaa                     <  1:1335067/76‑1 (MQ=255)
                                                       cgtcgCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACt                 >  1:212445/1‑76 (MQ=255)
                                                       cgtcgCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACt                 <  1:715918/76‑1 (MQ=255)
                                                          cgCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACa              <  2:1095714/76‑1 (MQ=255)
                                                           gCAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACAt             <  1:393068/76‑1 (MQ=255)
                                                           gCAAGTTCTGCCGATTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACAt             >  1:997055/1‑76 (MQ=255)
                                                            cAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATc            <  2:569261/76‑1 (MQ=255)
                                                            cAAGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATc            >  1:356738/1‑76 (MQ=255)
                                                             aaGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCa           <  1:128059/76‑1 (MQ=255)
                                                              aGTTCTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCAc          >  1:566970/1‑76 (MQ=255)
                                                                  cTGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaa      >  1:917148/1‑76 (MQ=255)
                                                                   tGCTGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaaa     <  2:213727/76‑1 (MQ=255)
                                                                   tGCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGaaa     >  2:956271/1‑76 (MQ=255)
                                                                    gCCGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAg    <  2:1481326/76‑1 (MQ=255)
                                                                     ccGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGc   >  2:416713/1‑76 (MQ=255)
                                                                      cGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg  >  2:1189482/1‑76 (MQ=255)
                                                                      cGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg  <  1:742390/76‑1 (MQ=255)
                                                                      cGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg  <  2:74372/76‑1 (MQ=255)
                                                                      cGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg  >  1:916965/1‑76 (MQ=255)
                                                                      cGTATCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCg  >  2:1211091/1‑76 (MQ=255)
                                                                         |                                                                        
GAGCAGATTGAATTGATAGATTCTGGAGACTAGCCATATGGCACGTTATTTCCGTCGTCGCAAGTTCTGCCGTTTCACCGCGGAAGGCGTTCAAGAGATCGACTATAAAGATATCGCTACGCTGAAAAACTACATCACCGAAAGCG  >  REL606/4403957‑4404102

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