Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,751,080 T→A 14.6% S147T (TCG→ACG) ‡ PP_1551 → replication protein O

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,751,0800TA14.6% 75.6 / 9.4 42S147T (TCG→ACG) ‡PP_1551replication protein O
Reads supporting (aligned to +/- strand):  ref base T (10/23);  new base A (4/2);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

GACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGC  >  NC_002947/1750939‑1751226
                                                                                                                                              |                                                                                                                                                   
gACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGa                                                                                                                                             <  1:5495/149‑1 (MQ=255)
 aCTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTACATCTTCCGAAGGAAGAt                                                                                                                                            <  1:281977/149‑1 (MQ=255)
   tGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTACGTCGGAAGATGt                                                                                                                                          >  1:202912/1‑149 (MQ=255)
     gATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATCCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGtt                                                                                                                                         <  1:267668/148‑1 (MQ=255)
       ggtccGTTGATCCGACACAGACCAAAACAGCCGACTCGGCTGAAGTTGTCCGAATCGTCTCAGAGTCGAACCAGACCAAAACCGCCGACTCCCTTCTTTATTCTCAGAAAGACACCCCCTATGTAAATCTTCCTTCGGAAGATTTACAt                                                                                                                                     <  5:294009/145‑1 (MQ=255)
                 tCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAGCCCCCCAATGTAAATCTTGCGTGGGGAGATGGTACATGCCCCCCCa                                                                                                                            >  8:348167/1‑149 (MQ=255)
                  ccGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCGCCTCTGGAAATCTGGCTGGGGAAGATGGTTCATGAGACCCCAg                                                                                                                           >  6:134177/1‑149 (MQ=255)
                    gAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCg                                                                                                                         <  5:321918/149‑1 (MQ=255)
                          gACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCCGGAAGATGTTACATGCCCCCCCAGCgagccg                                                                                                                    <  1:45561/149‑1 (MQ=255)
                             cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGACACGGCCTAGGTTAATCTTCCTGCGTAAGATGTTACATGGCCCGCCAGCgcgccggcgc                                                                                                                >  8:158715/1‑149 (MQ=255)
                             cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCGATGTAAATCTTACTGCGGAAGATGTTACATTCACACCCAGCgagccggcgc                                                                                                                >  2:139830/1‑149 (MQ=255)
                             aaaaaCCGCAGAGTCGGCTTAAGTGTTCAGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAATAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccggagc                                                                                                                <  8:302742/148‑1 (MQ=255)
                                       gACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATATTGCTTCGGAAGATGGTAGATGACCACCCACCGCGCCCGCGCCGACGCCCGc                                                                                                      >  2:224652/1‑149 (MQ=255)
                                            ggCTGAAGTGGTCCTAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAgg                                                                                                 <  3:216052/149‑1 (MQ=255)
                                            ggCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAg                                                                                                  <  5:263798/148‑1 (MQ=255)
                                              cTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAg                                                                                                  <  1:62498/146‑1 (MQ=255)
                                              cTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGAt                                                                                               <  3:312705/149‑1 (MQ=255)
                                                 atttgttCCGAATCGGATCAGAGTCTAAACATACAAAAACTGCAGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCACCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGa                                                                                             <  8:337732/142‑1 (MQ=255)
                                                 aaGTGGTCCGAATCGGCTATGAGTAGAAACAGTCCAAAAACGCCGACTCCCTTCTTTATTCTAATAAAGAAACCCCCAATGTAAATATTCCTTCGGAATATGTTACATGCCCACCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGa                                                                                             <  6:35453/148‑1 (MQ=255)
                                                                    agagTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCAAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCgg                                                                           <  7:348167/147‑1 (MQ=255)
                                                                    agagTCGAAACACACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCgg                                                                           <  4:262994/147‑1 (MQ=255)
                                                                      agTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCaaa                                                                       <  5:296310/149‑1 (MQ=255)
                                                                       gTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCaaa                                                                       <  5:134177/148‑1 (MQ=255)
                                                                             aCAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCa                                                                <  6:187210/149‑1 (MQ=255)
                                                                                  ccAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAGCGCCCGATGTAAATCTTGGTTGGGGAGATGTTACATGAGGCCCAAGCGAGGCGGGGCCGCGACCGGCGAAGGAGGAGCGCAAGGCGACCAACGGGAAAgccgcacagcgg                                                           >  4:177597/1‑142 (MQ=255)
                                                                                   cAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAGGCGCCCGATGTAAATCTTACTTCGGCAGATGTTAGATGCACCCCCAGCGACCCGGAGCCGACGCCGGCCAAGCAGGCGCGCAAGGCGCACTGCGCCCAAGCCGCCAAGCTgg                                                          >  4:203708/1‑149 (MQ=255)
                                                                                            cgACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATc                                                 <  7:158715/149‑1 (MQ=255)
                                                                                                  ccTTCTTTATTCTAAGAAAGAGGCGCCCTATGTAAATCTTCATTCGGAAGAGGTTACATGACCCCCCGGCGACCCGCCGCCGACACCGGCCCAGGAGGAGCGGCAGGCGCCCATCGGCAAAGCCGCCAGGCTGGCCCGCAAGCCGCCCg                                           >  2:124402/1‑149 (MQ=255)
                                                                                                  ccTTCTTTATTCTAAGAAAGAACCCCCCTATGTAAATCTTGCTTAGGAAGATGTTACAAGCCAACACAGCGAGCCGGAGCCGACTACGGCCAAGGGAGACCGAAAGGCACCCTACGGCAAACACCCCATGCAGGCCGACAAACCGCAcc                                           >  6:381288/1‑148 (MQ=255)
                                                                                                  ccTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGAt                                                                                                                                            <  4:362277/52‑1 (MQ=255)
                                                                                                  ccTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGAt                                                                                                                                            >  3:362277/1‑52 (MQ=255)
                                                                                                              tAAGAAAGAAGCCCCCTATGTAAATCTTCCGTGGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGGAAGGCGCCCTTCGGCGAAGCCGCCACGCAGGCCGACCATCCGCACGGGCGCGaggcat                               >  8:285810/1‑144 (MQ=255)
                                                                                                                        aCCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGc                     <  3:177597/149‑1 (MQ=255)
                                                                                                                            ccTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGAt                 >  5:211033/1‑149 (MQ=255)
                                                                                                                               aTGTAACATCTTCCGAAGGAAGATTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTAc              >  5:124671/1‑149 (MQ=255)
                                                                                                                                tgtaacaTCTTCCGAAGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTg                   >  8:286275/7‑144 (MQ=255)
                                                                                                                                tgtaacaTCTTCCGAAGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTg                   <  7:286275/138‑1 (MQ=255)
                                                                                                                                 gtaacaTCTTCCGAAGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACt             >  5:195193/6‑149 (MQ=255)
                                                                                                                                    aaTCTTCGTTCGGTAGATGTTACTTGGCCCCCCAGGGAGCGGGAGGCGAGTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGGAAAGCGGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTgg          <  5:307229/148‑1 (MQ=255)
                                                                                                                                    aaTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTgg          <  2:202912/148‑1 (MQ=255)
                                                                                                                                     caTCTTCCGAAGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGc         >  8:351939/2‑149 (MQ=255)
                                                                                                                                      aTCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCt        <  6:211033/149‑1 (MQ=255)
                                                                                                                                           ccttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGCCTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCGGATCGCTGATTACTTGGCTGTCCg   >  3:42284/1‑149 (MQ=255)
                                                                                                                                           ccttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGa                  >  2:162726/1‑134 (MQ=255)
                                                                                                                                           ccttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGa                  <  1:162726/134‑1 (MQ=255)
                                                                                                                                           ccttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCg   >  3:46499/1‑149 (MQ=255)
                                                                                                                                             ttcGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGc  >  8:208918/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAAGGATGAGCGCAAGGCTCCCTTCGGCAAAGCCGCCATGCTGGCCGACAATCCGCACGGGCTCGATGAATCCCTGATCGCTGATTACTTGGCTGTCCGC  >  NC_002947/1750939‑1751226

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