Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,836,583 G→A 9.5% L226L (TTG→TTA)  PP_5115 → zinc protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,836,5830GA 9.5% 117.2 / 2.5 42L226L (TTG→TTA) PP_5115zinc protease
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base A (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCAGCCGCAGCGCGTACCTGGTACAAGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCTGCCAGGCAGGCCTACGCGC  >  NC_002947/5836443‑5836728
                                                                                                                                            |                                                                                                                                                 
aCTCCAGCCGCAGCGCGTACCTGGTACAAGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCg                                                                                                                                           >  7:161102/1‑149 (MQ=255)
    cAGCCGCAGCGCGTACCTGGTACAAGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGcc                                                                                                                                        <  8:257686/148‑1 (MQ=255)
    cAGCCGCAGCGCGTACCTGGTACAAGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGc                                                                                                                                         <  5:104743/147‑1 (MQ=255)
               cgTACCTGGTACAAGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAgg                                                                                                                            <  5:234946/149‑1 (MQ=255)
                          cAAGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTccc                                                                                                                 <  8:161102/149‑1 (MQ=255)
                           aaGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCtcctcc                                                                                                                  >  7:7374/1‑147 (MQ=255)
                                cTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGcc                                                                                                           <  8:7374/149‑1 (MQ=255)
                                        accCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAg                                                                                                               <  1:85424/137‑1 (MQ=255)
                                        accCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAg                                                                                                               >  2:85424/1‑137 (MQ=255)
                                          cccAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGtttt                                                                                                   >  6:58186/1‑147 (MQ=255)
                                              acaaTGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACg                                                                                             >  8:246121/1‑149 (MQ=255)
                                              acaaTGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACg                                                                                             <  5:58186/149‑1 (MQ=255)
                                                  tGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTgc                                                                                         <  7:246121/149‑1 (MQ=255)
                                                        ccTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTc                                                                                                  <  1:267015/134‑1 (MQ=255)
                                                        ccTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTc                                                                                                  >  2:267015/1‑134 (MQ=255)
                                                                    cGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCcctc                                                                                                                      <  3:68551/102‑1 (MQ=255)
                                                                    cGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCcctc                                                                                                                      >  4:68551/1‑102 (MQ=255)
                                                                              aGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGa                                                              >  6:57534/1‑148 (MQ=255)
                                                                               gCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGAt                                                             >  7:120238/1‑148 (MQ=255)
                                                                                 cTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCAt                                                          <  8:120238/149‑1 (MQ=255)
                                                                                   gCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCAt                                                          >  5:282370/1‑147 (MQ=255)
                                                                                      acaactacaaACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCt                                                     >  2:36968/1‑149 (MQ=255)
                                                                                        aactacaaACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCtt                                                    >  1:199572/1‑148 (MQ=255)
                                                                                              aaaCCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTa                                             >  3:116145/1‑149 (MQ=255)
                                                                                                           gcgcCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTa                                >  1:192548/1‑149 (MQ=255)
                                                                                                            cgcCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTAc                               <  5:57534/149‑1 (MQ=255)
                                                                                                             gcCACTTTGCGGCGATCCCCGCCCATCGCTTACCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTAcc                              >  4:155775/1‑149 (MQ=255)
                                                                                                             gcCACTTTGCGGCGATCCCCGCACATCGCTTACCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTAcc                              >  4:150586/1‑149 (MQ=255)
                                                                                                                  tttGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACcgcg                           >  2:237709/1‑147 (MQ=255)
                                                                                                                      cggcgATCCCCGCCCATCGCTTACCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCt                     <  3:155775/149‑1 (MQ=255)
                                                                                                                      cggcgATCCCCGCCCATCGCTTACCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCt                     <  3:150586/149‑1 (MQ=255)
                                                                                                                       ggcgATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCTg                    >  8:83844/1‑149 (MQ=255)
                                                                                                                           aTCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCACCTGCCAGGCTTg                                                                           >  4:59323/1‑90 (MQ=255)
                                                                                                                           aTCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCACCTGCCAGGCTTg                                                                           <  3:59323/90‑1 (MQ=255)
                                                                                                                             ccccGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCTGCcagg               >  6:324198/1‑148 (MQ=255)
                                                                                                                                cGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCTGCcaggcagg           <  2:192548/149‑1 (MQ=255)
                                                                                                                                cGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCTGCcaggcagg           >  8:64624/1‑149 (MQ=255)
                                                                                                                                   ccATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAg                                                                                <  2:173156/77‑1 (MQ=255)
                                                                                                                                   ccATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAg                                                                                >  1:173156/1‑77 (MQ=255)
                                                                                                                                   ccATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCa                                                  >  2:226932/1‑107 (MQ=255)
                                                                                                                                   ccATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCa                                                  <  1:226932/107‑1 (MQ=255)
                                                                                                                                         cTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCTGCCAGGCAGGCCTACgcgc  <  2:199572/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
ACTCCAGCCGCAGCGCGTACCTGGTACAAGACCTGGTACCACCCCAACAATGCCACCCTCGCCGTGGCCGGCGACATCAGCCTGCTACAACTACAAACCTTGGTGACGCGCCACTTTGCGGCGATCCCCGCCCATCGCTTGCCGATGCGGCCCTTACCCACAGGGCCCTCCTCCCAGGGCCGGCGTTTTCAGACGCTGCGCCTGCCAGGCTTGCACAACGGCGTGATCATAAGCTTCAAACTACCCAGCCAGCGTACCGCGCAGTCTGCCAGGCAGGCCTACGCGC  >  NC_002947/5836443‑5836728

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