Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,618,250 A→G 9.7% H2R (CAC→CGC)  PP_5652 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,618,2500AG 9.7% 105.6 / 3.2 41H2R (CAC→CGC) PP_5652hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base G (2/2);  total (23/18)
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.95e-01

GAGGTCATGCAACTTTCGGAGCGAATGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCTATAAACTTTGCGGCGGCGCCCAGTTG  >  NC_002947/4618112‑4618398
                                                                                                                                          |                                                                                                                                                    
gagGTCATGCAACTTTCGGAGCGAATGTTTGACGTCGCAGGCGTTCCAGCCACCAGCCGAACCTCCTACTACAGTGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATaaaa                                                                                                                                            <  5:311810/149‑1 (MQ=255)
       tGCAACTTTCGGAGCGAATGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATAc                                                                                                                                     >  5:270460/1‑149 (MQ=255)
       tGCAACTTTCGGAGCGAATGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATAc                                                                                                                                     >  7:145072/1‑149 (MQ=255)
        gCAACTTTCGGAGCGAATGTTTGACGTCGCAGGAGTTCCAGCCACGAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACt                                                                                                                                    >  1:254514/1‑149 (MQ=255)
           aCTTTCGGAGCGAATGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTga                                                                                                                                 <  4:167644/149‑1 (MQ=255)
             tttCGGAGCGAATGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTgag                                                                                                                                >  5:364698/1‑148 (MQ=255)
                       aaTGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCGCACCATAAAACACATACTTGAGAAATTCAGCaa                                                                                                                     >  2:91980/1‑149 (MQ=255)
                       aaTGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCGCACCATAAAACACATACTTGAGAAATTCAGCaa                                                                                                                     >  2:87211/1‑149 (MQ=255)
                       aaTGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAacac                                                                                                                                         <  3:100582/129‑1 (MQ=255)
                       aaTGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAacac                                                                                                                                         >  4:100582/1‑129 (MQ=255)
                           ttagACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCAGACCATAAAACACATACTTGAGAAATTCAGCAATGAt                                                                                                                 >  8:119079/4‑149 (MQ=255)
                                cGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCga                                                                                                             <  2:254514/148‑1 (MQ=255)
                                       ggAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGAc                                                                                                     <  6:270460/149‑1 (MQ=255)
                                               aGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATc                                                                                                                >  3:137617/1‑130 (MQ=255)
                                               aGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATc                                                                                                                <  4:137617/130‑1 (MQ=255)
                                                gCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCt                                                                                            <  2:78608/149‑1 (MQ=255)
                                                   accaGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTcc                                                                                          >  7:231028/1‑148 (MQ=255)
                                                   accaGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTc                                                                                           >  5:21447/1‑147 (MQ=255)
                                                       gCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCt                                                                                            >  2:79219/1‑142 (MQ=255)
                                                       gCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCt                                                                                            <  1:79219/142‑1 (MQ=255)
                                                        ccGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTc                                                                                    <  7:43977/149‑1 (MQ=255)
                                                          gAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTcac                                                                                  <  6:21447/149‑1 (MQ=255)
                                                              ctcctACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTACACATGCACACCATAAAAGTCATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTcaccac                                                                               >  6:166399/1‑148 (MQ=255)
                                                                 ctactaCAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACcactacc                                                                           <  6:364698/149‑1 (MQ=255)
                                                                     tacAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACcactacccac                                                                        >  1:140139/1‑148 (MQ=255)
                                                                       cAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACcactacccacac                                                                      >  6:270585/1‑148 (MQ=255)
                                                                         gggCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCGCACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACcactacccacactac                                                                   <  1:91980/149‑1 (MQ=255)
                                                                         gggCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCGCACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACcactacccacactac                                                                   <  1:87211/149‑1 (MQ=255)
                                                                          ggCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACcactacccacactacc                                                                  <  2:308272/149‑1 (MQ=255)
                                                                           gCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACcactacccacactacc                                                                  >  3:306683/1‑148 (MQ=255)
                                                                              ttGAATAGATATTTACACACAGGCACTTCCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAAc                                                              >  4:102666/1‑149 (MQ=255)
                                                                                     gATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCtt                                                        <  4:306683/148‑1 (MQ=255)
                                                                                              acacaGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTc                                                           >  7:257218/1‑136 (MQ=255)
                                                                                              acacaGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTc                                                           <  8:257218/136‑1 (MQ=255)
                                                                                                                 cTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCt                            <  7:119079/148‑1 (MQ=255)
                                                                                                                  tAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCta                           >  8:289084/1‑148 (MQ=255)
                                                                                                                  tAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGACAATTGACGTTTTCCATCCAGACTACCAACAACCCACACAACCCAACTTCATTAGCTATCCTCCAGAACTCGATGAATACta                           >  3:73124/1‑148 (MQ=14)
                                                                                                                    aattGAGGCAAGCTAAACATGCACACCATAAAACACATACTTTAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCtat                          <  8:145072/145‑1 (MQ=255)
                                                                                                                                cTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCTATAAACTTTgcggcgg            >  1:238954/1‑149 (MQ=255)
                                                                                                                                    aCATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCTATAAACTTTGCGGCGGCGcc        >  3:112061/1‑149 (MQ=255)
                                                                                                                                      aTGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCTATAAACTTTGCGGCGGCGcc        >  8:188683/1‑147 (MQ=255)
                                                                                                                                          acacCATAAAACACATACTTGAGATATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCCCACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGATTTCTATAAACTTTGCGTCGGGGCCGAggtg  >  8:141246/1‑149 (MQ=14)
                                                                                                                                          acacCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCTATAAACTTTGCGGCGGCGCCCAgttg  >  8:41952/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GAGGTCATGCAACTTTCGGAGCGAATGTTTGACGTCGCAGGAGTTCCAGCCACCAGCCGAACCTCCTACTACAGGGCATTGAATAGATATTTACACACAGGCACTTTCTGATCCTAATTTGAGGCAAGCTAAACATGCACACCATAAAACACATACTTGAGAAATTCAGCAATGATCCGAGAATTGACGTTTTCCCTCCAGACTCACCACTACCCACACTACCCAACTTCCTTAGCTATCCTCCAGACCTCGCTGAATTCTATAAACTTTGCGGCGGCGCCCAGTTG  >  NC_002947/4618112‑4618398

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