Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,623,965 G→T 8.8% Q171H (CAG→CAT)  PP_2295 → antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,623,9650GT 8.8% 140.9 / 2.3 45Q171H (CAG→CAT) PP_2295antirestriction protein
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base T (2/2);  total (24/21)
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.97e-01

CCTATCATGGCGGTTCTGTTACAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGAACTACCGGGGGCTTTTTCTATCGCCGA  >  NC_002947/2623817‑2624107
                                                                                                                                                    |                                                                                                                                              
ccTATCATGGCGGTTCTGTTACAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGgcagc                                                                                                                                               >  1:1017801/1‑150 (MQ=255)
      aTGGCGGTTCTGTTACAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCtt                                                                                                                                         >  2:311962/1‑150 (MQ=255)
        ggCGGTTCTGTTACAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTgg                                                                                                                                       <  1:311962/150‑1 (MQ=255)
                 gTTACAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGAt                                                                                                                              >  2:384827/1‑150 (MQ=255)
                 gTTACAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGAt                                                                                                                              <  2:1180973/150‑1 (MQ=255)
                     cAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGt                                                                                                                          <  2:1017801/150‑1 (MQ=255)
                        aagaTGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTccc                                                                                                                       >  2:777939/1‑150 (MQ=255)
                           aTGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCcaca                                                                                                                    >  2:962053/1‑150 (MQ=255)
                           aTGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCcaca                                                                                                                    >  2:272455/1‑150 (MQ=255)
                             gACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGAAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCcacaca                                                                                                                  >  1:633528/1‑150 (MQ=255)
                                ttCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGc                                                                                                                                 <  2:988313/132‑1 (MQ=255)
                                ttCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGc                                                                                                                                 >  1:988313/1‑132 (MQ=255)
                                       ccgccTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATc                                                                                                        <  1:384827/150‑1 (MQ=255)
                                                 gggCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTa                                                                                              >  2:13199/1‑150 (MQ=255)
                                                    ccATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACgt                                                                                           <  1:962053/150‑1 (MQ=255)
                                                    ccATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACgt                                                                                           <  1:272455/150‑1 (MQ=255)
                                                       ttACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGa                                                                                                    >  2:1160681/1‑138 (MQ=255)
                                                       ttACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGa                                                                                                    <  1:1160681/138‑1 (MQ=255)
                                                            aCAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTa                                                                                   <  1:13199/150‑1 (MQ=255)
                                                             cAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAg                                                                                  >  2:601615/1‑150 (MQ=255)
                                                                          ttCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTg                                                                            >  1:954920/1‑143 (MQ=255)
                                                                          ttCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTg                                                                            <  2:954920/143‑1 (MQ=255)
                                                                             gCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGACAGGAGATGAATCTTTCTTTAACTGGATCGAGTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGa                                                                  >  2:134355/1‑150 (MQ=255)
                                                                                gaggagTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGa                                                                                                                               >  1:223165/1‑86 (MQ=255)
                                                                                gaggagTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGa                                                                                                                               <  2:223165/86‑1 (MQ=255)
                                                                                        gtgCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTg                                                                            >  1:924860/1‑129 (MQ=255)
                                                                                        gtgCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTg                                                                            <  2:924860/129‑1 (MQ=255)
                                                                                        gtgCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATtgt                                                                    >  2:48544/1‑137 (MQ=255)
                                                                                        gtgCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATtgt                                                                    <  1:48544/137‑1 (MQ=255)
                                                                                            tgcGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAAt                                                   >  2:87833/1‑150 (MQ=255)
                                                                                                   ttCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCATCAGGCTTg                                                                                                                                        >  2:224550/1‑58 (MQ=255)
                                                                                                   ttCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCATCAGGCTTg                                                                                                                                        <  1:224550/58‑1 (MQ=255)
                                                                                                   ttCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCATCAGGCTTg                                                                                                                                        <  1:825741/58‑1 (MQ=255)
                                                                                                   ttCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCATCAGGCTTg                                                                                                                                        >  2:825741/1‑58 (MQ=255)
                                                                                                      aCGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAgcg                                         >  1:1195208/1‑150 (MQ=255)
                                                                                                           cTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCt                                                                                                           >  2:329552/1‑79 (MQ=255)
                                                                                                           cTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCt                                                                                                           <  1:329552/79‑1 (MQ=255)
                                                                                                                     aTGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGAACt                          >  2:472891/1‑150 (MQ=255)
                                                                                                                          tctttcttTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGAACTACCgg                     <  1:472891/150‑1 (MQ=255)
                                                                                                                                   aaCTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGAACTACCGGGGGCttttt            <  1:601615/150‑1 (MQ=255)
                                                                                                                                      tGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAActcc                                                          >  1:227629/1‑101 (MQ=255)
                                                                                                                                      tGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAActcc                                                          <  2:227629/101‑1 (MQ=255)
                                                                                                                                         aTCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGa                             >  2:223860/1‑127 (MQ=255)
                                                                                                                                         aTCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGa                             <  1:223860/127‑1 (MQ=255)
                                                                                                                                             aCTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGAACTACCGGGGGCTTTTTCTATCGCcga  >  2:697580/1‑150 (MQ=255)
                                                                                                                                             aCTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGAACTACCGGGGGCTTTTTCTATCGCcga  <  1:87833/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCTATCATGGCGGTTCTGTTACAGAAGATGACTTCCAGGCCGCCTACCGGGGCCATTACGACAGCGAAAAAGCCTTCGCAGAGGAGTGGTGCTGCGAGCTTCACGGCCTGAAAGGAGATGAATCTTTCTTTAACTGGATCGACTGGCAGCAGGCTTGGGACTGCGATCTGTCCCACACATTCGTCTATCAGGACGGCTACGTGTTCCTTAGCTCCTGGTAATTGTGACGAACTCCTCAAAATCGTTTCAGCGCTGCCCTCGGGAACTACCGGGGGCTTTTTCTATCGCCGA  >  NC_002947/2623817‑2624107

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