Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,973,809 A→G 9.9% L22P (CTG→CCG)  PP_2600 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,973,8090AG 9.9% 102.6 / 2.5 41L22P (CTG→CCG) PP_2600hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base G (2/2);  total (21/20)
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.78e-01

CCTCGATGCCGGCCCTGCGCGCCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCGGGTGGTCCTGCAAGGTCAGGTTG  >  NC_002947/2973666‑2973950
                                                                                                                                               |                                                                                                                                             
ccTCGATGCCGGCCCTGCGCGCCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCcg                                                                                                                                          >  3:222710/1‑149 (MQ=255)
 cTCGATGCCGGCCCTGCGCGCCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCGGCCcgg                                                                                                                                         >  5:255415/1‑149 (MQ=255)
          ggCCCTGCGCGCCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGcc                                                                                                                                  >  4:77047/1‑147 (MQ=255)
          ggCCCTGCGCGCCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGcc                                                                                                                                  <  3:77047/147‑1 (MQ=255)
              cTGCGCGGCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAAcc                                                                                                                            <  4:6218/149‑1 (MQ=255)
                 cgcgcCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCGGCCCGGCGGCGCCAGAACCTAt                                                                                                                         <  6:255415/149‑1 (MQ=255)
                 cgcgcCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTAt                                                                                                                         <  6:332602/149‑1 (MQ=255)
                 cgcgcCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTAt                                                                                                                         <  8:372602/149‑1 (MQ=255)
                      cAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGAc                                                                                                                      >  5:2913/1‑147 (MQ=255)
                       aGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCCGCCCGGCGGCGCCAGAACCTATGACCAg                                                                                                                   >  7:168099/1‑149 (MQ=255)
                       aGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAg                                                                                                                   >  1:1027/1‑149 (MQ=255)
                                 gTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCAt                                                                                                         >  4:176482/1‑149 (MQ=255)
                                            gCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCGGGCTGgcggcggcc                                                                                                                                 >  2:102858/1‑111 (MQ=255)
                                            gCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCGGGCTGgcggcggcc                                                                                                                                 <  1:102858/114‑4 (MQ=255)
                                            gCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCgag                                                                                              <  7:51822/149‑1 (MQ=255)
                                              ccAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGagag                                                                                            >  4:285724/1‑149 (MQ=255)
                                                aGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAgagga                                                                                          <  2:1027/149‑1 (MQ=255)
                                                    gccgGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGCCCAGTACATCGCATTCGGTTTCGAGAGGAGCg                                                                                       >  5:170127/1‑148 (MQ=255)
                                                    gccgGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTa                                                                                                                 <  1:70162/122‑1 (MQ=255)
                                                    gccgGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTa                                                                                                                 >  2:70162/1‑122 (MQ=255)
                                                           cAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCAtt                                                                                                        <  3:310529/124‑1 (MQ=255)
                                                           cAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCAtt                                                                                                        >  4:310529/1‑124 (MQ=255)
                                                           cAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCAt                                                                               >  1:382516/1‑149 (MQ=255)
                                                               cGGTGGCACCGCCAAACACCGGTTCTTTC‑CCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAg                                                                                                                   <  6:124373/108‑1 (MQ=255)
                                                                                     ttctttctCCGCCAAGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGTAGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCTAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGa                                                     <  8:402471/149‑1 (MQ=255)
                                                                                          tctcCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAAt                                                >  1:395631/1‑149 (MQ=255)
                                                                                              cGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTc                                            >  3:182668/1‑149 (MQ=255)
                                                                                              cGCCACGATCACCTTCTGCCCATGCCAGGCAGCCATGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTc                                            <  8:266506/149‑1 (MQ=255)
                                                                                                     aTCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTc                                     <  3:176482/149‑1 (MQ=255)
                                                                                                            tCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTCCGGGATTCCTGTGCTTTTATTAGAACAAAGTTCCACATTCTTATcag                              >  7:3498/1‑149 (MQ=255)
                                                                                                               gCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATcagcag                           >  6:263726/1‑149 (MQ=255)
                                                                                                                cccATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCcaca                                        <  2:10950/135‑1 (MQ=255)
                                                                                                                cccATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCcaca                                        >  1:10950/1‑135 (MQ=255)
                                                                                                                 ccATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCg                         <  3:285724/149‑1 (MQ=255)
                                                                                                                   aTGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCggg                       <  4:222710/149‑1 (MQ=255)
                                                                                                                   aTGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCggg                       <  3:163869/149‑1 (MQ=255)
                                                                                                                       cAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCGggtggt                   >  2:189600/1‑149 (MQ=255)
                                                                                                                          gCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCGGGTGGTCCt                <  2:395631/149‑1 (MQ=255)
                                                                                                                             cccGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGAACAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCGTATCAGCAGCGGGTGGTCCTGc              >  5:339847/2‑148 (MQ=255)
                                                                                                                                gTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCATGTGCTTTTATTAGAACAATGTTCCACATTATTATCAGCATCCGGTGGTCCTGCaaag          >  4:225910/1‑147 (MQ=255)
                                                                                                                                        cgccgccAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCGGGTGGTCCTGCAAGGTCAGGTTg  <  4:338886/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CCTCGATGCCGGCCCTGCGCGCCAGCGGGTTGCGTGGCACCCAGGCCCAGCCGCCGGACCAGGCGGTGGCACCGCCAAACACCGGTTCTTTCTCCGCCACGATCACCTTCTGCCCATGCCAGGCAGCCGTGACGGCCGCCGCCAGCCCGGCGGCGCCAGAACCTATGACCAGTACATCGCATTCGGTTTCGAGAGGAGCGGAAGGCATTGCGGGTTTCCTGTGCTTTTATTAGAACAATGTTCCACATTCTTATCAGCAGCGGGTGGTCCTGCAAGGTCAGGTTG  >  NC_002947/2973666‑2973950

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: ‑