Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,666,747 A→T 13.9% intergenic (+21/‑45) acnA‑II → / → prpF aconitate hydratase/aconitate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,666,7470AT13.9% 76.5 / 6.2 36intergenic (+21/‑45)acnA‑II/prpFaconitate hydratase/aconitate isomerase
Reads supporting (aligned to +/- strand):  ref base A (12/19);  new base T (3/2);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

CAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGC  >  NC_002947/2666603‑2666889
                                                                                                                                                  |                                                                                                                                              
cAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCgg                                                                                                                                            <  3:160808/149‑1 (MQ=255)
 aaTGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGCGCGCCAAGCGGc                                                                                                                                           >  2:134892/1‑149 (MQ=255)
       gAGTGCCTGGAAGTACCGGTGACCTGCCGGCGGGACCCCGCCGAGGACGTGGCGATCAACGAAGCGGGCGGGGGGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAgg                                                                                                                                     <  8:290223/149‑1 (MQ=255)
        aGTGCCTGGAAGTACTGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCCCGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCTGGGGCCGCTTTGCGGcccccga                                                                                                                                     <  1:98491/149‑4 (MQ=255)
          tGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTa                                                                                                                                  <  5:267442/149‑1 (MQ=255)
             cTGGAAGTACCGGGGACCTGCCGTCTGGACACCGCCGAGTAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACGGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTAttt                                                                                                                               <  1:134892/149‑1 (MQ=255)
                gAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCa                                                                                                                                       <  5:249015/138‑1 (MQ=255)
                        ggTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTGAGCGGGCTAt                                                                                                                    >  3:26173/1‑149 (MQ=255)
                         gTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGa                                                                                                                        >  1:254923/1‑144 (MQ=255)
                         gTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGa                                                                                                                        <  2:254923/144‑1 (MQ=255)
                                         acacCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGCCCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCCGGGGATGTGCCCGGCCGAGCAGGATGCCCAAGaaaa                                                                                                   >  6:256284/1‑149 (MQ=255)
                                             cgccgAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGATCATGAATCGGGGGTCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGt                                                                                                <  1:278884/149‑1 (MQ=255)
                                                       gtgtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCCCAAAGCGGCCCCAGGGTATTTCACCGGCCGATAAGGATGTCCAACAAAATGTCACAATTAccc                                                                                     >  7:112110/1‑149 (MQ=255)
                                                         gtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGa                                                                                                      <  7:160246/130‑1 (MQ=255)
                                                         gtCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTAccc                                                                                     >  6:125673/1‑147 (MQ=255)
                                                                  cGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCaaaa                                                                           >  2:53598/1‑148 (MQ=255)
                                                                      gcgggcggGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATcccc                                                                      <  5:125673/149‑1 (MQ=255)
                                                                            ggggTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGGATTTCACCGGACTCTCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCAcc                                                                <  6:368914/149‑1 (MQ=255)
                                                                                          cTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTGCCCGGTCTATAAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGcggcgg                                                  >  7:91588/1‑149 (MQ=255)
                                                                                                  cAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCa                                                                              >  3:65684/1‑114 (MQ=255)
                                                                                                  cAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCa                                                                              <  4:65684/114‑1 (MQ=255)
                                                                                                          ccTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCaaa                                                                            >  5:315363/1‑108 (MQ=255)
                                                                                                          ccTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCaaa                                                                            <  6:315363/108‑1 (MQ=255)
                                                                                                            tCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTAACCGGACGATCAGGCTGCCCAAGAAAATGTCCCATGTACCCCAGCACAAAATCCCCGCCCCCTACATCCGCGGCGGCCCCAGCCAAGGcc                                    >  1:308013/1‑145 (MQ=255)
                                                                                                              gAAGCGACCGCCTGAGCATGACTCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCgtgt                                 <  8:248003/147‑1 (MQ=255)
                                                                                                              gAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGtt                                <  1:333084/148‑1 (MQ=255)
                                                                                                                      cGCCTGAGCATGAATCGGCGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCt                       <  4:26173/149‑1 (MQ=255)
                                                                                                                        ccTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCgcgg                                                     <  1:149590/117‑1 (MQ=255)
                                                                                                                        ccTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTAACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTg                      >  3:313415/1‑148 (MQ=255)
                                                                                                                                 gAATCGGGGGACGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCtcc            <  5:203277/149‑3 (MQ=255)
                                                                                                                                  aaTCGGGGGCCGCAAAGCGGACCCAGTGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGTCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGa                 <  5:298528/143‑1 (MQ=255)
                                                                                                                                    cctggggCCGCTTTGCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCt                       >  7:131355/4‑135 (MQ=255)
                                                                                                                                    cctggggCCGCTTTGCGGCCCCCGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCt                       <  8:131355/132‑1 (MQ=255)
                                                                                                                                    cctggggCCGCTTTGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGTACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGccc          >  7:337685/4‑148 (MQ=255)
                                                                                                                                    cctggggCCGCTTTGCGGCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTAGATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGccc          >  1:293151/4‑147 (MQ=255)
                                                                                                                                            ccGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGTCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCTGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGgc  <  8:112110/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CAATGGCGAGTGCCTGGAAGTACCGGTGACCTGCCGTCTGGACACCGCCGAGGAAGTGTCGATCTACGAAGCGGGCGGGGTGCTGCAGCGCTTTGCCCAGGACTTCCTCGAAGCGACCGCCTGAGCATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGC  >  NC_002947/2666603‑2666889

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: