Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,256,639 C→G 14.8% A325A (GCC→GCG)  mrp → ATP‑dependent Fe‑S cluster transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,256,6390CG14.8% 71.2 / 4.2 27A325A (GCC→GCG) mrpATP‑dependent Fe‑S cluster transferase
Reads supporting (aligned to +/- strand):  ref base C (12/11);  new base G (2/2);  total (14/13)
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.87e-01

CTCGAATTGCGGCCATGCCGAACACCTGTTCGGTGAGGGTGGCGGCGAGAAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGC  >  NC_002947/1256492‑1256782
                                                                                                                                                    |                                                                                                                                               
cTCGAATTGCGGCCATGCCGAACACCTGTTCGGTGAGGGTGGCGGCGAGAAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGa                                                                                                                                               >  1:292175/1‑150 (MQ=255)
          ggCCATGCCGAACACCTGTTCGGTGAGGGTGGCGGCGAGAAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGc                                                                                                                                     >  2:175915/1‑150 (MQ=255)
              tcaccGAACAGGTGTTCGGTGAGGGTGGCGGCGAGAAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGa                                                                                                                                 >  1:447620/4‑150 (MQ=255)
                          tGTTCGGTGAGGGTGGCGACGAGAAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATAAGTGCACAGGCCGACAGGGGCCCGTGTCTTTTTAACACCGAGCCGGAAAGCCAGAGTGCCATGGttt                                                                                                                     >  2:374620/1‑150 (MQ=255)
                               ggTGAGGGTGGCGACGAGAAGCTGGCCAGCCAGTATGCGAGT‑TACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAg                                                                                                                <  1:374620/150‑1 (MQ=255)
                                              gagaAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCTGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAg                                                                                                 <  2:121014/150‑1 (MQ=255)
                                              gagaAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAg                                                                                                 >  2:249345/1‑150 (MQ=255)
                                                         cagccagTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCac                                                                                                                          <  2:390553/114‑2 (MQ=255)
                                                         cagccagTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCac                                                                                                                          >  1:390553/1‑113 (MQ=255)
                                                             cagTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCa                                                                                                  <  1:215786/134‑1 (MQ=255)
                                                             cagTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCa                                                                                                  >  2:215786/1‑134 (MQ=255)
                                                               gTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACAAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGt                                                                                >  1:222276/1‑150 (MQ=255)
                                                                   g‑gCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCGGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGtgctg                                                                            >  1:353984/1‑150 (MQ=255)
                                                                   g‑gCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCGGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGtgctg                                                                            >  1:353905/1‑150 (MQ=255)
                                                                     gCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGtgctgc                                                                           >  2:276995/1‑150 (MQ=255)
                                                                      cGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCa                                                                          <  1:249345/150‑1 (MQ=255)
                                                                               gcgggCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTAACAGGAGCTGGCCCGCCCGGTAGGCGCGAGGATCGTTCTGCAGCAAGCTgc                                                                 <  2:264666/147‑1 (MQ=255)
                                                                                   gCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGaa                                                                      <  1:90205/141‑1 (MQ=255)
                                                                                   gCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGaa                                                                      >  2:90205/1‑141 (MQ=255)
                                                                                     gTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGcc                                                           >  2:269976/1‑150 (MQ=255)
                                                                                               ttGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATAccc                                                 <  1:269976/150‑1 (MQ=255)
                                                                                                tGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCa                                                <  2:292175/150‑1 (MQ=255)
                                                                                                           atgcGTGAGCGGGCCGACAGCGGAAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATc                                     <  2:222276/147‑1 (MQ=255)
                                                                                                                                       ccACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGtt         <  1:276995/150‑1 (MQ=255)
                                                                                                                                        cACGGCGATCGCGGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGttt        <  2:353905/150‑1 (MQ=255)
                                                                                                                                        cACGGCGATCGCGGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGttt        <  2:353984/150‑1 (MQ=255)
                                                                                                                                              gATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGc  >  1:23942/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTCGAATTGCGGCCATGCCGAACACCTGTTCGGTGAGGGTGGCGGCGAGAAGCTGGCCAGCCAGTATG‑GCGTCGACCTGCTGGCGTCGCTGCCGTTGTCGATGCTGATACGTGAGCAGGCCGACAGCGGCAAGCCCACGGCGATCGCCGAGCCGGAAAGCCAGATTGCCATGGTTTACCAGGAGCTGGCCCGCCAGGTAGGCGCGCGGATCGTGCTGCAGGAAGCTGCGGCGCCGGCAATGCCCAGCATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGC  >  NC_002947/1256492‑1256782

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: ‑