Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 25,774 A→C 16.1% intergenic (+73/‑52) lspA → / → fkpB prolipoprotein signal peptidase (signal peptidase II)/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091325,7740AC16.1% 103.9 / 11.5 37intergenic (+73/‑52)lspA/fkpBprolipoprotein signal peptidase (signal peptidase II)/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base A (16/15);  new base C (2/4);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGT  >  NC_000913/25627‑25921
                                                                                                                                                   |                                                                                                                                                   
gATACTGCCATCTGTGTCGGTGCGACACTGATTGTGCTGCAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAg                                                                                                                                                    <  3:1245587/149‑1 (MQ=255)
     tGCCATCTGTGTCGGTGCGGCACTGATAGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCGGCGta                                                                                                                                               >  3:1262757/1‑148 (MQ=255)
               gtCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                     <  4:264461/149‑5 (MQ=39)
               gtCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                     <  4:1262757/149‑5 (MQ=39)
                     gCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                               <  4:201761/149‑2 (MQ=38)
                                     tGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                   <  1:231240/145‑1 (MQ=255)
                                     tGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                   >  2:231240/1‑145 (MQ=255)
                                            ttttttGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAg                                                                                                                                                    >  4:590335/1‑105 (MQ=255)
                                            ttttttGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAg                                                                                                                                                    <  3:590335/105‑1 (MQ=255)
                                                           agCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGaa                                                                                         <  3:304934/149‑1 (MQ=255)
                                                                      aataatAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                >  1:651442/1‑99 (MQ=25)
                                                                      aataatAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                <  2:651442/99‑1 (MQ=25)
                                                                        taataaACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                          <  2:262639/87‑1 (MQ=255)
                                                                        taataaACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                          >  1:262639/1‑87 (MQ=255)
                                                                          ataaACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                        >  2:753230/1‑87 (MQ=255)
                                                                          ataaACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                        <  1:753230/87‑1 (MQ=255)
                                                                                                   aCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCa                                                    <  1:1603279/146‑1 (MQ=255)
                                                                                                   aCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCa                                                    >  2:1603279/1‑146 (MQ=255)
                                                                                                        tCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGcc                                                                    <  3:1808760/125‑1 (MQ=255)
                                                                                                        tCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGcc                                                                    >  4:1808760/1‑125 (MQ=255)
                                                                                                        tCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTaaaa                                            >  2:871223/1‑149 (MQ=255)
                                                                                                         cTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTc                                                                                             <  3:518254/99‑1 (MQ=255)
                                                                                                         cTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTGAAATATAAGAGCAAACCTGCATGTc                                                                                             >  4:518254/1‑99 (MQ=255)
                                                                                                          ttATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAAt                                                                                        <  3:1030112/103‑1 (MQ=255)
                                                                                                          ttATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAAt                                                                                        >  4:1030112/1‑103 (MQ=255)
                                                                                                                    tACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                   <  3:1837266/66‑1 (MQ=255)
                                                                                                                    tACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                   >  4:1837266/1‑66 (MQ=255)
                                                                                                                             ctgcGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcg                                        >  4:1830782/1‑132 (MQ=255)
                                                                                                                             ctgcGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcg                                        <  3:1830782/132‑1 (MQ=255)
                                                                                                                              tgcGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTaaaa                                            >  3:464714/1‑127 (MQ=255)
                                                                                                                              tgcGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTaaaa                                            <  4:464714/127‑1 (MQ=255)
                                                                                                                                       gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTAccc             <  1:871223/149‑1 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGca          >  2:276124/1‑149 (MQ=255)
                                                                                                                                               gATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTcc                                                                >  1:378184/1‑90 (MQ=255)
                                                                                                                                               gATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTcc                                                                <  2:378184/90‑1 (MQ=255)
                                                                                                                                                aTAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACg    >  1:74554/1‑149 (MQ=255)
                                                                                                                                                aTAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACg    >  3:1826456/1‑149 (MQ=255)
                                                                                                                                                  aaGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGt  >  3:571894/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGT  >  NC_000913/25627‑25921

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