Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 4,421,418 C→T R12Q (CGG→CAG)  yjfO ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,421,4180CT75.8% 233.9 / 81.5 124R12Q (CGG→CAG) yjfOconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (57/37);  total (72/52)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 124 total aligned reads displayed.

GCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTA  >  W3110S.gb/4421353‑4421466
                                                                 |                                                
gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGt                                               >  1:600870/1‑69 (MQ=255)
gctgTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTa                                              >  1:1768417/1‑70 (MQ=255)
gctgTGGCGTACCTTGCAGTGCGCTACAGGCAATTAACATCAACACCAATAATAAACTGGCAAACCGGTa                                              >  1:2838372/1‑70 (MQ=255)
    tgGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAAt                                           <  1:4107466/69‑1 (MQ=255)
      gCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAAc                                        >  1:2408376/1‑70 (MQ=255)
      gCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACg                                       >  1:3163305/1‑71 (MQ=255)
      gCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACg                                       >  1:2663826/1‑71 (MQ=255)
                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:3499162/64‑1 (MQ=255)
                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:3950957/64‑1 (MQ=255)
                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:1898894/64‑1 (MQ=255)
                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:1229226/64‑1 (MQ=255)
                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:2763931/64‑1 (MQ=255)
                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:3574999/64‑1 (MQ=255)
                 aGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:1000537/64‑1 (MQ=255)
                   tGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTa                                 >  1:2229599/1‑64 (MQ=255)
                   tGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTa                                 >  1:3891598/1‑64 (MQ=255)
                   tGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTa                                 >  1:83289/1‑64 (MQ=255)
                     ggcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTg                        >  1:1624113/2‑71 (MQ=255)
                     cgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCt                         >  1:3824555/1‑70 (MQ=255)
                     cgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTg                        >  1:3504233/1‑71 (MQ=255)
                     cgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTg                        >  1:3672001/1‑71 (MQ=255)
                     cgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTg                        >  1:556676/1‑71 (MQ=255)
                     cgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTg                        >  1:3376388/1‑71 (MQ=255)
                     cgcTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTg                        >  1:12992/1‑71 (MQ=255)
                        tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGttt                            <  1:1635595/64‑1 (MQ=255)
                        tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGttt                            <  1:1647211/64‑1 (MQ=255)
                        tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGttt                            <  1:1385545/64‑1 (MQ=255)
                        tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGttt                            <  1:2258200/64‑1 (MQ=255)
                        tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGttt                            <  1:63352/64‑1 (MQ=255)
                        tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGttt                            <  1:1446549/64‑1 (MQ=255)
                        tACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGttt                            <  1:2096956/64‑1 (MQ=255)
                               aCTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   <  1:3985197/50‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCt                                     <  1:3846070/47‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCt                                     <  1:454199/47‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCt                                     <  1:1338751/47‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCt                                     <  1:1361014/47‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:3464706/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:3042185/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:2310645/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:1842646/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:133733/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:1276284/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:60034/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:1000790/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTACCAtagtg            >  1:2762842/1‑71 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGc                                      <  1:3644007/46‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGc                                      <  1:2617361/46‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGc                                      <  1:466959/46‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGc                                      <  1:3293693/46‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGc                                      <  1:1112184/46‑1 (MQ=255)
                                cTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGc                                      <  1:286965/46‑1 (MQ=255)
                                cTTAACATCAACACAAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:976836/1‑71 (MQ=255)
                                cTTAACAACAACACCAATAATAAACTGGCAAACTGCTAAATAACGCTATTACGTTTCCTGCTAACCAtagt             >  1:2417281/1‑71 (MQ=255)
                                  tAACATCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTAt                                   >  1:1161044/1‑47 (MQ=255)
                                     cATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGt                              <  1:3119142/49‑1 (MQ=255)
                                     cATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGt                              <  1:3766693/49‑1 (MQ=255)
                                     cATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGt                              <  1:3221478/49‑1 (MQ=255)
                                     cATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGt                              <  1:1226486/49‑1 (MQ=255)
                                     cATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGt                              <  1:357758/49‑1 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCa                 <  1:2949667/60‑1 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCa                 <  1:299547/60‑1 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCa                 <  1:1760065/60‑1 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCa                 <  1:1213459/60‑1 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCa                 <  1:1397065/60‑1 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCa                 <  1:1267986/60‑1 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGttt                            >  1:1924853/1‑49 (MQ=255)
                                       tCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGttt                            >  1:356613/1‑49 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:133682/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:543794/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:2362042/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:2112767/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:1632822/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:1162680/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:3681240/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:868384/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:3160178/1‑60 (MQ=255)
                                         aaCACCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:3851941/1‑60 (MQ=255)
                                           cacCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta  >  1:1284904/1‑71 (MQ=255)
                                           cacCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta  >  1:3771018/1‑71 (MQ=255)
                                           cacCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta  >  1:3811611/1‑71 (MQ=255)
                                           cacCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta  >  1:2998958/1‑71 (MQ=255)
                                           cacCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCtta  >  1:3650070/1‑71 (MQ=255)
                                            acCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTc                           >  1:353110/1‑45 (MQ=255)
                                            acCAATAATAAACTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCa                 <  1:2898391/55‑1 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:3857994/1‑51 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:2047461/1‑51 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:775609/1‑51 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:2697767/1‑51 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:2686810/1‑51 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:3577742/1‑51 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:1410119/1‑51 (MQ=255)
                                                  aataaACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCAta               >  1:2079503/1‑51 (MQ=255)
                                                      aaCTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc      <  1:501577/56‑1 (MQ=255)
                                                      aaCTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc      <  1:1225255/56‑1 (MQ=255)
                                                      aaCTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc      <  1:4083199/56‑1 (MQ=255)
                                                      aaCTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc      <  1:4015970/56‑1 (MQ=255)
                                                       aCTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGc      <  1:351471/55‑1 (MQ=255)
                                                        cTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTaa                    >  1:229230/1‑40 (MQ=255)
                                                        cTGGCAAACTGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCt    >  1:2234299/1‑56 (MQ=255)
                                                        cTGCCAAACTGGTAAATAACGCTATTACGTTTcctcct                      >  1:3270132/1‑38 (MQ=37)
                                                                 |                                                
GCTGTGGCGTACCTTGCAGTGCGCTACAGGCACTTAACATCAACACCAATAATAAACTGGCAAACCGGTAAATAACGCTATTACGTTTCCTGCTAACCATAGTGTAGTGCCTTA  >  W3110S.gb/4421353‑4421466

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