Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,979,338 C→A 16.1% G316V (GGG→GTG)  PP_2605 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,979,3380CA16.1% 85.3 / 5.5 35G316V (GGG→GTG) PP_2605hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/16);  new base A (4/2);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTAGCAATCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGATCATGTGG  >  NC_002947/2979193‑2979486
                                                                                                                                                 |                                                                                                                                                    
gTGGTAGCAATCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGt                                                                                                                                                   <  3:211715/149‑1 (MQ=255)
    tAGCAATCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGgc                                                                                                                                               <  2:77908/149‑1 (MQ=255)
        aaTCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCt                                                                                                                                           >  2:299682/1‑149 (MQ=255)
          tCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGc                                                                                                                                         >  5:304550/1‑149 (MQ=255)
          tCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCACgggg                                                                                                                                                 >  1:112975/1‑141 (MQ=255)
          tCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCACgggg                                                                                                                                                 <  2:112975/141‑1 (MQ=255)
                     gACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGAcacc                                                                                                                               <  2:201757/148‑1 (MQ=255)
                            cgcCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCACGGGGGCGCCTgg                                                                                                                                         >  5:115911/1‑130 (MQ=255)
                            cgcCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCACGGGGGCGCCTgg                                                                                                                                         <  6:115911/131‑2 (MQ=255)
                               cGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGAcacccgcccacccgc                                                                                                                    <  6:304550/149‑1 (MQ=255)
                                  ccATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTg                                                                                                                 >  1:264820/1‑149 (MQ=255)
                                    aTCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCTCACCCGCCTGtt                                                                                                               <  8:206100/149‑1 (MQ=255)
                                     tCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTc                                                                                                              <  1:299682/149‑1 (MQ=255)
                                     tCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTc                                                                                                              >  8:82450/1‑149 (MQ=255)
                                                     tGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCa                                                                                              <  7:82450/149‑1 (MQ=255)
                                                         aCCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCAc                                                                                          <  2:264820/149‑1 (MQ=255)
                                                              gCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTa                                                                                     <  7:188885/149‑1 (MQ=255)
                                                                    cACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCt                                                                                                                  >  3:245701/1‑114 (MQ=255)
                                                                    cACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCt                                                                                                                  <  4:245701/114‑1 (MQ=255)
                                                                           gATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGgcg                                                                        >  2:247046/1‑149 (MQ=255)
                                                                                ggTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGATCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCg                                                                   <  1:247046/149‑1 (MQ=255)
                                                                                ggTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCACgg                                                                                                                                                   >  8:290807/1‑68 (MQ=255)
                                                                                           tggccgaaatggccATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCAc                                                                                          >  7:9667/1‑115 (MQ=255)
                                                                                           tggccgaaatggccATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCAc                                                                                          <  8:9667/115‑1 (MQ=255)
                                                                                                                        cTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCAt                                <  5:146133/144‑1 (MQ=255)
                                                                                                                        cTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCAt                                >  6:146133/1‑144 (MQ=255)
                                                                                                                         tGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGc                                  <  2:258762/141‑1 (MQ=255)
                                                                                                                         tGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGc                                  >  1:258762/1‑141 (MQ=255)
                                                                                                                         tGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTg                          >  2:62099/1‑149 (MQ=255)
                                                                                                                                 cGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGAt                    >  8:37244/1‑147 (MQ=255)
                                                                                                                                 cGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGAt                    <  7:37244/147‑1 (MQ=255)
                                                                                                                                   gcgcCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGc                >  1:165983/1‑149 (MQ=255)
                                                                                                                                      cccAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTgg             >  4:90356/1‑149 (MQ=255)
                                                                                                                                          ggCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGAt         <  2:343818/149‑1 (MQ=255)
                                                                                                                                             gcCCACGGGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGATCAt      >  3:213914/1‑149 (MQ=255)
                                                                                                                                              cccccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGa          >  4:333694/1‑144 (MQ=255)
                                                                                                                                              cccccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGa          >  4:181788/1‑144 (MQ=255)
                                                                                                                                              cccccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGa          <  3:333694/144‑1 (MQ=255)
                                                                                                                                              cccccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGa          <  3:181788/144‑1 (MQ=255)
                                                                                                                                               ccccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCga                       >  7:237810/1‑130 (MQ=255)
                                                                                                                                               ccccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCga                       <  8:237810/130‑1 (MQ=255)
                                                                                                                                                 ccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGATCATGTgg  >  6:59417/1‑149 (MQ=255)
                                                                                                                                                 ccGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGATCATGTgg  >  1:181522/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTGGTAGCAATCGGCTTCGTTGACGAAGCGCCGACCATCGCGGCGCACGGCGATGAAACCGGGCTTGGCACGGTCGATCAGGTGTGGGAAATGGCCGAAATGGCCATTGCGGCGCGGTACCTGCGATACCGGCGCCCAGGCGCCCCCGTGGGCGCCTGCCGCAGACACCCGCCCACCCGCCTGTTCGCCCAGGCGCAAACCATCACCGCTATTCTCCCTGGGCGCAGCGGAGTAATGCTGGTCACCGTCCGGCGCGTGCGGCATCAGTTGCGAGATACGCTGGCGATCATGTGG  >  NC_002947/2979193‑2979486

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