Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 709,200 C→A 10.6% intergenic (+60/‑59) murJ → / → ribF lipid II flippase/bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947709,2000CA10.6% 134.7 / 4.2 47intergenic (+60/‑59)murJ/ribFlipid II flippase/bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (20/22);  new base A (2/3);  total (22/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 = 9.22e-01

GATCCTGGCGGGCGTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTCACCGCGGC  >  NC_002947/709053‑709348
                                                                                                                                                   |                                                                                                                                                    
gtggccgcggtgaaacccgtcgaagttgccaatggtggcgacacagccccggtgctcggggcgcaggttgtgaagacctcgaaccagctgcataacgcgcttcttgctcataaagtggccgattataaccacacctGGGCGCTGCTGCCt                                                                                                                                                    >  2:256672‑M1/137‑150 (MQ=255)
  tCCTGGCGGGCGTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGt                                                                                                                                                  <  1:801711/150‑1 (MQ=255)
           gcgTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGGTGACAGGCAGCAgcgc                                                                                                                                         <  1:450810/150‑1 (MQ=255)
                          ttGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTtata                                                                                                                          >  1:200417/1‑150 (MQ=255)
                             gCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATc                                                                                                                       <  2:492168/150‑1 (MQ=255)
                             gCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATc                                                                                                                       <  2:232075/150‑1 (MQ=255)
                             gCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATc                                                                                                                       <  2:232188/150‑1 (MQ=255)
                              cTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCg                                                                                                                      >  2:1104026/1‑150 (MQ=255)
                              cTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCg                                                                                                                      <  2:258548/150‑1 (MQ=255)
                               tGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTtata                                                                                                                          <  1:1023999/145‑1 (MQ=255)
                               tGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTtata                                                                                                                          >  2:1023999/1‑145 (MQ=255)
                                        cTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTAt                                                                                                            <  1:464480/150‑1 (MQ=255)
                                        cTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTAt                                                                                                            >  2:686717/1‑150 (MQ=255)
                                         tgtgtgGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATg                                                                                                           >  2:1075888/1‑150 (MQ=255)
                                           tgtgGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCACAGCACACGTCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAg                                                                                                         >  1:852044/1‑150 (MQ=255)
                                              ggCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCaa                                                                                                      >  1:909244/1‑150 (MQ=255)
                                                cTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTCTAATCGGCCACTTTATGAGCaaga                                                                                                    >  2:119031/1‑150 (MQ=255)
                                                 ttCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCaagaa                                                                                                   >  1:554778/1‑150 (MQ=255)
                                                      gCCCCGGCCTTTCGCGCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTCTTTGAGCGTCTCCGCCCTTGGGCGCTGCTGCTTGTCAAAAGCGCGGAGGTGGGGGTATAATCGGCCACTTTATGAGCAAGAAGcgcg                                                                                              <  1:810114/150‑1 (MQ=255)
                                                        cccGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGtt                                                                                            <  1:297641/150‑1 (MQ=255)
                                                        cccGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGtt                                                                                            >  1:434773/1‑150 (MQ=255)
                                                              aTTTCGCTCGCGTTGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGGTTTTTGCATTCCACGCCGCCCTTGGGCGCTGGTGCTTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAg                                                                                      <  1:498623/150‑1 (MQ=255)
                                                                  cgctcgcAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTgg                                                                                  <  1:555494/150‑1 (MQ=255)
                                                                     tcgcAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTc                                                                               >  1:541470/1‑150 (MQ=255)
                                                                        cAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAg                                                                            >  1:934861/1‑150 (MQ=255)
                                                                           ggCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTc                                                                         >  1:150544/1‑150 (MQ=255)
                                                                             ccTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCtt                                                                       <  2:434773/150‑1 (MQ=255)
                                                                               ttGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTca                                                                     <  2:540494/150‑1 (MQ=255)
                                                                               ttGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTca                                                                     <  2:957838/150‑1 (MQ=255)
                                                                                   aCTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAAc                                                                 <  2:909244/150‑1 (MQ=255)
                                                                                         gCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTgcgc                                                           >  1:239751/1‑150 (MQ=255)
                                                                                           ttttGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGccc                                                         >  1:1110823/1‑150 (MQ=255)
                                                                                                           gTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCtgt                                         <  2:92095/150‑1 (MQ=255)
                                                                                                              ggTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGc                                                                                                        >  2:149419/1‑84 (MQ=255)
                                                                                                              ggTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGc                                                                                                        <  1:149419/84‑1 (MQ=255)
                                                                                                              ggTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTc                                      >  1:373072/1‑150 (MQ=255)
                                                                                                                   ttCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGccac                                 <  1:372655/150‑1 (MQ=255)
                                                                                                                        gattataaccacacccGGGCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGa                                                                             >  2:145953‑M1/17‑101 (MQ=255)
                                                                                                                        gattataaccacacccGGGCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGa                                                                             <  1:145953‑M1/85‑1 (MQ=255)
                                                                                                                               accacacccGGGCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGc                                    >  1:87377‑M1/10‑135 (MQ=255)
                                                                                                                               accacacccGGGCGCTGGTGACAGGCAGCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGc                                    <  2:87377‑M1/126‑1 (MQ=255)
                                                                                                                                  gccCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGa                  <  1:1075888/150‑1 (MQ=255)
                                                                                                                                   cccTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGAc                 >  2:144588/1‑150 (MQ=255)
                                                                                                                                        GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCAtt                             <  2:589084‑M1/133‑1 (MQ=255)
                                                                                                                                        GGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCAtt                             >  1:589084‑M1/1‑133 (MQ=255)
                                                                                                                                          GCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCgacgggggtg          <  1:1084579‑M1/150‑4 (MQ=255)
                                                                                                                                          GCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTc          <  1:34817‑M1/150‑1 (MQ=255)
                                                                                                                                                  gcCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTCACCGCGGc  >  1:363593/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATCCTGGCGGGCGTCGTGACGTATTTTGGCTGCCTGTACCTGTGTGGCTTCCGGCCCCGGCATTTCGCTCGCAAGGCCTTGCACTGAGGCTTTTGGCGGGTCAGGCGTCGGTTTTTCGCATTCCACGCCGCCCTTGGGCGCTGCTGCCTGTCACCAGCGCCCGGGTGTGGTTATAATCGGCCACTTTATGAGCAAGAAGCGCGTTATGCAGCTGGTTCGAGGTCTTCACAACCTGCGCCCCGAGCACCGGGGCTGTGTCGCCACCATTGGCAACTTCGACGGGGTTCACCGCGGC  >  NC_002947/709053‑709348

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