Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,849,629 T→C L265L (TTA→CTA)  sppA → protease IV, a signal peptide peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,849,6290TC96.2% 89.8 / ‑6.3 26L265L (TTA→CTA) sppAprotease IV, a signal peptide peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/13);  minor base A (0/1);  total (12/14)
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.88e-01

CCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTC  >  NC_000913/1849387‑1849866
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
ccaACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCTCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCTCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCaaa                                                                                                                                                                                                                                       <  1:19730/251‑1 (MQ=255)
                      gTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAat                                                                                                                                                                                                                 >  1:113834/1‑251 (MQ=255)
                      gTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAat                                                                                                                                                                                                                 <  1:74957/251‑1 (MQ=255)
                       ttGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTAGCTTGAGGGTATAACCAAAACCGGTGGCGATCCCGCGaaa                                                                                                                                                                                                                  <  2:83002/249‑1 (MQ=255)
                         gCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGAATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAatat                                                                                                                                                                                                               <  2:167532/250‑1 (MQ=255)
                                  gCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGTTGCCGATACCGCGAAATTTGCACTgaaa                                                                                                                                                                                                      >  1:40648/1‑247 (MQ=255)
                                     gAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAAc                                                                                                                                                                                                   >  2:83097/1‑250 (MQ=255)
                                                            gcgTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAGTATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCg                                                                                                                                                                           <  2:89554/251‑1 (MQ=255)
                                                              gTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCg                                                                                                                                                                           >  1:161172/1‑249 (MQ=255)
                                                                                      gAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGAc                                                                                                                                                 <  1:105874/251‑1 (MQ=255)
                                                                                       aaCCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGAcc                                                                                                                                                <  2:124497/251‑1 (MQ=255)
                                                                                                 tCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTc                                                                                                                                       <  1:29904/250‑1 (MQ=255)
                                                                                                  cGTGATGAGATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAAGATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAAGCGAAAAAGCACTGACCAAAGAATTCgg                                                                                                                                     <  2:113879/251‑1 (MQ=255)
                                                                                                              tCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGCTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGAc                                                                                                                         >  1:23255/1‑251 (MQ=255)
                                                                                                                                          gCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGttat                                                                                             <  2:213845/251‑1 (MQ=255)
                                                                                                                                                 tGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacgat                                                                                       <  1:105054/250‑1 (MQ=255)
                                                                                                                                                   gTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAACCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTattacgattac                                                                                    <  2:113834/251‑1 (MQ=255)
                                                                                                                                                                      tCTGAATACTGTTGCCGATAAACGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGCAAAAGCACTGACCAACGAATTCGGCTGGCGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGc                                                                   >  2:96900/1‑249 (MQ=255)
                                                                                                                                                                        tGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGc                                                                   >  1:162271/1‑247 (MQ=255)
                                                                                                                                                                                         aaCCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATc                                                <  2:79849/249‑1 (MQ=255)
                                                                                                                                                                                          acaggCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGgt                                             >  1:178854/4‑251 (MQ=255)
                                                                                                                                                                                                aGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTGTCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTgtcgtc                                       >  1:391/1‑251 (MQ=255)
                                                                                                                                                                                                    cccTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCttt                                    <  2:161172/250‑1 (MQ=255)
                                                                                                                                                                                                            aGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGgc                           >  2:76295/1‑251 (MQ=255)
                                                                                                                                                                                                                            cgcgcAAGGGTTGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCg           >  2:105413/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ttGCTTGAGGGTCTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTc  >  2:236469/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             
CCAACGGTCTGTACTACAAATCGTTGCTGGATAAGCTGAAAGTTTCCACCCATGTGTTCCGCGTGGGTACGTATAAATCTGCCGTTGAACCGTTTATTCGTGATGATATGTCACCGGCAGCCCGCGAAGCTGACAGCCGCTGGATTGGTGAGCTGTGGCAAAACTATCTGAATACTGTTGCCGCTAACCGGCAGATCCCTGCTGAGCAGGTATTCCCTGGCGCGCAAGGGTTGCTTGAGGGTTTAACCAAAACCGGTGGCGATACCGCGAAATATGCACTGGAAAACAAGCTGGTCGATGCACTGGCATCGAGTGCGGAAATCGAAAAAGCACTGACCAAAGAATTCGGCTGGAGTAAGACTGATAAAAATTATCGCGCCATCAGTTATTACGATTACGCATTGAAAACGCCGGCAGATACCGGTGACAGCATCGGTGTCGTCTTTGCTAATGGCGCAATTATGGATGGCGAGGAAACTC  >  NC_000913/1849387‑1849866

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