Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,519 T→G 25.9% intergenic (+44/‑204) ygaM → / → nrdH putative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,5190TG25.9% 103.2 / 50.3 73intergenic (+44/‑204)ygaM/nrdHputative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base G (9/10);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGT  >  NC_000913/2800420‑2800618
                                                                                                   |                                                                                                   
cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGt                                                                                                     >  2:402588/1‑100 (MQ=255)
 agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTc                                                                                                    >  1:697682/1‑100 (MQ=255)
  gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                   >  1:282255/1‑100 (MQ=255)
    agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTc                                                                                                 >  2:1019698/1‑100 (MQ=255)
     gcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                                >  2:2483272/1‑100 (MQ=255)
      cTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgg                                                                                               >  2:703750/1‑100 (MQ=255)
       tGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCggg                                                                                              <  2:2308365/100‑1 (MQ=255)
       tGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCggg                                                                                              <  1:495445/100‑1 (MQ=255)
        gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCCGGCGGAAATTGTAGTGCAAAAATGGTAATAAATACGCAACTTTGACGCCGAAGCCTGTc                                                                                                    >  1:661417/1‑93 (MQ=255)
        gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgggg                                                                                             >  2:650420/1‑100 (MQ=255)
        gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgggg                                                                                             >  2:1399416/1‑100 (MQ=255)
        gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgggg                                                                                             <  1:260738/100‑1 (MQ=255)
        gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgggg                                                                                             <  2:1642646/100‑1 (MQ=255)
         cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGTAAAAATGAGAATAAGTACGCGTCTTTGACCCCGAAGGCTGTCTTCgggg                                                                                             >  2:1162323/1‑99 (MQ=255)
         cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                                            <  2:1505879/100‑1 (MQ=255)
         cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                            <  1:1467126/100‑1 (MQ=255)
         cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                            <  2:2208440/100‑1 (MQ=255)
         cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                            <  2:2054867/100‑1 (MQ=255)
         cGGTAGGTATTTTCATCAGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTAACCCCGAAGACAGGCTTCGGGGt                                                                                            <  1:2014393/100‑1 (MQ=255)
          ggtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGtt                                                                                           <  1:2321418/100‑1 (MQ=255)
          ggtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                            <  1:1043947/99‑1 (MQ=255)
               gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttt                                                                                      >  2:1832878/1‑100 (MQ=255)
               gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtcaaaga                                                                                      <  2:1460267/100‑7 (MQ=38)
                taTTTTTATCGGCGCACTGTTAAGTATGCGCAAATCGTAGTGCAAACATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                            <  1:1824601/93‑1 (MQ=255)
                taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtttt                                                                                     >  1:272964/1‑100 (MQ=255)
                 aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                    >  1:895970/1‑100 (MQ=255)
                 aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGAAAAAATGATAATAAATACGAGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttttt                                                                                    >  1:1410589/1‑100 (MQ=255)
                  tttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                                   <  2:1059656/100‑1 (MQ=255)
                  tttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                                   <  2:91532/100‑1 (MQ=255)
                     ttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAATTGATAATAAATACGCGTCTTTGGCCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGcc                                                                                 <  2:282258/99‑1 (MQ=255)
                      tATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATATTAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGcc                                                                                 <  1:2083809/98‑1 (MQ=255)
                       aTCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                              <  1:412542/100‑1 (MQ=255)
                          ggCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                            <  1:1107721/99‑1 (MQ=255)
                          ggCGCACTGATAAGCATGCGCAAATCGTAGTGGAAAACTGATAATAAATACGCGTCTTTGATCGCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                             >  1:61663/1‑98 (MQ=255)
                             gcACTGTTAAGCATGAGCAAATCGTAGTGCAAAAATTATAATAAATACGCGTCTTTGACCGCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAAt                                                                         <  1:2370225/99‑1 (MQ=255)
                                cTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCTTGTTTTATcaca                                                                     >  1:1281467/1‑100 (MQ=255)
                                cTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcaca                                                                     >  2:2035759/1‑100 (MQ=255)
                                       gCATGCGCAAATCGTAGTGCAAAAATGATAATAAATAGGCGTCTTTGACCCCGAAGCCTGTCTTCTGGGTTTCTTTTTGCCTGGTGAATCACaaaa                                                                  >  1:1384136/1‑96 (MQ=255)
                                        cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGAGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCGGGTGAATCACAAAAATccc                                                             <  1:1282949/100‑1 (MQ=255)
                                          tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGTGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccccc                                                           <  2:697685/100‑1 (MQ=255)
                                          tGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccccc                                                           >  2:584028/1‑100 (MQ=255)
                                           gcgcAAATCCTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAACAATccccc                                                           <  2:2162006/99‑1 (MQ=255)
                                            cgcAACTCGTAGTGCAAAAATGATAATAAATACAAGTCTTTGACCCCGAAACCTGTCTTCGGTGTTTCTTTTTGCCAGGTGAATCACAAAAATCACCCt                                                          >  1:181231/1‑99 (MQ=255)
                                                aaTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcccc                                                     <  2:458867/100‑1 (MQ=255)
                                                    gTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                 >  2:1206098/1‑100 (MQ=255)
                                                      aatacAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGTCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                               <  1:1669101/96‑1 (MQ=255)
                                                        tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                             >  2:2062373/1‑100 (MQ=255)
                                                         gCAAAAATGATAATAAATACGCGTCTTTGACCCCGACGCCTGTCTTCGGGGATTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                            >  1:452774/1‑100 (MQ=255)
                                                         gCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTc                                            >  1:1080515/1‑100 (MQ=255)
                                                         gCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGAGGTTTCTTTTTGCCTGGTGAATCACACAAATCCCCCTACCCCGTCACGCTc                                            >  1:2394027/1‑100 (MQ=255)
                                                            aaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCGTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCata                                         <  2:1816733/100‑1 (MQ=255)
                                                            aaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCata                                         <  2:2054604/100‑1 (MQ=255)
                                                                 gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                    <  2:353649/100‑1 (MQ=255)
                                                                 gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                    <  1:2435020/100‑1 (MQ=255)
                                                                 gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                    <  1:2175682/100‑1 (MQ=255)
                                                                  ataataAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAgg                                   <  2:2482146/100‑1 (MQ=255)
                                                                   taataaATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAggg                                  <  1:1758725/100‑1 (MQ=255)
                                                                    aataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGt                                 >  1:857233/1‑100 (MQ=255)
                                                                      taaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTaa                               >  2:350852/1‑100 (MQ=255)
                                                                             caaaaagaaaCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                        >  2:236381/10‑100 (MQ=37)
                                                                              cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTTAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCCGGGTAATTTCGAcc                       >  1:585850/1‑100 (MQ=255)
                                                                               gTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAcc                       <  2:2132413/99‑1 (MQ=255)
                                                                                     aaCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAttt                >  2:103316/2‑100 (MQ=255)
                                                                                      aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg               >  1:1615958/1‑100 (MQ=255)
                                                                                      aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg               >  2:191290/1‑100 (MQ=255)
                                                                                      aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg               >  2:1122268/1‑100 (MQ=255)
                                                                                      aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCCTCACGCTCATATCCCGGGTAATTTCGACCACTATTTg               >  1:1743761/1‑100 (MQ=255)
                                                                                       ccccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc              <  2:2405909/100‑1 (MQ=255)
                                                                                       ccccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt             >  1:1497378/1‑100 (MQ=255)
                                                                                       ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc              >  2:146372/1‑100 (MQ=255)
                                                                                        cccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTTGTGAATCACAAAAATCCCCCTACCCCGTCACGATCATATCCAGGGTAATTTCGAACACTATTTCCt             >  1:2071787/1‑100 (MQ=255)
                                                                                        cccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt             >  1:881504/1‑100 (MQ=255)
                                                                                               ccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATAtt      >  2:321037/1‑100 (MQ=255)
                                                                                                cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtg     >  1:140636/1‑100 (MQ=255)
                                                                                                   tCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCACATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgt  <  1:1863829/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGT  >  NC_000913/2800420‑2800618

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