Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,471,518 T→A I325N (ATC→AAC)  yjgQ → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,471,5180TA98.4% 227.6 / NA 61I325N (ATC→AAC) yjgQconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (44/16);  minor base C (0/1);  total (44/17)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTT  >  REL606/4471447‑4471590
                                                                       |                                                                        
tagcgtaCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTc                                                                      >  1:2112048/1‑76 (MQ=255)
 agcgtaCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCg                                                                     <  1:282006/76‑1 (MQ=255)
   cgtaCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGc                                                                   >  1:318235/1‑76 (MQ=255)
    gtaCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGcc                                                                  >  2:1894739/1‑76 (MQ=255)
     taCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTTTCTTCTACGTACTGGACCAGAACTTCCGccc                                                                 >  1:617325/1‑76 (MQ=255)
     taCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAAATTCGGccc                                                                 >  1:1974562/1‑76 (MQ=255)
      aCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGcccc                                                                >  1:1381032/1‑75 (MQ=255)
      aCCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCg                                                                <  1:593291/76‑1 (MQ=255)
        cGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCt                                                              >  1:3960/1‑76 (MQ=255)
          aTGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGa                                                            >  1:1962390/1‑76 (MQ=255)
           tGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCCCTGAc                                                           >  2:1118043/1‑76 (MQ=255)
              gcgtgcgtGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGtt                                                        >  1:2136934/1‑76 (MQ=255)
                  gcgtgTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGtt                                                    >  2:995882/1‑76 (MQ=255)
                  gcgtgTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGtt                                                    <  2:536516/76‑1 (MQ=255)
                    gtgtGGTCACCGGTATCAGTTTCGGTTTTTTCTTCTACGTACTGGACCAGAACTTCGGCCCGCCGCCGTTGGTTTa                                                  >  1:987330/1‑76 (MQ=255)
                    gtgtGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTa                                                  <  2:1969631/76‑1 (MQ=255)
                    gtgtGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTa                                                  >  2:1505171/1‑76 (MQ=255)
                     tgtgGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTAt                                                 >  1:100401/1‑76 (MQ=255)
                      gtgGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCCGACGTTGGTTTTTg                                                >  1:2169349/1‑76 (MQ=255)
                        ggTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGc                                              >  1:2106409/1‑76 (MQ=255)
                         gacACCGGTATCAATTTCGGTGTTCTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCa                                             <  1:1283817/74‑1 (MQ=255)
                         gTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCCGACGTTGGTTTATGGCa                                             >  1:207970/1‑76 (MQ=255)
                           cACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATc                                           <  2:775496/76‑1 (MQ=255)
                            aCCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATcc                                          >  2:501306/1‑76 (MQ=255)
                            aCCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATcc                                          <  2:1803225/76‑1 (MQ=255)
                              cgttATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCccg                                        <  2:231274/73‑1 (MQ=255)
                              cGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTTTGGCATcccc                                        >  1:1068682/1‑75 (MQ=255)
                              cGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCccg                                        >  1:64127/1‑76 (MQ=255)
                                 tATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCccgccg                                     >  1:644644/1‑76 (MQ=255)
                                 tATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCCGACGTTGGTTTATGGCATCcccccg                                     >  1:1891981/1‑76 (MQ=255)
                                     aGTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGatca                                 <  2:1954875/76‑1 (MQ=255)
                                      gTTTCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGatcat                                >  2:164185/1‑76 (MQ=255)
                                       tttCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGgtcatc                               >  1:114941/1‑76 (MQ=255)
                                        ttCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGGTCATCg                              >  1:532345/1‑76 (MQ=255)
                                        ttCGGTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGCTCATCg                              >  1:957336/1‑76 (MQ=255)
                                           ggTTTTTTCTTCTCCGTACTGGACCAGAACTTCGGCCCGCTGTCGTTGGTTTATGGCATCCCGCCGATCATCGgcg                           <  1:404898/76‑1 (MQ=255)
                                           ggTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGgcg                           >  2:364738/1‑76 (MQ=255)
                                           ggTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGgcg                           >  1:566825/1‑76 (MQ=255)
                                           ggTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGgcg                           >  1:1736002/1‑76 (MQ=255)
                                           ggTTTTGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGgcg                           >  1:2059855/1‑76 (MQ=255)
                                             ttttGTCTTCTACGTACTGGACCAGACCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCa                         <  2:1340342/76‑1 (MQ=255)
                                             ttttGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGgcgcc                         >  1:279218/1‑75 (MQ=255)
                                             ttttGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCa                         >  1:1839336/1‑76 (MQ=255)
                                               ttGTCTTCTACCTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACt                       >  1:781991/1‑76 (MQ=255)
                                                tGTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTg                      >  1:1023732/1‑76 (MQ=255)
                                                 gTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGt                     >  1:1823432/1‑76 (MQ=255)
                                                 gTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGt                     >  1:2112853/1‑76 (MQ=255)
                                                 gTCTTCTACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGt                     >  1:2153001/1‑76 (MQ=255)
                                                  tcttctACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGtt                    <  1:1802507/76‑1 (MQ=255)
                                                     tctACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGcc                 >  1:10469/1‑76 (MQ=255)
                                                     tctACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGcc                 >  2:870503/1‑76 (MQ=255)
                                                      ctACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCa                <  2:200484/76‑1 (MQ=255)
                                                       tACGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCaa               >  1:1957217/1‑76 (MQ=255)
                                                        aCGTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAg              >  1:1785512/1‑76 (MQ=255)
                                                          gTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAgcg            >  2:1930576/1‑76 (MQ=255)
                                                          gTACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAgcg            <  1:833128/76‑1 (MQ=255)
                                                           tACTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCCAgcgc           >  1:1215886/1‑76 (MQ=255)
                                                            aCTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGcc          <  2:1035709/76‑1 (MQ=255)
                                                             cTGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCa         <  1:2119082/76‑1 (MQ=255)
                                                              tGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAg        >  1:1294491/1‑76 (MQ=255)
                                                              tGGACCAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCCGCGCACTGTTGCCAAGCGCCAg        >  1:865043/1‑76 (MQ=255)
                                                                  ccAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttc    <  2:2139497/76‑1 (MQ=255)
                                                                  ccAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttc    <  2:69525/76‑1 (MQ=255)
                                                                  ccAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCCGCGCACTGTTGCCCAGCGCCAGcttc    >  2:931721/1‑76 (MQ=255)
                                                                   cAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttct   >  2:1189524/1‑76 (MQ=255)
                                                                   cAGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttct   <  1:648533/76‑1 (MQ=255)
                                                                    aGAACTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttctt  >  1:61378/1‑76 (MQ=255)
                                                                    aGAACTTCGGCCCGCTGACGTTAGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGcttctt  >  1:2219852/1‑76 (MQ=255)
                                                                       |                                                                        
TAGCGTACCGATGGGCGTGCGTGTGGTCACCGGTATCAGTTTCGGTTTTGTCTTCTACGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTT  >  REL606/4471447‑4471590

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