Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,517 T→A 25.0% intergenic (+42/‑206) 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,5170TA25.0% 126.0 / 48.4 76intergenic (+42/‑206)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 (31/26);  new base A (9/10);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

ACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTG  >  NC_000913/2800419‑2800615
                                                                                                  |                                                                                                  
aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTg                                                                                                   <  2:446268/100‑1 (MQ=255)
aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTg                                                                                                   <  2:651574/100‑1 (MQ=255)
 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                                                                                          <  2:1642646/100‑1 (MQ=255)
         gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgggg                                                                                          <  1:260738/100‑1 (MQ=255)
          cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGTAAAAATGAGAATAAGTACGCGTCTTTGACCCCGAAGGCTGTCTTCgggg                                                                                          >  2:1162323/1‑99 (MQ=255)
          cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                                         <  2:1505879/100‑1 (MQ=255)
          cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                         <  2:2208440/100‑1 (MQ=255)
          cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGt                                                                                         <  1:1467126/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:91532/100‑1 (MQ=255)
                   tttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                                <  2:1059656/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                                 <  1:2175682/100‑1 (MQ=255)
                                                                  gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                 <  1:2435020/100‑1 (MQ=255)
                                                                  gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                 <  2:353649/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            >  2:191290/1‑100 (MQ=255)
                                                                                       aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg            >  2:1122268/1‑100 (MQ=255)
                                                                                       aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg            >  1:1615958/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)
                                                                                                  |                                                                                                  
ACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTG  >  NC_000913/2800419‑2800615

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: