Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 902,621 G→T 10.1% intergenic (+47/+29) trxB → / ← nadC thioredoxin reductase/quinolinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947902,6210GT10.1% 163.9 / 3.6 54intergenic (+47/+29)trxB/nadCthioredoxin reductase/quinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (22/26);  new base T (3/3);  total (25/29)
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

CACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTGACCCCG  >  NC_002947/902476‑902767
                                                                                                                                                 |                                                                                                                                                  
cACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgcgc                                                                                                                                                 >  4:203863/1‑149 (MQ=255)
cACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGCGCTGCATGGCGCCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAgcac                                                                                                                                                 >  2:266527/1‑147 (MQ=255)
 aCGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCTTGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAATGCGCCCCCAATAAAAAACCGGCCAAGCGcc                                                                                                                                                <  1:28781/149‑1 (MQ=255)
   gTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCgg                                                                                                                                              >  7:133517/1‑149 (MQ=255)
      tACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                           >  4:197274/1‑149 (MQ=255)
      tACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                           >  4:5222/1‑149 (MQ=255)
       aCCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGttt                                                                                                                                           <  7:276796/148‑1 (MQ=255)
         cGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGtttttt                                                                                                                                        <  8:133517/149‑1 (MQ=255)
                ccATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCt                                                                                                                                 <  3:203863/149‑1 (MQ=255)
                  aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGcc                                                                                                                                  <  2:206652/146‑1 (MQ=255)
                  aTCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGcc                                                                                                                                  >  1:206652/1‑146 (MQ=255)
                     aCTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCTAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAAt                                                                                                                            <  3:5222/149‑1 (MQ=255)
                     aCTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAAt                                                                                                                            <  3:197274/149‑1 (MQ=255)
                              ggTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTcaca                                                                                                                   <  2:212172/149‑1 (MQ=255)
                                    ggCTGCTTGGGGGGGCTGTTGTTCGTGGGCTACCTGTCCGGTTTGGCCAACGCCGCGTTCTTATATGATAAAAAGGGGTCCAAAGCGCCCCAAATACAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCa                                                                                                             <  1:266527/149‑1 (MQ=255)
                                      cTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAgg                                                                                                           <  7:227593/149‑2 (MQ=255)
                                                 cgcTGGATGGCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGCGCGCATAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAg                                                                                                <  5:271053/149‑1 (MQ=255)
                                                    tGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATAGACAGGTc                                                                                             >  2:175889/1‑149 (MQ=255)
                                                      gATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGc                                                                                                                        >  6:235757/1‑120 (MQ=255)
                                                      gATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGc                                                                                                                        <  5:235757/120‑1 (MQ=255)
                                                       aTGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCAc                                                                                          <  5:119024/149‑1 (MQ=255)
                                                           cGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGtttttt                                                                                                                                        <  6:195154/99‑1 (MQ=38)
                                                           cGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCGCTTGGCCGGtttttt                                                                                                                                        >  5:195154/1‑99 (MQ=38)
                                                             aGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCttt                                                                                    <  1:175889/149‑1 (MQ=255)
                                                               cgcTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTc                                                                                                          <  4:229366/125‑1 (MQ=255)
                                                               cgcTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTc                                                                                                          >  3:229366/1‑125 (MQ=255)
                                                                        gACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTg                                                                          <  3:134394/148‑1 (MQ=255)
                                                                                             ttCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATg                                                     >  4:262389/1‑148 (MQ=255)
                                                                                                  atctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCAt                                                                     >  8:146515/1‑127 (MQ=255)
                                                                                                  atctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCAt                                                                     <  7:146515/127‑1 (MQ=255)
                                                                                                   tctgaGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCAt                                                                     >  8:72225/1‑126 (MQ=255)
                                                                                                      gagaAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAg                                                  >  6:760/1‑142 (MQ=255)
                                                                                                      gagaAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAg                                                  <  5:760/142‑1 (MQ=255)
                                                                                                         aaaaaGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCg                                        >  5:251746/1‑149 (MQ=255)
                                                                                                         aaaaaGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCg                                        >  4:31552/1‑149 (MQ=255)
                                                                                                             aGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATcc                                                                             <  4:332145/108‑1 (MQ=255)
                                                                                                               gggTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCACTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATg                                                     >  3:334118/1‑130 (MQ=255)
                                                                                                               gggTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCACTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATg                                                     <  4:334118/130‑1 (MQ=255)
                                                                                                                  tGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCaa                               <  3:262389/149‑1 (MQ=255)
                                                                                                                           cccccAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCa                      <  6:251746/149‑1 (MQ=255)
                                                                                                                              ccAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCg                                   <  2:149350/133‑1 (MQ=255)
                                                                                                                              ccAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCg                                   >  1:149350/1‑133 (MQ=255)
                                                                                                                              ccAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAggg                   >  3:58355/1‑149 (MQ=255)
                                                                                                                              ccAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAggg                   >  3:177737/1‑149 (MQ=255)
                                                                                                                               cAATAAAAAACCGGCGCTTGGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGt                  >  3:203341/1‑149 (MQ=39)
                                                                                                                                aaTAAAAAACCGGCGCTTGGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtc                 >  3:69249/1‑149 (MQ=255)
                                                                                                                                aaTAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCg                                   <  3:61931/131‑1 (MQ=255)
                                                                                                                                aaTAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCg                                   >  4:61931/1‑131 (MQ=255)
                                                                                                                                   aaaaaaCCGGCCAAGCGCCGGTTTTTTATTGACTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGTTCCACCGATTTCACATACTTGGTCATGGAAACAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCacg                    >  6:348572/1‑141 (MQ=255)
                                                                                                                                   aaaaaaCCGGACAATCGCCGGTTTTTTATTGCCTGCAATCATCTCACAGGCTCAGTCGCATCGACAGTTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCGGTCTCGGCAATCACCCGCAgg                    <  5:348572/143‑1 (MQ=255)
                                                                                                                                    aaaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgttt             >  3:42907/1‑149 (MQ=255)
                                                                                                                                     aaaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttc            >  4:259130/1‑149 (MQ=255)
                                                                                                                                      aaaCCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGtcgtttcg           <  4:42907/149‑1 (MQ=255)
                                                                                                                                         ccGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTg        <  3:259130/149‑1 (MQ=255)
                                                                                                                                              cAATCGGCGGTTTTTTATTGCCTGCAATCAGCTCACATGCTCATTCGCATCGACAGGTGCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTGAcccc   <  2:220283/149‑1 (MQ=255)
                                                                                                                                               aaGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTGACCccg  >  1:286906/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CACGTTTACCGCCAGGCCATCACTTCGGCCGGTGCGGGCTGCATGGCGGCGCTGGATGTCGAGCGCTACCTGGACGGTTTGGCCAACGCCTCGTTCTGATCTGAGAAAAAGGGGTGCAAAGCGCCCCCAATAAAAAACCGGCCAAGCGCCGGTTTTTTATTGCCTGCAATCAGCTCACAGGCTCAGTCGCATCGACAGGTCCACCGCTTTCACATCCTTGGTCATGGCACCAATCGAGATGTAGTCCACCCCCGTCTCGGCAATCACCCGCAGGGTCGTTTCGTTGACCCCG  >  NC_002947/902476‑902767

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