Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,201,881 A→G 11.9% D47G (GAC→GGC)  PP_4583 → putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,201,8810AG11.9% 98.8 / 5.9 44D47G (GAC→GGC) PP_4583putative Peptidase
Reads supporting (aligned to +/- strand):  ref base A (17/21);  new base G (2/3);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

CCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGCCCAC  >  NC_002947/5201736‑5202030
                                                                                                                                                 |                                                                                                                                                     
ccTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAg                                                                                                                                                   <  1:439625/150‑1 (MQ=255)
 cTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAgg                                                                                                                                                  >  2:797429/1‑150 (MQ=255)
  ttCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGc                                                                                                                                                 >  2:73921/1‑150 (MQ=255)
  ttCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGc                                                                                                                                                 >  1:580745/1‑150 (MQ=255)
    cAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCag                                                                                                                                               <  2:582133/150‑2 (MQ=255)
                  ccccAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAgcgc                                                                                                                                 <  1:421823/150‑1 (MQ=255)
                        cccccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTc                                                                                                                           <  2:550746/150‑1 (MQ=255)
                           ccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCgaa                                                                                                                        <  2:580745/150‑1 (MQ=255)
                                  cccAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCa                                                                                                                 >  1:179533/1‑150 (MQ=255)
                                   ccAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACTCCCGAGGTACTCGCTTACCCGCGAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCaa                                                                                                                <  1:274513/150‑1 (MQ=255)
                                     aaGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCa                                                                                                                                                <  2:181594/116‑1 (MQ=255)
                                     aaGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCa                                                                                                                                                >  1:181594/1‑116 (MQ=255)
                                               ggCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCt                                                                                                    <  1:208496/150‑1 (MQ=255)
                                                gCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTg                                                                                                   <  1:73921/150‑1 (MQ=255)
                                                  cACCGACCCGTATGCCTGTCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGga                                                                                                 <  2:381350/150‑1 (MQ=255)
                                                  cACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGga                                                                                                 <  1:608387/150‑1 (MQ=255)
                                                  cACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCAGCCTGga                                                                                                 <  2:488125/150‑1 (MQ=255)
                                                    ccGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGgaga                                                                                               >  1:321703/1‑150 (MQ=255)
                                                               gCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTcc                                                                                    <  2:640695/150‑2 (MQ=255)
                                                               gCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCa                                                                                    <  2:324446/150‑1 (MQ=255)
                                                                 cTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCAtt                                                                                  >  2:787242/1‑150 (MQ=255)
                                                                   ggCTGCAACAGCGCGACACCCCCGGGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGc                                                                                <  2:588758/150‑1 (MQ=255)
                                                                     cTGCAACAGCGCGACCCCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCg                                                                              <  1:415127/150‑1 (MQ=255)
                                                                         aaCAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGgcgc                                                                          >  1:491735/1‑150 (MQ=255)
                                                                                 gACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCTTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGc                                                                  >  2:379184/1‑150 (MQ=255)
                                                                                      ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTggg                                                                    >  1:237252/1‑143 (MQ=255)
                                                                                      ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTggg                                                                    <  2:237252/143‑1 (MQ=255)
                                                                                       cccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTAc                                                            >  1:169183/1‑150 (MQ=255)
                                                                                                  cGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTacac                                                 <  1:739755/150‑1 (MQ=255)
                                                                                                       aCCTGCAAGCTGAAAACGCCTACCAGGCAGCGTAGAAGTCCCTCCCGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGCCCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTa                                            <  2:285658/150‑1 (MQ=255)
                                                                                                         cTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTacc                                          >  2:687774/1‑150 (MQ=255)
                                                                                                              aGCTGACAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAAGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCAccgc                                     <  2:360071/150‑1 (MQ=255)
                                                                                                                   aaaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcg                                >  1:632234/1‑150 (MQ=255)
                                                                                                                    aaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcga                               >  1:697786/1‑150 (MQ=255)
                                                                                                                              cAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcgacga                            <  1:71137/143‑1 (MQ=255)
                                                                                                                              cAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcgacga                            >  2:71137/1‑143 (MQ=255)
                                                                                                                                 gAAGCCTGCCTGGCCGACCAGGCAGCGTTTCGTGAGCGCCTGGTCGAAGAAAGCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGt                  >  1:447903/1‑150 (MQ=255)
                                                                                                                                 gAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGt                  >  1:687595/1‑150 (MQ=255)
                                                                                                                                 gAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGt                  <  2:179533/150‑1 (MQ=255)
                                                                                                                                      gtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGta                                                       >  2:595613/2‑108 (MQ=255)
                                                                                                                                      gtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGta                                                       <  1:595613/107‑1 (MQ=255)
                                                                                                                                       tgcctgGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTAc            <  1:3148/150‑1 (MQ=255)
                                                                                                                                             gCCGACCAGGCACCGTTGCGTGAGCGACTGTTCGAACAAATCATGGGCCGCATCCTGGAGTCCGAAAATGGGGCGGGGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGc      <  1:787242/150‑1 (MQ=255)
                                                                                                                                                gACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGAGACGAATATCCGCGTCACTATCGCTGCCCa   >  1:816422/1‑150 (MQ=255)
                                                                                                                                                 aCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGCCCAc  <  2:447903/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGCCCAC  >  NC_002947/5201736‑5202030

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