Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,850,090 A→C 16.0% Q604P (CAG→CCG)  PP_3397 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,850,0900AC16.0% 51.3 / 5.0 25Q604P (CAG→CCG) PP_3397hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base C (2/2);  total (13/12)
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.41e-01

CATCAACGTCGCACCTAACCAGGCACCGACTGCCAACCCGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCACCTACCAGGCACGTGACTCACTGGGGGGTCTGTCCACCCCAAGTACGGTC  >  NC_002947/3849950‑3850235
                                                                                                                                            |                                                                                                                                                 
cATCAACGTCGCACCTAACCAGGCACCGACTGCCAACCCGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAgccggccgc                                                                                                                                          >  2:85742/1‑150 (MQ=255)
  tCAACGTCGCACCTAACCAGGCACCGACTGCCAACCCGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAgccggccgccg                                                                                                                                        >  2:601397/1‑150 (MQ=255)
        tCGCACCTAACCAGGCACCGACTGCCAACCCGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTg                                                                                                                                  >  1:262087/1‑150 (MQ=255)
        tCGCACCTAACCAGGCACCGACTGCCAACCCGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTg                                                                                                                                  <  2:262087/150‑1 (MQ=255)
          gCACCTAACCAGGCACCGACTGCCAACCCGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGccggccggcggccggctg                                                                                                                                <  2:596354/150‑5 (MQ=255)
                                     ccGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCAGCCAGCCGGCCGCCGGCCGTg                                                                                                                                  >  1:535804/1‑121 (MQ=255)
                                     ccGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTg                                                                                                                                  <  2:535804/121‑1 (MQ=255)
                                       ggAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGCCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGAc                                                                                                   <  1:95504/150‑1 (MQ=255)
                                                 gCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCa                                                                                         <  1:85742/150‑1 (MQ=255)
                                                               gCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAAc                                                                           <  1:36425/150‑1 (MQ=255)
                                                                  aTTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACcac                                                                        >  2:94754/1‑150 (MQ=255)
                                                                        cATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAAc                                                                           >  2:97177/1‑141 (MQ=255)
                                                                        cATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAAc                                                                           <  1:97177/141‑1 (MQ=255)
                                                                         aTCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCCGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTc                                                                 >  2:64375/1‑150 (MQ=255)
                                                                         aTCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCCGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTc                                                                 >  2:64084/1‑150 (MQ=255)
                                                                           cAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGc                                                                                    <  1:243876/129‑1 (MQ=255)
                                                                           cAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGc                                                                                    >  2:243876/1‑129 (MQ=255)
                                                                                      ccAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGTCGTGGCACGGTCAGCACCGACGACAACCCAGCTGGGTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCa                                                    >  2:476421/1‑150 (MQ=255)
                                                                                       cAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCCGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCAc                                                   <  1:64084/150‑1 (MQ=255)
                                                                                       cAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCCGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCAc                                                   <  1:64375/150‑1 (MQ=255)
                                                                                       cAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCAc                                                   <  1:94754/150‑1 (MQ=255)
                                                                                                                               gCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCACCTACCAGGCACGTGACTCACTGGGGGGTCTGTCCACCCCa           >  2:92902/1‑150 (MQ=255)
                                                                                                                                 cAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCACCTACCAGGCACGTGACTCACTGGGGGGTCTGTCCACCCCAAg         >  2:336377/1‑150 (MQ=255)
                                                                                                                                        aCCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCACCTACCAGGCACGTGACTCACTGGGGGGTCTGTCCCCCCCAAGTACGGTc  >  2:297009/1‑150 (MQ=255)
                                                                                                                                        aCCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCACCTACCAGGCACGTGACTCACTGGGGGGTCTGTCCACCCCAAGTACGGTc  <  1:601397/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CATCAACGTCGCACCTAACCAGGCACCGACTGCCAACCCGGAAACCGTCGCCACCTTGGGCGTGCCATTGACCATCAACGTGCTGGCCAACGACACCGATCCAGAAGGCAACGTACCGCTGAGCGTGGCCAACATCACCCAGCCGGCCGCCGGCCGTGGCACGGTCAGCACCGACGGCACAACCGTGACCTATACCCCACCGGCCACTGTAACCACAGCCTTCACCACCACCTTCACCTACCAGGCACGTGACTCACTGGGGGGTCTGTCCACCCCAAGTACGGTC  >  NC_002947/3849950‑3850235

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