Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,816,579 G→A intergenic (‑30/‑97) rph ← / → yicC truncated RNase PH/UPF0701 family protein YicC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,816,5790GA93.4% 684.1 / 24.2 199intergenic (‑30/‑97)rph/yicCtruncated RNase PH/UPF0701 family protein YicC
Reads supporting (aligned to +/- strand):  ref base G (3/10);  new base A (97/87);  total (101/98)
Fisher's exact test for biased strand distribution p-value = 4.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 201 total aligned reads displayed.

CTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGC  >  NC_000913/3816530‑3816626
                                                 |                                               
cTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTgg                                                 >  2:5227903/1‑50 (MQ=255)
cTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTgg                                                 <  2:1534343/50‑1 (MQ=255)
  aCGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                               >  2:2824777/1‑50 (MQ=255)
   cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGc                                                <  2:516273/48‑1 (MQ=255)
   cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                              >  2:5103216/1‑50 (MQ=255)
    ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                             >  1:1894419/1‑50 (MQ=255)
     gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                            <  1:1324721/50‑1 (MQ=255)
      ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcgc                                           <  1:5296033/50‑1 (MQ=255)
      ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                            >  2:1922368/1‑49 (MQ=255)
        tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                         >  2:4211408/1‑50 (MQ=255)
        tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                         >  2:1151312/1‑50 (MQ=255)
        tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                         >  1:4040325/1‑50 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCAtt                                        <  2:1179413/49‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                        <  1:2195197/49‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                        <  2:4882676/49‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                        <  1:4533369/49‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                        <  1:362337/49‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                       <  2:3249700/50‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                       <  1:2720057/50‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                       <  2:1234235/50‑1 (MQ=255)
          ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                       >  1:5357043/1‑50 (MQ=255)
           tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                       >  1:1367579/1‑49 (MQ=255)
           tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                       >  1:322172/1‑49 (MQ=255)
            ggACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                     >  2:3210783/1‑50 (MQ=255)
            ggACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                     >  1:1143485/1‑50 (MQ=255)
            ggACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                     >  1:1572513/1‑50 (MQ=255)
             gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                     <  1:3761031/49‑1 (MQ=255)
             gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                     <  1:3545370/49‑1 (MQ=255)
             gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                    >  2:4967230/1‑50 (MQ=255)
              aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                   <  2:2338173/50‑1 (MQ=255)
              aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                   >  2:4303730/1‑50 (MQ=255)
              aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                   <  1:5135135/50‑1 (MQ=255)
              aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                   >  1:2011349/1‑50 (MQ=255)
               cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACgg                                  <  2:500640/50‑1 (MQ=255)
               cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                   >  2:549851/1‑49 (MQ=255)
                gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                 <  2:5660020/50‑1 (MQ=255)
                gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                 <  2:3392279/50‑1 (MQ=255)
                gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                 >  1:5451726/1‑50 (MQ=255)
                 cATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGAc                                <  2:697466/50‑1 (MQ=255)
                 cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                                <  2:2833359/50‑1 (MQ=255)
                 cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                                <  2:2556079/50‑1 (MQ=255)
                 cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                                >  2:1723255/1‑50 (MQ=255)
                  atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                               >  2:3414107/1‑50 (MQ=255)
                  atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                               <  2:8663/50‑1 (MQ=255)
                  atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                               <  1:2867030/50‑1 (MQ=255)
                  atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                               <  1:2980811/50‑1 (MQ=255)
                  atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                               >  1:85468/1‑50 (MQ=255)
                   tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                              >  1:1936174/1‑50 (MQ=255)
                   tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTAGACGGACtt                              <  2:3586008/50‑1 (MQ=255)
                    atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                             >  2:2608826/1‑50 (MQ=255)
                    atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                             >  2:2506945/1‑50 (MQ=255)
                     ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                            <  2:1008842/50‑1 (MQ=255)
                     ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                            >  1:735844/1‑50 (MQ=255)
                     ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                             <  1:625376/49‑1 (MQ=255)
                       gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                          >  2:4500653/1‑50 (MQ=255)
                       gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                          >  1:1692123/1‑50 (MQ=255)
                        aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                         <  1:941428/50‑1 (MQ=255)
                        aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                         >  2:5411520/1‑50 (MQ=255)
                        aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                         <  1:1843514/50‑1 (MQ=255)
                        aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                          <  2:2163855/49‑1 (MQ=255)
                         aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                        <  1:3920155/50‑1 (MQ=255)
                          aTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcgg                       >  2:1825098/1‑50 (MQ=255)
                            ctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                      >  2:3455353/1‑49 (MQ=255)
                            ctcCGGCTTGAAACAAATGTGACTGCGCATTAGACGGACTTCCGGCGGt                      >  2:1805246/1‑49 (MQ=255)
                             tcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                    >  1:2189883/1‑50 (MQ=255)
                              ccGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                   >  2:3838101/1‑50 (MQ=255)
                               cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGtt                     <  1:4489947/47‑1 (MQ=255)
                               cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGtt                     <  2:3715533/47‑1 (MQ=255)
                               cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAtt                  <  2:4382077/50‑1 (MQ=255)
                               cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAt                   <  1:2593300/49‑1 (MQ=255)
                               cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTACCGGCGGtt                     <  2:3822034/47‑1 (MQ=255)
                                ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc                 >  1:3838052/1‑50 (MQ=255)
                                ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc                 >  1:459122/1‑50 (MQ=255)
                                 gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc                <  1:4958956/50‑1 (MQ=255)
                                    tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt             >  2:3498009/1‑50 (MQ=255)
                                    tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt             >  2:3107766/1‑50 (MQ=255)
                                     gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTaa             >  2:1912863/1‑48 (MQ=255)
                                     gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc            >  1:433617/1‑50 (MQ=255)
                                      aaaCAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc           >  2:4740827/1‑50 (MQ=255)
                                      aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc           >  1:3386579/1‑50 (MQ=255)
                                      aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc           <  2:5386768/50‑1 (MQ=255)
                                      aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc           >  2:528258/1‑50 (MQ=255)
                                       aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt          <  2:2986643/50‑1 (MQ=255)
                                       aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt          >  2:2910641/1‑50 (MQ=255)
                                       aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt          >  2:2060636/1‑50 (MQ=255)
                                       aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt          >  1:1058901/1‑50 (MQ=255)
                                        aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg         <  2:1280069/50‑1 (MQ=255)
                                         cAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa        >  1:2908919/1‑50 (MQ=255)
                                          aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa        >  1:3350746/1‑49 (MQ=255)
                                          aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa        >  1:5646752/1‑49 (MQ=255)
                                          aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc       <  2:5812268/50‑1 (MQ=255)
                                           aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc       >  1:5783754/1‑49 (MQ=255)
                                           aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc       >  1:3534579/1‑49 (MQ=255)
                                           aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACa      <  2:5260716/50‑1 (MQ=255)
                                            aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc       >  1:552467/1‑48 (MQ=255)
                                             tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg    >  1:4676546/1‑50 (MQ=255)
                                             tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg    >  2:2878647/1‑50 (MQ=255)
                                             tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg    <  1:3110688/50‑1 (MQ=255)
                                              gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg   >  1:4841409/1‑50 (MQ=255)
                                               tgGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGc  <  1:787104/50‑1 (MQ=255)
                                                 |                                               
CTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGGC  >  NC_000913/3816530‑3816626

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