Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,865,592 C→A 12.5% R109L (CGC→CTC)  smc ← chromosome partition protein Smc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,865,5920CA12.5% 84.2 / 3.3 32R109L (CGC→CTC) smcchromosome partition protein Smc
Reads supporting (aligned to +/- strand):  ref base C (15/13);  new base A (2/2);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

GGGCTTGGCCTCGATCAGCTTGGAGATCATGCCCTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGCGGCCACTGGAGCCGTTGAAGATGACGTCGGTCATCGAC  >  NC_002947/4865447‑4865741
                                                                                                                                                 |                                                                                                                                                     
gggCTTGGCCTCGATCAGCTTGGAGATCATGCCCTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTTATCTCACGCCGCCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTATGCCCGTCGCGGGt                                                                                                                                                   >  2:347287/1‑150 (MQ=255)
         cTCGATCAGCTTGGAGATCATGCCCTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGAGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTgcg                                                                                                                                          >  1:83143/1‑150 (MQ=255)
            gATCAGCTTGGAGATCATGCCCTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTgcggcg                                                                                                                                       >  2:443172/1‑150 (MQ=255)
              tCAGCTTGGAGATCATGCCCTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGa                                                                                                                                     <  1:443172/150‑1 (MQ=255)
                               cccTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAg                                                                                                                    <  1:134101/150‑1 (MQ=255)
                                ccTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCg                                                                                                                          <  1:78873/143‑1 (MQ=255)
                                ccTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCg                                                                                                                          >  2:78873/1‑143 (MQ=255)
                                     tCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGCGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCa                                                                                                              >  1:341507/1‑150 (MQ=255)
                                              gAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGAGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCg                                                                                                     >  1:178126/1‑150 (MQ=255)
                                               aGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGa                                                                                                    >  1:230980/1‑150 (MQ=255)
                                               aGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGa                                                                                                    <  1:456672/150‑1 (MQ=255)
                                                   gCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAg                                                                                                >  2:514249/1‑150 (MQ=255)
                                                    cTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGAGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAgg                                                                                               <  2:83143/150‑1 (MQ=255)
                                                       cgcgGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGAGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGggtg                                                                                            <  2:178126/150‑1 (MQ=255)
                                                            ccccAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTGGAGGTAATAGCGGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGACATCTCGGCGTAGGCGGCATACTCGCCGACCACGGGGGTTTc                                                                                       >  2:114061/1‑150 (MQ=255)
                                                                     cGTCCCAAGGAAAATATCCGTGACGTACCCCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTc                                                                              <  1:564974/150‑1 (MQ=255)
                                                                          cAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATTCTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAaca                                                                         >  2:497202/1‑150 (MQ=255)
                                                                                        gTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATActg                                                           <  2:230980/150‑1 (MQ=255)
                                                                                                       cggcACTTGGTGCCGTTGAGGGAATATCTGTTCTGCGCGGCGCGTGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCATGGTTGTTTCGCTGTTGTCGAACACCATCGCGATACTGGATTAGCTGCCCgcc                                            >  2:344743/1‑148 (MQ=255)
                                                                                                                gTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGAACGGCTTGCGGcgc                                   >  2:679515/1‑148 (MQ=255)
                                                                                                                 tGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGCGGCCAc                                  >  2:374457/1‑150 (MQ=255)
                                                                                                                  gCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGc                                        >  1:305569/1‑143 (MQ=255)
                                                                                                                  gCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGc                                        <  2:305569/143‑1 (MQ=255)
                                                                                                                     gatgcggTAAGAGCTGTTCTGCCAGTCGCGTTTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTGGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGCGGCCACTGGa                              <  1:344743/145‑1 (MQ=255)
                                                                                                                     gTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTCGCCGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGCGGCCACTGGa                              <  1:347287/150‑1 (MQ=255)
                                                                                                                             aaTAGCTGTTCTGCCCGTCGAGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCtg                                                                                   >  1:432716/1‑89 (MQ=255)
                                                                                                                             aaTAGCTGTTCTGCCCGTCGAGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCtg                                                                                   <  2:432716/89‑1 (MQ=255)
                                                                                                                              aTAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGtt                                                                                         >  2:91520/1‑82 (MQ=255)
                                                                                                                              aTAGCTGTTCTGCCCGTCGCGGGTAACCTGGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGtt                                                                                         <  1:91520/82‑1 (MQ=255)
                                                                                                                                 gCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGCGGCCACTGGAGCCGTTGAAGAt                  <  1:374457/150‑1 (MQ=255)
                                                                                                                                              tCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCtgtt                                                                                 <  2:100677/74‑1 (MQ=255)
                                                                                                                                              tCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCtgtt                                                                                 >  1:100677/1‑74 (MQ=255)
                                                                                                                                                 cgGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGCGGCCACTGGAGCCGTTGAAGATGACGTCGGTCATCGAc  <  1:514249/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGGCTTGGCCTCGATCAGCTTGGAGATCATGCCCTGCTCGATGATCGAGTAGCTGCGCGGCCCCAGACCCGTCCCAAGGAAAATATCCGTGATGTCACGCCGGCGGCACTTGGTGCCGTTGAGGTAATAGCTGTTCTGCCCGTCGCGGGTAACCTTGCGGCGGATCGAAATCTCGGCGTAGGCGGCATACTCGCCGACCAGGGTGGTTTCGCTGTTGTCGAACACCAGCTCGATACTGGCCTGGCTGACCGGCTTGCGGCCACTGGAGCCGTTGAAGATGACGTCGGTCATCGAC  >  NC_002947/4865447‑4865741

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