Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 529,932 A→G 11.4% intergenic (+26/‑158) PP_5SC → / → birA 5S ribosomal RNA/bifunctional biotin‑[acetyl‑CoA‑carboxylase] ligase/biotin operon repressor BirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947529,9320AG11.4% 90.8 / 3.5 35intergenic (+26/‑158)PP_5SC/birA5S ribosomal RNA/bifunctional biotin‑[acetyl‑CoA‑carboxylase] ligase/biotin operon repressor BirA
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base G (2/2);  total (17/18)
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.71e-01

GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGG  >  NC_002947/529792‑530054
                                                                                                                                            |                                                                                                                          
taaggcaggcatgatagccatgggctgttggcttgcataggaaaagggcctgcaggagaggccgcgtatcgcaatcacagggaggagcattggcaaggccaggatttttaagcccacaaaggCAAAACCCCTATCTGCGCATGCAGGTAg                                                                                                                   >  1:332409‑M1/123‑150 (MQ=255)
  aggcaggcatgatagccatgggctgttggcttgcataggaaaagggcctgcaggagaggccgcgtatcgcaatcacagggaggagcattggcaaggccaggatttttaagcccacaaagaCAAAACCCCTATCTGCGCATGCAGGTAggg                                                                                                                 >  1:159458‑M1/121‑150 (MQ=255)
                   aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGc                                                                                                >  1:379608/1‑150 (MQ=32)
                        ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTaaa                                                                                           >  1:482877/1‑150 (MQ=37)
                         cTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCACAATCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTaaaa                                                                                          >  1:346447/1‑150 (MQ=255)
                          tGATCCCATCCCGATCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTaaaaa                                                                                         <  2:253563/150‑1 (MQ=38)
                             tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATcc                                                                                      <  1:254883/150‑1 (MQ=40)
                                    ccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCtt                                                                               <  2:307969/150‑1 (MQ=255)
                                              tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGctc                                                                     <  2:308535/150‑1 (MQ=255)
                                                      cGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtgt                                                             <  2:379608/150‑1 (MQ=255)
                                                        aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGa                                                           <  2:482877/150‑1 (MQ=255)
                                                         tGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGAt                                                          >  2:105145/1‑150 (MQ=255)
                                                                ccGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATa                                                   >  1:225651/1‑150 (MQ=255)
                                                                     ggTAGTGTGGGGTTTCCCCATGTGAGAGTAAGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCgg                                              <  2:408823/150‑1 (MQ=255)
                                                                              gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAAt                                                                         <  2:239675/114‑1 (MQ=255)
                                                                              gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAAt                                                                         >  1:239675/1‑114 (MQ=255)
                                                                                 tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAg                                  >  2:148855/1‑150 (MQ=255)
                                                                                         tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCtttt                          <  2:346447/150‑1 (MQ=255)
                                                                                          gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc                         <  1:105145/150‑1 (MQ=255)
                                                                                              gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt                     <  2:341253/150‑1 (MQ=255)
                                                                                                    gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTg                                                        >  1:117684/1‑109 (MQ=255)
                                                                                                    gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTg                                                        <  2:117684/109‑1 (MQ=255)
                                                                                                    gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTcc                        <  2:372756/141‑1 (MQ=255)
                                                                                                    gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTcc                        >  1:372756/1‑141 (MQ=255)
                                                                                                    aggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc                         >  1:512191‑M1/23‑140 (MQ=255)
                                                                                                    aggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc                         <  2:512191‑M1/118‑1 (MQ=255)
                                                                                                           tCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGc        >  1:577710/1‑150 (MQ=255)
                                                                                                           tCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCACACAGc        >  2:95925/1‑150 (MQ=255)
                                                                                                              aGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCa     >  2:502787/1‑150 (MQ=255)
                                                                                                               gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCAt    <  2:577710/150‑1 (MQ=255)
                                                                                                                 ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATgg  <  1:148855/150‑1 (MQ=255)
                                                                                                                 ccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCt                                      >  2:162320‑M1/10‑114 (MQ=255)
                                                                                                                 ccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCt                                      <  1:162320‑M1/105‑1 (MQ=255)
                                                                                                                     tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt                     <  1:59986/127‑1 (MQ=255)
                                                                                                                     tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt                     >  2:59986/1‑127 (MQ=255)
                                                                                                                                            |                                                                                                                          
GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGG  >  NC_002947/529792‑530054

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