Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,410,778 G→C 11.1% H121Q (CAC→CAG PP_3894 ← replication protein O

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,410,7780GC11.1% 128.6 / 7.2 45H121Q (CAC→CAGPP_3894replication protein O
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base C (3/2);  total (25/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.19e-01

AGGGGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCTCG  >  NC_002947/4410631‑4410925
                                                                                                                                                    |                                                                                                                                                   
aGGGGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgt                                                                                                                                                    <  6:271220/149‑1 (MQ=255)
 ggggTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtg                                                                                                                                                   >  6:263870/1‑149 (MQ=255)
    gTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgtgt                                                                                                                                                >  1:19223/1‑149 (MQ=255)
        ccTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGa                                                                                                                                             <  2:223771/148‑1 (MQ=255)
             tCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCa                                                                                                                                        >  8:284496/1‑148 (MQ=255)
                          cGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCaa                                                                                                                                       >  6:108367/1‑136 (MQ=255)
                          cGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCaa                                                                                                                                       <  5:108367/136‑1 (MQ=255)
                          cGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTg                                                                                                                          >  5:256843/1‑149 (MQ=255)
                          cGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTg                                                                                                                          <  3:137442/149‑1 (MQ=255)
                                ccTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTc                                                                                                                     <  2:19223/148‑1 (MQ=255)
                                       gcCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCaa                                                                                                                                       <  2:318107/123‑1 (MQ=255)
                                       gcCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCaa                                                                                                                                       >  1:318107/1‑123 (MQ=255)
                                             cGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAg                                                                                                       >  3:115932/1‑149 (MQ=255)
                                             cGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACACTGTGTAGaaaaca                                                                                                                                        <  1:161006/116‑6 (MQ=255)
                                             cGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACACTGTGTAGaaaaca                                                                                                                                        >  2:161006/1‑111 (MQ=255)
                                                   ggatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGCTGTTTTCTACACAGTGTGTAGATTTCAACACGGATGGCTTGCTGTCGCTTCTCCACTCAGGTTCCg                                                                                                 >  8:117841/1‑149 (MQ=255)
                                                    gatgcgTCGACAACCAATTCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCgg                                                                                                <  6:256843/149‑1 (MQ=255)
                                                          tCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCgg                                                                                                <  6:32147/143‑1 (MQ=255)
                                                          tCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCgg                                                                                                >  5:32147/1‑143 (MQ=255)
                                                            gACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGaaa                                                                                         <  5:145517/148‑1 (MQ=255)
                                                                        gTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGAACAGGGGTGTTTTATACACAGTGTGTAGATTTAAACACGGTTGGCTTGATGTCGATTATCCACTCAGGTACCGGTAGTAAAGTGATTGGGTcg                                                                             >  1:86566/1‑149 (MQ=255)
                                                                                        gtatttCTGTCTTATATGTGTGTAATTTTCAACACAGGGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgc                                                             <  7:117841/148‑1 (MQ=255)
                                                                                              cTGGCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGaaaa                                                       <  5:263870/148‑1 (MQ=255)
                                                                                               tGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGaaaa                                                       <  8:276365/147‑1 (MQ=255)
                                                                                               tGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGaaaa                                                       >  7:276365/1‑147 (MQ=255)
                                                                                                   ttttATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgcg                                                            >  8:137079/1‑138 (MQ=255)
                                                                                                   ttttATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgcg                                                            <  7:137079/138‑1 (MQ=255)
                                                                                                          tgtgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCAc                                                   >  2:216793/1‑140 (MQ=255)
                                                                                                          tgtgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCAc                                                   <  1:216793/140‑1 (MQ=255)
                                                                                                               aaTTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATc                                      >  3:288873/1‑148 (MQ=255)
                                                                                                                 ttttCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgc                                                             >  6:132774/1‑123 (MQ=255)
                                                                                                                 ttttCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgc                                                             <  5:132774/123‑1 (MQ=255)
                                                                                                                       acacAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTAGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAAtt                              >  7:141899/1‑148 (MQ=255)
                                                                                                                                  aGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTgg                                                                                 >  3:74989/1‑87 (MQ=255)
                                                                                                                                  aGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTgg                                                                                 <  4:74989/87‑1 (MQ=255)
                                                                                                                                       tgtTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAATGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAAGTGCGCGGCCgagccgcg              >  8:127260/1‑144 (MQ=255)
                                                                                                                                         aaatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgcg                                                            >  6:219233/4‑101 (MQ=255)
                                                                                                                                         aaatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgcg                                                            <  5:219233/98‑1 (MQ=255)
                                                                                                                                           ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTTGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCga                    <  7:181849/139‑1 (MQ=255)
                                                                                                                                           ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCgaga                  <  5:216193/141‑1 (MQ=255)
                                                                                                                                           ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCgaga                  >  6:216193/1‑141 (MQ=255)
                                                                                                                                           ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCga                    <  7:177165/139‑1 (MQ=255)
                                                                                                                                           ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCga                    >  8:177165/1‑139 (MQ=255)
                                                                                                                                           ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCga                    >  8:181849/1‑139 (MQ=255)
                                                                                                                                           atCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCTAGACACGTTgg          >  3:240840/2‑149 (MQ=255)
                                                                                                                                                 acaGTGTGTGGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCt    >  3:264847/1‑149 (MQ=255)
                                                                                                                                                 acaGTGTGTAGATTTCAACACGGTTGGATTGCTTTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCt    >  4:48542/1‑149 (MQ=255)
                                                                                                                                                   aGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGAAGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCTCg  >  1:137316/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGGGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCTCG  >  NC_002947/4410631‑4410925

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: