Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,409,146 G→A 11.1% intergenic (+84/‑191) PP_1230 → / → nadA hypothetical protein/quinolinate synthase iron‑sulfur cluster subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,409,1460GA11.1% 138.1 / 7.6 54intergenic (+84/‑191)PP_1230/nadAhypothetical protein/quinolinate synthase iron‑sulfur cluster subunit
Reads supporting (aligned to +/- strand):  ref base G (28/20);  new base A (3/3);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAGGGCCGGGTGATAC  >  NC_002947/1408999‑1409295
                                                                                                                                                   |                                                                                                                                                     
gACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGgg                                                                                                                                                     >  1:1141627‑M2/1‑148 (MQ=255)
     ttGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGcccac                                                                                                                                                >  2:606328/1‑148 (MQ=255)
     ttGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGc                                                                                                                                                >  2:1032164/1‑150 (MQ=255)
     ttGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGGGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGGTCACCCGCGAAGGGGCCCGc                                                                                                                                                >  1:738875/1‑150 (MQ=255)
          gaaTTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAg                                                                                                                                           <  2:940965/150‑1 (MQ=255)
            aTTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCa                                                                                                                                         <  1:473794/150‑1 (MQ=255)
                          tGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCt                                                                                                                           >  2:1196502/1‑150 (MQ=255)
                                    tCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCtt                                                                                                                          <  2:367203/141‑1 (MQ=255)
                                    tCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCtt                                                                                                                          >  1:367203/1‑141 (MQ=255)
                                                   gccgcAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACaaaaa                                                                                                                                <  2:997792/120‑1 (MQ=255)
                                                   gccgcAGGGGTGATGGTTTCTCCTGTTGCAGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACaaaaa                                                                                                                                >  1:997792/1‑120 (MQ=255)
                                                         ggggTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTGCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGATGGATTGTCTGTAGCCTTCg                                                                                            >  1:880322/1‑150 (MQ=255)
                                                            gggATTGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGcc                                                                                         >  1:347640/3‑150 (MQ=255)
                                                             tGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCt                                                                                        >  2:901785/1‑150 (MQ=255)
                                                               aTGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGTGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTcc                                                                                      >  2:1061908/1‑150 (MQ=255)
                                                                      ctcCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGataata                                                                               >  1:28171/1‑150 (MQ=255)
                                                                        ccTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATaataat                                                                             <  1:1196502/150‑1 (MQ=255)
                                                                         cTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATc                                                                            >  2:366670/1‑150 (MQ=255)
                                                                          tGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGCGGATTTCACAGTGTACGGGCGGGTTCACCCGCGAAGGGGCGCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTACGGCCCCCGCTACTAATcg                                                                           >  2:532815/1‑149 (MQ=38)
                                                                                ggCCTCTGCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGccc                                                                     <  2:1113341/150‑1 (MQ=255)
                                                                                 gCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCg                                                                    >  2:451510/1‑150 (MQ=255)
                                                                                   ctctTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCa                                                                  <  1:901785/150‑1 (MQ=255)
                                                                                    tctTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCag                                                                 >  2:177512/1‑150 (MQ=255)
                                                                                      ttCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttcaacg                                                                      <  2:385462‑M2/143‑82 (MQ=255)
                                                                                      ttCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttcaacg                                                                      >  1:385462‑M2/1‑62 (MQ=255)
                                                                                      ttCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCagag                                                               >  2:395038/1‑150 (MQ=255)
                                                                                             gTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCACGAAGGGGCCCCCGCAAGCCCTCCCACAACCCCTTTAGTAATCGCGGAATTGTCTGTAGCGTTCCGCCTCCGATAATATTCCGCGCCCGCAGACGGCCGta                                                        >  2:359854/1‑149 (MQ=255)
                                                                                                 aCCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGAggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttcaacggtcatttttc                                                            <  1:235935‑M2/142‑92 (MQ=255)
                                                                                                 aCCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGAggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttcaacggtcatttttc                                                            >  2:235935‑M2/1‑51 (MQ=255)
                                                                                                   ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttcaa                                                                        <  2:1159130‑M2/128‑80 (MQ=255)
                                                                                                   ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttcaa                                                                        >  1:1159130‑M2/1‑49 (MQ=255)
                                                                                                   ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgc                                                                                                    <  2:743882‑M2/100‑52 (MQ=255)
                                                                                                   ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgc                                                                                                    >  1:743882‑M2/1‑49 (MQ=255)
                                                                                                   ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGAggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgc                                                                                                    <  2:665039‑M2/100‑52 (MQ=255)
                                                                                                   ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGAggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgc                                                                                                    >  1:665039‑M2/1‑49 (MQ=255)
                                                                                                     gCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCaa                                                <  1:1032164/150‑1 (MQ=255)
                                                                                                     gCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGAggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttc                                                                          <  2:615024‑M2/124‑78 (MQ=255)
                                                                                                     gCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGAggccgcaacaggagaaaccatcacccctgcggcgacagggcccccagatgcgccgccatgaacatggcaaattcttc                                                                          >  1:615024‑M2/1‑47 (MQ=255)
                                                                                                     gCTCCGACGAGGATTTCACAGTGTAGGGGCGGGTCCACCCGCGCAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCaa                                                <  2:738875/150‑1 (MQ=255)
                                                                                                      cTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGTGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAAc                                               <  2:581630/150‑1 (MQ=255)
                                                                                                          tACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTTTGTACCCTTTTGTCTCGGGTAAGAAGGCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGc                                           <  1:511238/150‑1 (MQ=255)
                                                                                                             gAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGt                                        >  2:103947/1‑150 (MQ=255)
                                                                                                               ggATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTcc                                      <  2:44165/150‑1 (MQ=255)
                                                                                                                gATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCg                                     >  2:225517/1‑150 (MQ=255)
                                                                                                                        agagtAGGGGCCGGTTCACCCGCCAGGGGGCCCCCACAAGCACGACAAAAACCCCTTCAGTAATCGCTGGATGGCCTGTAGCCTTCGCCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCtt                             >  2:199521/4‑150 (MQ=255)
                                                                                                                          tgtAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCg                                                                                            <  2:174253/85‑1 (MQ=255)
                                                                                                                          tgtAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCg                                                                                            >  1:174253/1‑85 (MQ=255)
                                                                                                                           gtAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGc                          <  1:532815/150‑1 (MQ=255)
                                                                                                                              ggggCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGa                       >  1:226669/1‑150 (MQ=255)
                                                                                                                                       tCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAgg              <  1:395038/150‑1 (MQ=255)
                                                                                                                                           ccGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAGGGCCg          >  2:244456/1‑150 (MQ=255)
                                                                                                                                           ccGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAGGGCCg          >  1:611228/1‑150 (MQ=255)
                                                                                                                                                aaGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAGGGCCGGGTGa     <  1:177512/150‑1 (MQ=255)
                                                                                                                                                 aGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAGGGCCGGGTGAt    <  1:606328/150‑1 (MQ=255)
                                                                                                                                                   gggCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAGGGCCGGGTGATac  <  1:103947/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGATCTGCCAGGGCCGGGTGATAC  >  NC_002947/1408999‑1409295

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