Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,410,778 G→C 10.3% 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,7780GC10.3% 102.8 / 4.0 39H121Q (CAC→CAG) PP_3894replication protein O
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base C (2/2);  total (19/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 = 1.21e-01

CCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCG  >  NC_002947/4410639‑4410919
                                                                                                                                           |                                                                                                                                             
ccTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTCTCAACACAGTGGCAGGGGTGTTTTCTACACACTGTGTAGaa                                                                                                                                      <  1:305052/149‑2 (MQ=255)
ccTTGGCGTTGCCGTGACCGGCTTCCTCGGCGCCGGTCGATGCGGATCCGTCGACAACGAATGCCTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACACTGTGTAGaa                                                                                                                                      <  6:257999/149‑2 (MQ=255)
                  cGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTg                                                                                                                    >  7:230714/1‑149 (MQ=255)
                         cTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGa                                                                                                             >  6:325157/1‑149 (MQ=255)
                                ccGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCa                                                                                                      >  6:293898/1‑149 (MQ=255)
                                           ggatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCg                                                                                           <  8:230714/149‑1 (MQ=255)
                                            gatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTg                                                                                                                    >  8:21942/1‑123 (MQ=255)
                                            gatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTg                                                                                                                    <  7:21942/123‑1 (MQ=255)
                                                 gTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGa                                                                                     <  4:130169/149‑1 (MQ=255)
                                                     aCAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTCTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGt                                                                                 >  6:14231/1‑149 (MQ=255)
                                                         ccAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTTATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTTTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCCCTCAGGTCCCGGTAGGAAAGTGAtt                                                                             >  8:300536/1‑149 (MQ=255)
                                                         ccAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGAtt                                                                             <  5:325157/149‑1 (MQ=255)
                                                          cAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTg                                                                            >  6:215385/1‑149 (MQ=255)
                                                           aaTGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTgg                                                                           <  2:71810/149‑1 (MQ=255)
                                                                             tttgtatttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATc                                                         >  8:245649/1‑149 (MQ=255)
                                                                                       tGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAg                                                >  7:87547/1‑148 (MQ=255)
                                                                                        gTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCa                                              >  1:28551/1‑149 (MQ=255)
                                                                                                    tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACgg                                   <  8:87547/148‑1 (MQ=255)
                                                                                                    tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGa                                  >  4:221842/1‑149 (MQ=255)
                                                                                                      tAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGAt                                 <  3:221842/148‑1 (MQ=255)
                                                                                                       aaTTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATc                                <  3:181354/148‑1 (MQ=255)
                                                                                                              aaCACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACgg                                   <  5:336431/138‑1 (MQ=255)
                                                                                                              aaCACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACgg                                   >  6:336431/1‑138 (MQ=255)
                                                                                                              aaCACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAAtt                        >  8:233081/1‑149 (MQ=255)
                                                                                                               acacAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTg                       <  5:293898/149‑1 (MQ=255)
                                                                                                                  caGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATc                    <  7:233081/149‑1 (MQ=255)
                                                                                                                   aGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCg                   <  5:14231/149‑1 (MQ=255)
                                                                                                                       gCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCg               <  5:215385/149‑1 (MQ=255)
                                                                                                                           gggTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGaa                          >  3:91489/1‑134 (MQ=255)
                                                                                                                           gggTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGaa                          <  4:91489/134‑1 (MQ=255)
                                                                                                                             ttgaaatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGacac         >  1:7642/7‑149 (MQ=255)
                                                                                                                             ttgaaatCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGacac         >  1:2347/7‑149 (MQ=255)
                                                                                                                               gtTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCCGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGt       <  7:300536/149‑1 (MQ=255)
                                                                                                                                  ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCg                   <  4:192850/134‑1 (MQ=255)
                                                                                                                                  ttCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCg                   >  3:192850/1‑134 (MQ=255)
                                                                                                                                    cTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGgcg  >  7:277542/1‑149 (MQ=255)
                                                                                                                                    cTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGgc   <  7:245649/148‑1 (MQ=255)
                                                                                                                                     tACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTg                                                                            <  4:107238/74‑1 (MQ=255)
                                                                                                                                     tACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTg                                                                            >  3:107238/1‑74 (MQ=255)
                                                                                                                                           |                                                                                                                                             
CCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCG  >  NC_002947/4410639‑4410919

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