Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,042,173 G→C 9.3% E26Q (GAG→CAG) ‡ PP_5537 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,042,1730GC 9.3% 129.1 / 3.6 43E26Q (GAG→CAG) ‡PP_5537membrane protein
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base C (2/2);  total (23/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 = 8.42e-01

GACCTCACGGGTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGTGC  >  NC_002947/3042027‑3042320
                                                                                                                                                  |                                                                                                                                                   
gACCTCACGGGTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAg                                                                                                                                                   >  8:137546/1‑149 (MQ=255)
         ggTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCgctg                                                                                                                                          <  2:237893/149‑1 (MQ=255)
                 aGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTACGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATcagcat                                                                                                                                         <  1:305738/142‑6 (MQ=255)
                 aGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTACGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATcagcat                                                                                                                                         >  2:305738/1‑137 (MQ=255)
                  gTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTc                                                                                                                                 <  6:67898/149‑1 (MQ=255)
                      aagcaaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTc                                                                                                                                 <  4:224865/145‑1 (MQ=255)
                      aagcaaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTc                                                                                                                                 >  3:224865/1‑145 (MQ=255)
                        gcaaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTg                                                                                                                           <  5:312991/149‑1 (MQ=255)
                         caaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTgg                                                                                                                          <  8:165180/149‑1 (MQ=255)
                              cTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCAct                                                                                                                      >  3:96743/1‑148 (MQ=255)
                               tATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGATTGCGTATCTGCGCGCTGGCAATACTGGTGATGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCAGAGCAATCGCTGCTGGCATTCGGGCTGGCActgc                                                                                                                    >  1:149869/1‑149 (MQ=255)
                                aTCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCActgct                                                                                                                   >  2:322283/1‑149 (MQ=255)
                                 tCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCActgc                                                                                                                    <  6:117711/147‑1 (MQ=255)
                                   actactGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGgc                                                                                                                >  5:146430/1‑149 (MQ=255)
                                       ctGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCa                                                                                                             <  5:43644/148‑1 (MQ=255)
                                          tgCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGc                                                                                                          <  6:146430/148‑1 (MQ=255)
                                               gCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGacac                                                                                                    <  3:139904/149‑1 (MQ=255)
                                                 aaaCCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACgg                                                                                                  <  1:322283/149‑1 (MQ=255)
                                                     ccTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATca                                                                                                                                             >  8:175098/1‑101 (MQ=255)
                                                     ccTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATca                                                                                                                                             <  7:175098/102‑2 (MQ=255)
                                                       tATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAgg                                                                                            >  3:366477/1‑149 (MQ=255)
                                                             ggAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCa                                                                                      >  6:20244/1‑149 (MQ=255)
                                                                gagCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATc                                                                                    <  4:96743/148‑1 (MQ=255)
                                                                  gCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGAc                                                                                 <  7:133490/149‑1 (MQ=255)
                                                                   cTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGAcc                                                                                <  7:137546/149‑1 (MQ=255)
                                                                         gggCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTAcc                                                                          >  5:212422/1‑149 (MQ=255)
                                                                                   aTCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCgg                                                                >  3:42988/1‑149 (MQ=255)
                                                                                       gcgcgcTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTg                                                            >  7:241796/1‑149 (MQ=255)
                                                                                           gcTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACCCGGCCCAGGTGGTCATAGACCTGTACCTGATGGCCGGGCTGGCCg                                                        >  8:333192/1‑149 (MQ=255)
                                                                                                 aaTACTTTTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGATGGCGCAGGCGGACACGGCCAAGCTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCt                                                  >  4:243801/1‑149 (MQ=255)
                                                                                                      tGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACg                                                                                                   <  7:89292/95‑1 (MQ=255)
                                                                                                      tGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACg                                                                                                   >  8:89292/1‑95 (MQ=255)
                                                                                                        gTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGt                                           <  4:366477/149‑1 (MQ=255)
                                                                                                                  cGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACaa                                  <  5:20244/148‑1 (MQ=255)
                                                                                                                                       tGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGCCCCAGGTGGTCATCGACCTGTACCTGATCGCCGGGCTGGCCGGTTGCTGGATGGTCAGAGACAAACGCATGCGCGGCAGAGACgg             >  6:272443/1‑148 (MQ=255)
                                                                                                                                           gATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGCACCGCGATcc        >  2:162272/1‑148 (MQ=255)
                                                                                                                                           gATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATcg        >  1:169297/1‑149 (MQ=255)
                                                                                                                                             ttGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATcgcg      <  4:42988/149‑1 (MQ=255)
                                                                                                                                              tGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCgg     >  3:365925/1‑149 (MQ=255)
                                                                                                                                                 cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGc                                                   >  4:262165/1‑100 (MQ=255)
                                                                                                                                                 cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGc                                                   <  3:262165/100‑1 (MQ=255)
                                                                                                                                                 cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGtgc  >  4:6622/1‑149 (MQ=255)
                                                                                                                                                 cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGtgc  >  3:338577/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACCTCACGGGTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGTGC  >  NC_002947/3042027‑3042320

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