Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,800,287 C→T 5.7% intergenic (‑85/‑137) yicS ← / → nlpA hypothetical protein/cytoplasmic membrane lipoprotein‑28

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,800,2870CT 5.7% 270.6 / 3.7 87intergenic (‑85/‑137)yicS/nlpAhypothetical protein/cytoplasmic membrane lipoprotein‑28
Reads supporting (aligned to +/- strand):  ref base C (43/39);  new base T (3/2);  total (46/41)
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.78e-01

ATGTGTATCGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCTTATCGATAG  >  W3110S.gb/3800221‑3800354
                                                                  |                                                                   
aTGTGTATCGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGc                                                                 >  1:316110/1‑71 (MQ=255)
aTGTGTATCGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGc                                                                 >  1:409363/1‑71 (MQ=255)
   tgtATCGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGt                                                               <  1:2870075/70‑1 (MQ=255)
       tCGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCg                                                                <  1:522366/65‑1 (MQ=255)
       tCGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGt                                                            >  1:1049503/1‑69 (MQ=255)
        cGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAg                                                         >  1:3105488/1‑71 (MQ=255)
         gAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAg                                                         >  1:2868012/1‑70 (MQ=255)
         gAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGt                                                        <  1:1977201/71‑1 (MQ=255)
         gAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGt                                                        <  1:3850805/71‑1 (MQ=255)
          aaTGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAg                                                         >  1:1218865/1‑69 (MQ=255)
          aaTGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGt                                                        >  1:2778290/1‑70 (MQ=255)
          aaTGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGt                                                        >  1:1200446/1‑70 (MQ=255)
          aaTGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTa                                                       <  1:725962/71‑1 (MQ=255)
          aaTGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACAAAACGCGCTTGCGTGGTCAGt                                                        >  1:1429298/1‑70 (MQ=255)
           aTGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaa                                                      <  1:2353945/71‑1 (MQ=255)
            tGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGt                                                        >  1:2589740/1‑68 (MQ=255)
            tGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaag                                                     <  1:1621941/71‑1 (MQ=255)
                 cATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaag                                                     <  1:665239/66‑1 (MQ=255)
                 cATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTg                                                <  1:118867/71‑1 (MQ=255)
                   tAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACgcgc                                                                     <  1:10324/48‑1 (MQ=255)
                   tAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACgcgc                                                                     <  1:598022/48‑1 (MQ=255)
                   tAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaag                                                     <  1:2290572/64‑1 (MQ=255)
                   tAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGc                                               >  1:1297763/1‑70 (MQ=255)
                   tAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGc                                               >  1:3357174/1‑70 (MQ=255)
                     aGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGt                                                            >  1:2404199/1‑55 (MQ=255)
                     aGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCa                                             <  1:3143648/70‑1 (MQ=255)
                     aGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCa                                             <  1:1557789/70‑1 (MQ=255)
                      gtgtAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaa                                                      <  1:842328/60‑1 (MQ=255)
                       tgtAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTTGTCAGTaa                                                      <  1:2146663/59‑1 (MQ=255)
                         tAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTg                                                              <  1:1687767/49‑1 (MQ=255)
                         tAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTg                                                              <  1:2481608/49‑1 (MQ=255)
                         tAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTg                                                              <  1:1386442/49‑1 (MQ=255)
                          aGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACtt                                        >  1:239928/1‑70 (MQ=255)
                           gTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTta                                      >  1:1246369/1‑71 (MQ=255)
                           gTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTta                                      <  1:2046156/71‑1 (MQ=255)
                              tGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTatat                                   >  1:2962208/1‑71 (MQ=255)
                               gCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTatat                                   <  1:2540975/70‑1 (MQ=255)
                               gCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATAtt                                  >  1:962416/1‑71 (MQ=255)
                                cTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTta                                      <  1:552566/66‑1 (MQ=255)
                                 tttCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaaga                                                    <  1:3819484/51‑1 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCg                                                                <  1:2788243/37‑1 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaagaa                                                   >  1:2145469/1‑50 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTaaga                                                    >  1:2465910/1‑49 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:2979543/1‑65 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:2020314/1‑65 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTtat                                     >  1:3041539/1‑64 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTcc                                <  1:3823316/69‑1 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCAc                              >  1:3171028/1‑71 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCAc                              >  1:1005199/1‑71 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCAc                              >  1:2986797/1‑71 (MQ=255)
                                   tCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCAc                              >  1:211293/1‑71 (MQ=255)
                                     aTCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACtt                            >  1:235164/1‑71 (MQ=255)
                                     aTCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACtt                            <  1:3516929/71‑1 (MQ=255)
                                     aTCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACt                             >  1:1164403/1‑70 (MQ=255)
                                       cGCTCTCATCATGCTGTACGAAACGCGCGTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACttt                           >  1:2516935/1‑70 (MQ=255)
                                           ctcATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTc                      >  1:2859766/1‑71 (MQ=255)
                                           ctcATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTc                      >  1:2450010/1‑71 (MQ=255)
                                            tcatcaTGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACttt                           <  1:3501183/65‑1 (MQ=255)
                                             catcatGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCt                    >  1:851024/1‑71 (MQ=255)
                                              atcatGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCt                    <  1:2811989/70‑1 (MQ=255)
                                               tcatGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCt                    <  1:1193403/69‑1 (MQ=255)
                                               tcatGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCt                    <  1:794216/69‑1 (MQ=255)
                                               tcatGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCt                    <  1:580436/69‑1 (MQ=255)
                                               tcatGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCt                    <  1:3294931/69‑1 (MQ=255)
                                                catGCTGTACGAAAGGTGTTTGCGTGGTCAGTaagaa                                                   >  1:1115128/1‑37 (MQ=25)
                                                catGCTGTACGAAAGGTGTTTGCGTGGTCAGTaagaa                                                   >  1:2860657/1‑37 (MQ=25)
                                                catGCTGTACGAAAGGTGTTTGCGTGGTCAGTaagaa                                                   >  1:3337020/1‑37 (MQ=25)
                                                 atGCTGTACGAAACGCGCTTGCGTGGTCAGTaagaa                                                   >  1:1607563/1‑36 (MQ=255)
                                                 atGCTGTACGAAACGCGCTTGCGTGGTCAGTaagaa                                                   >  1:2349598/1‑36 (MQ=255)
                                                 atGCTGTACGAAACGCGCTTGCGTGGTCAGTaagaa                                                   >  1:2888612/1‑36 (MQ=255)
                                                 atGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCtttt                 <  1:16189/70‑1 (MQ=255)
                                                 atGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCtttt                 >  1:1663205/1‑70 (MQ=255)
                                                    cTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTatat                                   <  1:2970633/49‑1 (MQ=255)
                                                    cTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTAt              <  1:1290852/70‑1 (MQ=255)
                                                      gTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCt           >  1:3440533/1‑71 (MQ=255)
                                                       tACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCtt          <  1:3353033/71‑1 (MQ=255)
                                                       tACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCtt          <  1:2051529/71‑1 (MQ=255)
                                                       tACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCt           >  1:318418/1‑70 (MQ=255)
                                                        aCGAAAGGTGTTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCtt          <  1:3110260/70‑1 (MQ=255)
                                                        aCGAAAGGTGTTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCtt          <  1:3009999/70‑1 (MQ=255)
                                                        aCGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACtt                                        <  1:2257174/40‑1 (MQ=255)
                                                        aCGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:569279/1‑44 (MQ=255)
                                                         cGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCa                               >  1:1026098/1‑48 (MQ=255)
                                                         cGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTa                         >  1:613871/1‑54 (MQ=255)
                                                            aaCGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTta                                      <  1:478146/38‑1 (MQ=255)
                                                            aaCGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTta                                      <  1:1235939/38‑1 (MQ=255)
                                                            aaCGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCa                               >  1:1913572/1‑45 (MQ=255)
                                                            aaCGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCTTATc       <  1:1686998/69‑1 (MQ=255)
                                                             aCGCGCTTGCGTGGTCAGTAAGAAGTGACAGACTTTATATTCCACTTTTATTCCTTTTTATTCTTATcag     >  1:3699883/1‑68 (MQ=255)
                                                               gcgcTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCTTATCGATAg  <  1:475186/71‑1 (MQ=255)
                                                                cgcTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:2566049/1‑36 (MQ=255)
                                                                cgcTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:3326896/1‑36 (MQ=255)
                                                                cgcTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:1347972/1‑36 (MQ=255)
                                                                cgcTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:1724980/1‑36 (MQ=255)
                                                                cgcTTGCGTGGTCAGTAAGAAGTGCCAGACTTtata                                    >  1:2808182/1‑36 (MQ=255)
                                                                  |                                                                   
ATGTGTATCGAATGAAGCATAAGTGTAGTTTGCTTTCATCGCTCTCATCATGCTGTACGAAACGCGCTTGCGTGGTCAGTAAGAAGTGCCAGACTTTATATTCCACTTTTATTCCTTTTTATTCTTATCGATAG  >  W3110S.gb/3800221‑3800354

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: