Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,410,822 G→C 21.0% 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,8220GC21.0% 36.4 / 5.5 19H121Q (CAC→CAGPP_3894replication protein O
Reads supporting (aligned to +/- strand):  ref base G (9/6);  new base C (2/2);  total (11/8)
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.72e-01

GGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCG  >  NC_002947/4410678‑4410965
                                                                                                                                                |                                                                                                                                               
ggTTGCCTTGTCGTTGCCGTGACCTCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgtgt                                                                                                                                            <  2:191406/150‑1 (MQ=255)
      cTTGTCGTTGCCATGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGAttt                                                                                                                                      >  1:425627/1‑150 (MQ=255)
                        gCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCt                                                                                                                    <  2:366921/150‑1 (MQ=255)
                                        gTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTCTTTATGTTTCTCATTCAGAAGGACGGAAGTGTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTTCTGTCGATTCTCCACTCa                                                                                                    <  1:227398/150‑1 (MQ=255)
                                                       tCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGaaa                                                                                     >  1:462476/1‑150 (MQ=255)
                                                       tCGACAACCAATGCGTTAGCATTGGAATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATGTCAACACGGTTGGCTTGCTGTCGATTCTCCACGCAGGTACCGGTAGGaaa                                                                                     >  1:153891/1‑150 (MQ=255)
                                                              ccAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTg                                                                              >  2:41512/1‑150 (MQ=255)
                                                                   gCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgt                                                                                                                                              <  2:149467/81‑1 (MQ=255)
                                                                   gCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgt                                                                                                                                              >  1:149467/1‑81 (MQ=255)
                                                                             tGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTgccgc                                                               <  1:41512/150‑1 (MQ=255)
                                                                                  tttgtatttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcg                                                          >  1:425379/1‑150 (MQ=255)
                                                                                   ttgtatttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgc                                                         <  2:462476/150‑1 (MQ=255)
                                                                                                       tatgtAATTTTCAACACAGTGGAAGGGGTGTTTTCTAAACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTACACTCAGATAACGGTAGGAAAGTGATTGGCTCGCGACTGCCGCCATAGCGGAAAAGCACGCGCTGACgg                                     <  2:33903/148‑1 (MQ=255)
                                                                                                         tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGAt                                   >  1:5262/1‑150 (MQ=255)
                                                                                                         tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGAt                                   >  1:591080/1‑150 (MQ=255)
                                                                                                                    acacAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCGGCACGGGTGGCTTGCTTAAAAATCTAAACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGa                        >  2:235080/1‑150 (MQ=255)
                                                                                                                          tGGCAGGGGTTTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGccc                  <  2:434074/150‑1 (MQ=255)
                                                                                                                                      aatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGAtt                                                                               <  2:208586/75‑1 (MQ=255)
                                                                                                                                      aatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGAtt                                                                               <  2:208625/75‑1 (MQ=255)
                                                                                                                                      aatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGAtt                                                                               >  1:208625/3‑77 (MQ=255)
                                                                                                                                      aatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGAtt                                                                               >  1:208586/3‑77 (MQ=255)
                                                                                                                                          tACACAGTGTGTAGATTTCAGTACTGGTGGCTTCCTTTCAATACTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGgcgcg  >  2:254611/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCG  >  NC_002947/4410678‑4410965

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