Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 253,256 C→G 7.5% intergenic (‑41/+130) PP_0202 ← / ← PP_0203 CBS domain‑containing protein/dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947253,2560CG 7.5% 166.6 / 4.6 53intergenic (‑41/+130)PP_0202/PP_0203CBS domain‑containing protein/dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (21/28);  new base G (2/2);  total (23/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 = 8.60e-01

TTTCTCCGCCAGCAGCTTCAGGGCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTG  >  NC_002947/253108‑253403
                                                                                                                                                    |                                                                                                                                                   
tttCTCCGCCAGCAGCTTCAGGGCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAAc                                                                                                                                                     >  1:189190/1‑149 (MQ=255)
         cagcagCTTCAGGGCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATc                                                                                                                                            <  2:170293/149‑1 (MQ=255)
               cTTCAGGGCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGTTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTAc                                                                                                                                      <  2:245629/149‑1 (MQ=255)
                ttCAGGGCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACg                                                                                                                                     >  1:262750/1‑149 (MQ=255)
                   aGGGCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGgcg                                                                                                                                  <  4:114675/149‑1 (MQ=255)
                     ggCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGct                                                                                                                                >  3:8582/1‑149 (MQ=255)
                        gTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTg                                                                                                                             <  4:243227/149‑1 (MQ=255)
                             gCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACgg                                                                                                                        >  2:225618/1‑149 (MQ=255)
                                    gTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATg                                                                                                                   >  5:60951/1‑147 (MQ=255)
                                       aTCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTaaa                                                                                                                                                 <  1:162223/114‑1 (MQ=255)
                                       aTCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTaaa                                                                                                                                                 >  2:162223/1‑114 (MQ=255)
                                          cGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAAGTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAAc                                                                                                           >  8:92161/1‑149 (MQ=255)
                                          cGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAAGTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAAc                                                                                                           >  8:87357/1‑149 (MQ=255)
                                          cGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAAGTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAAc                                                                                                           <  7:92161/149‑1 (MQ=255)
                                          cGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAAGTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAAc                                                                                                           <  7:87357/149‑1 (MQ=255)
                                            ggCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACgg                                                                                                         <  2:262750/149‑1 (MQ=255)
                                                  tGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGa                                                                                                   <  4:8582/149‑1 (MQ=255)
                                                       tACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATg                                                                                                     >  3:285135/1‑142 (MQ=255)
                                                       tACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATg                                                                                                     <  4:285135/142‑1 (MQ=255)
                                                       tACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGa                                                                                              <  2:7161/149‑1 (MQ=255)
                                                         cacGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAgg                                                                                                 >  3:130376/1‑144 (MQ=255)
                                                         cacGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAgg                                                                                                 <  4:130376/144‑1 (MQ=255)
                                                         cacGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGAtt                                                                                            <  1:64463/149‑1 (MQ=255)
                                                                  gTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGTAGTAACTAAAGCATCCTGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCTATTGAAAAACGt                                                                                   >  5:351249/1‑149 (MQ=255)
                                                                        ggACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATgg                                                                              <  8:299748/148‑1 (MQ=255)
                                                                                  ttGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGc                                                                                                                                               <  6:67704/73‑1 (MQ=255)
                                                                                  ttGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGc                                                                                                                                               >  5:67704/1‑73 (MQ=255)
                                                                                    gatgATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGAtt                                                                                            <  6:259190/119‑1 (MQ=255)
                                                                                    gAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGAtt                                                                                            >  5:259190/1‑122 (MQ=255)
                                                                                    gAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCtt                                                                    <  1:225618/146‑1 (MQ=255)
                                                                                    gAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTg                                                                   >  1:8277/1‑147 (MQ=255)
                                                                                        aTTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCgg                                                             <  2:189190/149‑1 (MQ=255)
                                                                                            gTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCtt                                                         >  4:346234/1‑149 (MQ=255)
                                                                                             ttCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCttt                                                        >  4:199830/1‑149 (MQ=255)
                                                                                             ttCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCttt                                                        <  3:346234/149‑1 (MQ=255)
                                                                                             ttCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCtt                                                         <  6:60951/148‑1 (MQ=255)
                                                                                                gACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTaaa                                                                                                                                                 <  7:299889/57‑1 (MQ=255)
                                                                                                gACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTaaa                                                                                                                                                 >  8:299889/1‑57 (MQ=255)
                                                                                                gACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATg                                                                               <  4:322924/123‑1 (MQ=255)
                                                                                                gACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATg                                                                               >  3:322924/1‑123 (MQ=255)
                                                                                                gACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTcg                                                     >  4:154198/1‑149 (MQ=255)
                                                                                                 aCGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTcgc                                                    <  2:8277/149‑1 (MQ=255)
                                                                                                       ttCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTa                                               >  8:117861/1‑148 (MQ=255)
                                                                                                                   tctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAagag                                  >  8:324283/1‑149 (MQ=255)
                                                                                                                       ggggTGAAAAAAGCACTGTCTCGGAGAACCTACACCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCg                              <  2:104091/148‑1 (MQ=255)
                                                                                                                       cGGGTGAAAAAAGCACTGGCTCTGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                <  4:302690/147‑1 (MQ=255)
                                                                                                                          gTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTa                           <  3:199830/149‑1 (MQ=255)
                                                                                                                              aaaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTg                        <  7:117861/148‑1 (MQ=255)
                                                                                                                                           tcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGc             <  7:162271/146‑1 (MQ=255)
                                                                                                                                           tcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACAGAGAAGAGTCAGGTACTGCGCATCAAAga             >  8:162271/1‑145 (MQ=255)
                                                                                                                                                   aCTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCg         <  4:226777/142‑1 (MQ=255)
                                                                                                                                                   aCTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCg         >  3:226777/1‑142 (MQ=255)
                                                                                                                                                   aCTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCt   <  3:154198/148‑1 (MQ=255)
                                                                                                                                                   aCTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTg  >  3:165770/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTCTCCGCCAGCAGCTTCAGGGCGTCGAGCACCGAGTCATCCGGGCCGATGGTGTACACGGTCTGGTGCTGGGACTTGGTCTTGAGGATTTGTTCGACGGTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTG  >  NC_002947/253108‑253403

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