Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,477,260 G→A 12.1% intergenic (‑45/+236) pyrL ← / ← yjgH pyrBI operon leader peptide/predicted mRNA endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,477,2600GA12.1% 168.3 / 12.2 58intergenic (‑45/+236)pyrL/yjgHpyrBI operon leader peptide/predicted mRNA endoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (4/3);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTA  >  W3110S.gb/4477190‑4477326
                                                                      |                                                                  
aCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg                                                                    <  1:2987456/71‑1 (MQ=255)
aCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg                                                                    <  1:3260861/71‑1 (MQ=255)
aCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg                                                                    <  1:1713715/71‑1 (MQ=255)
 cAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                   >  1:3859729/1‑71 (MQ=255)
 cAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg                                                                    <  1:598693/70‑1 (MQ=255)
  aaaaTGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcg                                                                  >  1:4179631/1‑71 (MQ=255)
  aaaaTGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcg                                                                  >  1:1733894/1‑71 (MQ=255)
     aTGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcggc                                                                >  1:107669/1‑70 (MQ=255)
         cGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAg                                                           >  1:967119/1‑71 (MQ=255)
          gAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAg                                                           >  1:2408684/1‑70 (MQ=255)
          gAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAg                                                           >  1:3382938/1‑70 (MQ=255)
             cacacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAtt                                                       >  1:1063209/1‑71 (MQ=255)
             cacacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAtt                                                       >  1:3484954/1‑71 (MQ=255)
             cacacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAtt                                                       >  1:589269/1‑71 (MQ=255)
              acacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                      <  1:3490389/71‑1 (MQ=255)
              acacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                      <  1:309578/71‑1 (MQ=255)
               cacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTg                                                     <  1:4148331/71‑1 (MQ=255)
                  tGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGc                                                          <  1:2531557/63‑1 (MQ=255)
                       cATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCa                                                             <  1:297773/55‑1 (MQ=255)
                       cATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCa                                                             <  1:2762312/55‑1 (MQ=255)
                         tACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATAtt                                           >  1:4268373/1‑71 (MQ=255)
                           cATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAtt                                                       >  1:1133483/1‑57 (MQ=255)
                           cATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATAtttt                                         >  1:67857/1‑71 (MQ=255)
                           cATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATAtttt                                         >  1:1066989/1‑71 (MQ=255)
                            aTCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                      <  1:947274/57‑1 (MQ=255)
                            aTCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                      <  1:779490/57‑1 (MQ=255)
                            aTCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                      <  1:680971/57‑1 (MQ=255)
                            aTCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTg                                        >  1:2784959/1‑71 (MQ=255)
                             tcctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGt                                       <  1:3066685/71‑1 (MQ=255)
                               ctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCa                                     <  1:2325773/71‑1 (MQ=255)
                                     gCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                                >  1:3543030/1‑70 (MQ=255)
                                       aaaTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                      >  1:1281316/1‑46 (MQ=255)
                                       aaaTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAAc                             >  1:2683637/1‑71 (MQ=255)
                                        aaTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACg                            >  1:143509/1‑71 (MQ=255)
                                         aTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAtt                                                       >  1:796171/1‑43 (MQ=255)
                                         aTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTg                                                     <  1:2618290/45‑1 (MQ=255)
                                         aTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACgg                           <  1:2050006/71‑1 (MQ=255)
                                         aTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACgg                           >  1:3917846/1‑71 (MQ=255)
                                         aTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACgg                           >  1:389089/1‑71 (MQ=255)
                                         aTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACgg                           >  1:471436/1‑71 (MQ=255)
                                         aTTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACgg                           >  1:2419202/1‑71 (MQ=255)
                                          ttGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGc                                                          <  1:3392463/39‑1 (MQ=255)
                                          ttGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAata                                             <  1:4248189/52‑1 (MQ=255)
                                          ttGTCCGGCATTATACTCATCGTAAGAAGCGGCGCAAGCATTTGa                                                    >  1:2738704/1‑45 (MQ=255)
                                             tCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGttt                        <  1:3584028/70‑1 (MQ=255)
                                              ccGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGa                                                    >  1:2225478/1‑41 (MQ=255)
                                               cGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCa                                     <  1:3498248/55‑1 (MQ=255)
                                               cGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAtt                     >  1:3059416/1‑71 (MQ=255)
                                                  cATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTg                   >  1:1643510/1‑70 (MQ=255)
                                                    ttATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATTTTTCGt                                       <  1:1832591/48‑1 (MQ=255)
                                                    ttATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGt                                       <  1:1867300/48‑1 (MQ=255)
                                                    ttATACTCATCGTCAGAAACTGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTa                       <  1:2280541/64‑1 (MQ=255)
                                                    ttATACTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTa                       <  1:1429580/64‑1 (MQ=255)
                                                     tataCTCATCGTCAGAAGCGGCGCAAGCATTTGATGc                                                 <  1:4249947/37‑1 (MQ=255)
                                                     tataCTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGt                                       <  1:1080976/47‑1 (MQ=255)
                                                      ataCTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAgcgc                                 <  1:617724/52‑1 (MQ=255)
                                                       taCTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaa             <  1:1514488/71‑1 (MQ=255)
                                                          tCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAg          >  1:2024250/1‑71 (MQ=255)
                                                             tCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAAc                             <  1:3884665/49‑1 (MQ=255)
                                                             tCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTc        <  1:2479052/70‑1 (MQ=255)
                                                             tCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa       <  1:3186584/71‑1 (MQ=255)
                                                             tCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                                <  1:4011402/46‑1 (MQ=255)
                                                               gTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAgc                                   >  1:2250016/1‑41 (MQ=255)
                                                               gTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGt         >  1:4133516/1‑67 (MQ=255)
                                                               gTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGt         >  1:3551531/1‑67 (MQ=255)
                                                                  aGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGta  >  1:1300226/1‑71 (MQ=255)
                                                                      |                                                                  
ACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTA  >  W3110S.gb/4477190‑4477326

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