Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,572,253 A→T 13.3% intergenic (‑72/+20) yicG ← / ← orn inner membrane protein/oligoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,572,2530AT13.3% 79.2 / 5.3 30intergenic (‑72/+20)yicG/orninner membrane protein/oligoribonuclease
Reads supporting (aligned to +/- strand):  ref base A (10/16);  new base T (2/2);  total (12/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.60e-01

TGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGG  >  NC_002947/5572106‑5572397
                                                                                                                                                   |                                                                                                                                                
tGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAat                                                                                                                                                <  2:729793/150‑1 (MQ=255)
               cAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCACCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACC‑AAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGa                                                                                                                                <  1:561847/150‑1 (MQ=255)
                    aGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTca                                                                                                                            <  2:450897/150‑1 (MQ=255)
                           cAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGGTACAATATCTGCGCCGTTTCGACTCACAGCTtg                                                                                                                     >  2:740350/1‑150 (MQ=255)
                                 tGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATAACCCACCAGATTCAGGTCGGTGGAGTAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                               <  1:542737/150‑1 (MQ=255)
                                            gTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATa                                                                                                    <  1:740087/150‑1 (MQ=255)
                                             tAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAg                                                                                                   <  1:258492/150‑1 (MQ=255)
                                                    gATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCa                                                                                            >  2:118186/1‑150 (MQ=255)
                                                          gTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                      >  1:633925/1‑150 (MQ=255)
                                                            agagCATCAACCACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcg                                                                                    <  2:621386/150‑1 (MQ=255)
                                                              agCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgat                                                                                  <  1:534976/150‑1 (MQ=255)
                                                                     acaacaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTgcgc                                                                                                                                        >  2:639873/1‑89 (MQ=255)
                                                                     acaacaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTgcgc                                                                                                                                        <  1:639873/89‑1 (MQ=255)
                                                                      caacaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTcgc                                                                          >  1:582352/1‑150 (MQ=255)
                                                                          aTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGat                                                                      >  2:33973/1‑150 (MQ=255)
                                                                                cGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGATGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTc                                                                <  1:245872/150‑1 (MQ=255)
                                                                                       gTGCAGGAAAGGCGCGCAGTATAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAgc                                                            >  2:325450/1‑147 (MQ=255)
                                                                                       gTGCAGGAAAGGCGCGCAGTATAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAgc                                                            <  1:325450/147‑1 (MQ=255)
                                                                                                         gTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                       <  2:582352/150‑1 (MQ=255)
                                                                                                            tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTg                                    >  1:378910/1‑150 (MQ=255)
                                                                                                              aCCAGTTGCCCAAGGAACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                  <  1:271646/150‑1 (MQ=255)
                                                                                                               ccAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCctt                                        <  2:256265/143‑1 (MQ=255)
                                                                                                               ccAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCctt                                        >  1:256265/1‑143 (MQ=255)
                                                                                                                 aGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCt                               <  2:659916/150‑1 (MQ=255)
                                                                                                                    tGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTc                            <  1:33973/150‑1 (MQ=255)
                                                                                                                              aCCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTcc                  >  1:613648/1‑150 (MQ=255)
                                                                                                                                        ggCGCAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                  >  2:218883/1‑124 (MQ=255)
                                                                                                                                        ggCGCAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                  <  1:218883/124‑1 (MQ=255)
                                                                                                                                         gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                       <  2:278913/118‑1 (MQ=255)
                                                                                                                                         gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttc                                       >  1:278913/1‑118 (MQ=255)
                                                                                                                                              gagtgaaTATCCGCTCCCTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTCGGAGCTGTCACGCACGTCCGGCGCCCAGCGCGCgg  <  2:80920/145‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGG  >  NC_002947/5572106‑5572397

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