Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 902,619 A→C 7.0% intergenic (+45/+31) trxB → / ← nadC thioredoxin reductase/quinolinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947902,6190AC 7.0% 166.6 / 2.3 55intergenic (+45/+31)trxB/nadCthioredoxin reductase/quinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base C (2/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTGAC  >  NC_002947/902479‑902763
                                                                                                                                            |                                                                                                                                                
gTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCgg                                                                                                                                          <  3:291321/149‑1 (MQ=255)
     ccGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttttt                                                                                                                                     >  7:15557/1‑149 (MQ=255)
       gCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTa                                                                                                                                   <  1:11085/149‑1 (MQ=255)
        ccAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTAt                                                                                                                                  >  6:189242/1‑149 (MQ=255)
        ccAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTAt                                                                                                                                  >  8:40151/1‑149 (MQ=255)
        ccAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTAt                                                                                                                                  >  2:313029/1‑149 (MQ=255)
         cAGGCCATCCCTTCGGCCTGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTGTGGCCAACGCCTCGTTCTGCGCTGAGCAACAGGGGTGCAACGCGCCCCCAATCAAACACCGGCGCTGGGCCGGTTTTTTAtt                                                                                                                                 <  1:91932/149‑1 (MQ=255)
         cAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGACAAAGGGGTGCAAAGCGCCCCCAATAAAAACCCGGCGCTTGGCCGGTTTTTTAtt                                                                                                                                 <  1:87802/149‑1 (MQ=255)
               aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                       >  4:319405/1‑137 (MQ=255)
               aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                       <  3:319405/137‑1 (MQ=255)
                tCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGt                                                                                                                                         <  4:250492/134‑1 (MQ=255)
                tCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGt                                                                                                                                         >  3:250492/1‑134 (MQ=255)
                  aCTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAAt                                                                                                                        <  5:36728/149‑1 (MQ=255)
                      cGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCa                                                                                                                      <  1:313029/147‑1 (MQ=255)
                            gTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                              <  8:15557/149‑1 (MQ=255)
                            gTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                              <  1:322598/149‑1 (MQ=255)
                               cGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCt                                                                                                           <  5:189242/149‑1 (MQ=255)
                                gggCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCt                                                                                                           >  7:9597/1‑148 (MQ=255)
                                      cATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTc                                                                                                      <  5:332424/147‑1 (MQ=255)
                                      cATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGc                                                                                                    <  8:9597/149‑1 (MQ=255)
                                                cTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGt                                                                                          <  5:72371/149‑1 (MQ=255)
                                                          aGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCtt                                                                                 >  1:203249/1‑148 (MQ=255)
                                                                aCCTGGACGGTTTGGCCAACGCCTCGTTGTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATc                                                                          >  1:186073/1‑149 (MQ=255)
                                                                 ccTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGCCGCATCGACAGGTCCACCGCTTTCACATcc                                                                         >  8:255945/1‑149 (MQ=255)
                                                                  cTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                              >  5:184410/1‑111 (MQ=255)
                                                                  cTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAg                                                                                                              <  6:184410/111‑1 (MQ=255)
                                                                   tGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCtt                                                                       >  3:276222/1‑149 (MQ=255)
                                                                                cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGt                                                                                                       >  8:289780/1‑104 (MQ=255)
                                                                                cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGt                                                                                                       <  7:289780/104‑1 (MQ=255)
                                                                                cAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCAcc                                                           >  2:152822/1‑148 (MQ=255)
                                                                                    gCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCg                                                      <  2:203249/149‑1 (MQ=255)
                                                                                      cTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATc                                                       >  4:166721/1‑146 (MQ=255)
                                                                                               atctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGt                                             <  4:276222/147‑1 (MQ=255)
                                                                                                tctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCa                                          >  7:89417/1‑149 (MQ=255)
                                                                                                tctgaGAAAAAGGGGTGCAAAGCCCCCCCAATAAAAAACCGGCAAAGCGCCCGTATTTTATTGCCTGGAATCAGCTCACAGGCTCAGTCGCATCGACACGTACACAGCTTTAAACTCCATCCTCAACGCACCAAACGAGATGTAGTcc                                           >  8:253819/1‑148 (MQ=255)
                                                                                                   gagaAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGt                                                                    >  5:241766/1‑120 (MQ=255)
                                                                                                   gagaAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGt                                                                    <  6:241766/120‑1 (MQ=255)
                                                                                                   gagaAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCAccc                                       >  8:292200/1‑149 (MQ=255)
                                                                                                      aaaaaGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGCGATGTAGTCCACCCCCg                                    >  4:271372/1‑149 (MQ=255)
                                                                                                        aaaGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGtc                                  >  4:175486/1‑149 (MQ=255)
                                                                                                            gggTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCgg                              <  7:292200/149‑1 (MQ=255)
                                                                                                             ggTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTa                                               >  3:210142/1‑131 (MQ=255)
                                                                                                             ggTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTa                                               <  4:210142/131‑1 (MQ=255)
                                                                                                                    aGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGGTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGCCAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCCCCCCCGTCTCGGCAATcccc                      >  7:36988/1‑149 (MQ=255)
                                                                                                                           gcAATAAAAAACCGGCGCTTGGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAggg               >  6:244010/2‑149 (MQ=39)
                                                                                                                             aaTAAAAAACCGTCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGAt                                                  <  5:1110/112‑1 (MQ=255)
                                                                                                                             aaTAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGAt                                                  >  6:1110/1‑112 (MQ=255)
                                                                                                                               tAAAAAACCGGCGCTTGGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgt           >  4:24550/1‑149 (MQ=255)
                                                                                                                                aaaaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATc                                                       <  4:106626/104‑1 (MQ=255)
                                                                                                                                aaaaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATc                                                       >  3:106626/1‑104 (MQ=255)
                                                                                                                                 aaaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgttt         >  7:151646/1‑149 (MQ=255)
                                                                                                                                 aaaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgttt         >  4:97679/1‑149 (MQ=255)
                                                                                                                                   aaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttcg       <  3:166721/149‑1 (MQ=255)
                                                                                                                                    aaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttcgt      <  8:89417/149‑1 (MQ=255)
                                                                                                                                        ggCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTGAc  >  6:247427/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTGAC  >  NC_002947/902479‑902763

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