Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,908,672 G→A 13.8% A193T (GCC→ACC)  PP_3450 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,908,6720GA13.8% 60.2 / 3.7 29A193T (GCC→ACC) PP_3450TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (12/13);  new base A (2/2);  total (14/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 = 4.94e-01

TCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGACCGGCGACCT  >  NC_002947/3908524‑3908820
                                                                                                                                                    |                                                                                                                                                    
tCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGc                                                                                                                                                     >  2:477976/1‑150 (MQ=255)
   aCTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCg                                                                                                                                                  >  1:121194/1‑150 (MQ=255)
       cagcaACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGCCACATCAAGACCCGGACCTGGCCCGGGaa                                                                                                                                              <  1:286231/150‑1 (MQ=255)
          caACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCCGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGCCCTGGCCCGGGAACAg                                                                                                                                           >  1:451681/1‑150 (MQ=255)
             cccGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTc                                                                                                                                        >  2:369274/1‑150 (MQ=255)
               cGAGCAGCCCCAGCAACTTAGGCACCTGGGCCGCGAAGGCCACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTAACATCAAGACCCGGACCTGGCACGGGAACAGTTCag                                                                                                                                      <  2:462510/150‑2 (MQ=255)
                    agcGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCt                                                                                                                                 <  1:369274/150‑1 (MQ=255)
                         cAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGcccc                                                                                                                            >  1:217002/1‑150 (MQ=255)
                                                  aaGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTAAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATg                                                                                                   <  1:96365/150‑1 (MQ=255)
                                                   aGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGAGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGc                                                                                                  <  1:342744/150‑1 (MQ=255)
                                                      ccacgaTGCAGCGCTGGTCCGGCGCCTCAGCAACGACCGGCAGCGCCCCTGCCTTGAACCGGTCCTCTGGCTGTCACATCAAGCCCCGGACCTGGCCCGGGAACCGGTCAAGCCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGt                                                                                               <  2:79278/148‑1 (MQ=255)
                                                                    tGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGcc                                                                                 <  1:477976/150‑1 (MQ=255)
                                                                              ccTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAg                                                                                                                                           >  1:383132/1‑82 (MQ=255)
                                                                              ccTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAg                                                                                                                                           <  2:383132/82‑1 (MQ=255)
                                                                                                      cTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGACCCGGGAACAGTTCAAGGCCTGCCCCGCCaa                                                                                                                       >  1:3694/1‑78 (MQ=255)
                                                                                                      cTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGACCCGGGAACAGTTCAAGGCCTGCCCCGCCaa                                                                                                                       >  1:3565/1‑78 (MQ=255)
                                                                                                      cTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGACCCGGGAACAGTTCAAGGCCTGCCCCGCCaa                                                                                                                       <  2:3565/78‑1 (MQ=255)
                                                                                                      cTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGACCCGGGAACAGTTCAAGGCCTGCCCCGCCaa                                                                                                                       <  2:3694/78‑1 (MQ=255)
                                                                                                                     cGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATg                                <  2:121194/150‑1 (MQ=255)
                                                                                                                         tGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGaa                            <  1:481358/150‑1 (MQ=255)
                                                                                                                         tGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGaa                            <  2:404996/150‑1 (MQ=255)
                                                                                                                         tGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCCGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGaa                            >  1:189673/1‑150 (MQ=255)
                                                                                                                                   aTCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATc                  >  1:580326/1‑150 (MQ=255)
                                                                                                                                   aTCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATc                  >  1:446084/1‑150 (MQ=255)
                                                                                                                                      aaGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGc               >  1:159965/1‑150 (MQ=255)
                                                                                                                                       aGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGc               <  1:505978/149‑1 (MQ=255)
                                                                                                                                       aGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGc               >  2:505978/1‑149 (MQ=255)
                                                                                                                                           ccGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCa                                                              >  1:610976/1‑98 (MQ=255)
                                                                                                                                           ccGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCa                                                              <  2:610976/98‑1 (MQ=255)
                                                                                                                                                   ggCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGACCGGCGACCt  >  2:523281/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGACCGGCGACCT  >  NC_002947/3908524‑3908820

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