Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,112,775 C→A 18.2% P518P (CCG→CCT) ‡ PP_1882 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,112,7750CA18.2% 50.2 / 4.2 22P518P (CCG→CCT) ‡PP_1882membrane protein
Reads supporting (aligned to +/- strand):  ref base C (11/7);  new base A (2/2);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

AGGGTCTAGCAGTGAATTGAACTGGTACTGGGTGGTGGTTTCCAGCTTGTTGCGTTTCACCACCTCCTGATAGCGACGGCCGTGCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCTGAGTGAGAACAGGGGCTTGCCCCCGGCCGCTGCTCGCCGGG  >  NC_002947/2112631‑2112919
                                                                                                                                                |                                                                                                                                                
agGGTCTAGCAGTGAATTGAACTGGTACTGGGTGGTGGTTTCCAGCTTGTTGCGTTTCACCACCTCCTGATAGCGACGGCCGTGCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATc                                                                                                                                             <  1:377920/150‑1 (MQ=255)
                 tGAACTGGTACTGGGTGGTGGTTTCCAGCTTGTTGCGTTTCACCACCTCCTGATAGCGACGGCCGTGCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAggg                                                                                                                            >  1:63316/1‑150 (MQ=255)
                                tggtggtTTCCAGCTTGTTGCGTTTCACCACCTCCTGATAGCGACGGCCGTGCCGCAAGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGa                                                                                                             >  1:198983/1‑150 (MQ=255)
                                         ccAGCTTGTTGCGTTTCACCACCTCCTGATAGCGACGGCCGTGCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCg                                                                                                    <  2:168421/150‑1 (MQ=255)
                                                              cctcctGATAGCGACGGCCGTGCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTtgc                                                                               >  2:449336/1‑150 (MQ=255)
                                                                                   gCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCg                                                                                                    >  2:175759/1‑108 (MQ=255)
                                                                                   gCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCg                                                                                                    <  1:175759/108‑1 (MQ=255)
                                                                                      gCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGt                                                       >  2:362387/1‑150 (MQ=255)
                                                                                          gCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCg                                                   >  2:129091/1‑150 (MQ=255)
                                                                                          gCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCg                                                   >  2:363445/1‑150 (MQ=255)
                                                                                          gCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCg                                                   >  1:223792/1‑150 (MQ=255)
                                                                                                      cgcATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCtgagt                                       <  1:362387/150‑1 (MQ=255)
                                                                                                         aTCGCGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGTTTCAGGGTTTCACTGTGCTCGACAACCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCtgagtgag                                    >  2:438157/1‑150 (MQ=255)
                                                                                                              ggTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCTGAGTGAGAACAg                               <  1:363445/150‑1 (MQ=255)
                                                                                                               gTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCTGAGTGAGAACAgg                              <  2:63316/150‑1 (MQ=255)
                                                                                                                     gTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCTGAGTGAGAACAGGGGCTTg                        <  2:198983/150‑1 (MQ=255)
                                                                                                                                  cctgacgcATGCGCACGATCCGGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGgc                                                                    >  1:386946‑M1/9‑93 (MQ=255)
                                                                                                                                  cctgacgcATGCGCACGATCCGGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGgc                                                                    <  2:386946‑M1/85‑1 (MQ=255)
                                                                                                                                   ctgacgcATGCGCACGATCCGGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGtt                                                                                 <  2:205950‑M1/72‑1 (MQ=255)
                                                                                                                                   ctgacgcATGCGCACGATCCGGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGtt                                                                                 >  1:205950‑M1/8‑79 (MQ=255)
                                                                                                                                     aGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCTGAGTGGGAACAGGGGCTTGCCCCCGGCCGCTGCTc        >  2:101957/1‑150 (MQ=255)
                                                                                                                                           TGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCTGAGTGAGAACAGGGGCTTGCCCCCGGCCGCTGCTCGCCggg  >  1:387345‑M1/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGGGTCTAGCAGTGAATTGAACTGGTACTGGGTGGTGGTTTCCAGCTTGTTGCGTTTCACCACCTCCTGATAGCGACGGCCGTGCCGCAGGCTGGACTTGGGCGCATCGAGGTAATAGTAAATTTTCTCTTCCAGGGGATGCGCCGGATCGTGCGCATGCGTCAGGGTTTCACTGTGCTCGACAATCATCGGGTTGAGGGTGATGGCGTTGCCTGCCAGTGGCGCCAGCGCCTTGTAGCGGTAGTGCTGAGTGAGAACAGGGGCTTGCCCCCGGCCGCTGCTCGCCGGG  >  NC_002947/2112631‑2112919

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