Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,515 C→A 26.2% intergenic (+40/‑208) 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,5150CA26.2% 129.5 / 45.9 74intergenic (+40/‑208)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 C (28/25);  new base A (9/10);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

ACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATT  >  NC_000913/2800419‑2800614
                                                                                                |                                                                                                   
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:1399416/1‑100 (MQ=255)
         gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgggg                                                                                         >  2:650420/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                                                                                        <  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: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           >  1:1615958/1‑100 (MQ=255)
                                                                                       aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg           >  2:1122268/1‑100 (MQ=255)
                                                                                       aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTg           >  2:191290/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)
                                                                                                |                                                                                                   
ACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATT  >  NC_000913/2800419‑2800614

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: