Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,098,519 A→G 6.5% intergenic (+116/+39) yggW → / ← yggM HemN family putative oxidoreductase/DUF1202 family putative secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,098,5190AG 6.5% 203.5 / 2.8 61intergenic (+116/+39)yggW/yggMHemN family putative oxidoreductase/DUF1202 family putative secreted protein
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (2/2);  total (31/30)
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.80e-01

GGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTA  >  NC_000913/3098394‑3098667
                                                                                                                             |                                                                                                                                                    
ggCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGc                                                                                                                                        <  3:234354/140‑1 (MQ=255)
ggCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGc                                                                                                                                        >  4:234354/1‑140 (MQ=255)
          aCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                            >  3:160967/1‑126 (MQ=21)
          aCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                            <  4:160967/126‑1 (MQ=21)
                            ggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCAcc                                                                                                   >  4:1093442/1‑149 (MQ=255)
                               gTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATc                                                                                                      >  4:969158/1‑143 (MQ=255)
                               gTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATc                                                                                                      <  3:969158/143‑1 (MQ=255)
                                   aCGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTgc                                                                                            >  1:792663/1‑149 (MQ=255)
                                        ttATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAc                                                                                                                                                    >  1:490568/1‑88 (MQ=255)
                                        ttATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAc                                                                                                                                                    <  3:449781/88‑1 (MQ=255)
                                        ttATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAc                                                                                                                                                    >  4:449781/1‑88 (MQ=255)
                                        ttATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAc                                                                                                                                                    <  2:490568/88‑1 (MQ=255)
                                         tATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCa                                                                                      <  3:1093442/149‑1 (MQ=255)
                                             cAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgttg                                                                                  >  4:1238320/1‑149 (MQ=255)
                                              aGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgttgt                                                                                 <  3:1238320/149‑1 (MQ=255)
                                               gCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATc                                                                                                                                            >  3:675794/1‑89 (MQ=255)
                                               gCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATc                                                                                                                                            <  4:675794/89‑1 (MQ=255)
                                               gCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgttgtt                                                                                >  3:17122/1‑149 (MQ=255)
                                                                      ggTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCttttt                                                                    >  4:1311832/1‑138 (MQ=255)
                                                                      ggTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCttttt                                                                    <  3:1311832/138‑1 (MQ=255)
                                                                        tCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCa                                                                >  3:1390378/1‑140 (MQ=255)
                                                                        tCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCa                                                                <  4:1390378/140‑1 (MQ=255)
                                                                           gATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCa                                                                >  4:190947/1‑137 (MQ=255)
                                                                           gATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCa                                                                <  3:190947/137‑1 (MQ=255)
                                                                            aTAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTg                                                                                                                             <  1:1334912/75‑1 (MQ=255)
                                                                            aTAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTg                                                                                                                             >  2:1334912/1‑75 (MQ=255)
                                                                            aTAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTc                                                          <  4:361579/142‑1 (MQ=255)
                                                                            aTAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTc                                                          >  3:361579/1‑142 (MQ=255)
                                                                            aTAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTa                                                      >  4:241850/1‑146 (MQ=255)
                                                                            aTAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTa                                                      <  3:241850/146‑1 (MQ=255)
                                                                              aaGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               >  1:1375770/1‑119 (MQ=255)
                                                                              aaGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                               <  2:1375770/119‑1 (MQ=255)
                                                                               aGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGc                                                                                                                            <  3:1431574/73‑1 (MQ=255)
                                                                               aGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGc                                                                                                                            >  4:1431574/1‑73 (MQ=255)
                                                                                 aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                              <  4:17122/149‑1 (MQ=255)
                                                                                     gcAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCa                                                                                                                           <  4:1459526/68‑1 (MQ=255)
                                                                                     gcAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCa                                                                                                                           >  3:1459526/1‑68 (MQ=255)
                                                                                                aTCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACtata                                                                                                                      >  3:1034939/1‑62 (MQ=255)
                                                                                                aTCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACtata                                                                                                                      <  4:1034939/62‑1 (MQ=255)
                                                                                                aTCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgttgtt                                                                                <  3:931655/100‑1 (MQ=255)
                                                                                                aTCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgttgtt                                                                                >  4:931655/1‑100 (MQ=255)
                                                                                                 tCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCg                                                                              >  1:717621/1‑101 (MQ=255)
                                                                                                 tCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCg                                                                              <  2:717621/101‑1 (MQ=255)
                                                                                                 tCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCttttt                                                                    >  4:2822/1‑111 (MQ=255)
                                                                                                 tCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCttttt                                                                    <  3:2822/111‑1 (MQ=255)
                                                                                                    gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttact                                                                                                           <  1:992490/69‑1 (MQ=255)
                                                                                                    gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttact                                                                                                           >  2:992490/1‑69 (MQ=255)
                                                                                                       gTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCg                                                                              <  1:1362512/95‑1 (MQ=255)
                                                                                                       gTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCg                                                                              >  2:1362512/1‑95 (MQ=255)
                                                                                                       gTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaaa                        >  1:1111676/1‑149 (MQ=255)
                                                                                                        ttGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                              <  3:118909/126‑1 (MQ=255)
                                                                                                        ttGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                              >  4:118909/1‑126 (MQ=255)
                                                                                                            ttGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAAt                                                                                                       <  3:1055148/65‑1 (MQ=255)
                                                                                                            ttGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAAt                                                                                                       >  4:1055148/1‑65 (MQ=255)
                                                                                                              gCCGGATTCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCCAACTTACTTAATCACCCCATTGCGCCGAAGTTGTTCGGTATCTTTTTTGCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGa                 >  1:570378/1‑149 (MQ=255)
                                                                                                                  gATGCGGCGTTGACCCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTaa                                                                                                        <  4:138147/58‑1 (MQ=25)
                                                                                                                  gATGCGGCGTTGACCCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTaa                                                                                                        >  3:138147/1‑58 (MQ=25)
                                                                                                                    tGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttttt           <  2:1111676/149‑1 (MQ=255)
                                                                                                                       ggggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttt            >  1:517907/1‑145 (MQ=255)
                                                                                                                       ggggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttt            <  2:517907/145‑1 (MQ=255)
                                                                                                                       ggggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTAc        <  2:792663/149‑1 (MQ=255)
                                                                                                                             aaCGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCATCGTTCATTATCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTa  >  1:337767/1‑149 (MQ=255)
                                                                                                                             aaCGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTa  >  3:40644/1‑149 (MQ=255)
                                                                                                                             aaCGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTa  >  2:817611/1‑149 (MQ=255)
                                                                                                                             |                                                                                                                                                    
GGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTA  >  NC_000913/3098394‑3098667

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