Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,750,733 C→G 8.1% A31G (GCC→GGC)  PP_1551 → replication protein O

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,750,7330CG 8.1% 148.9 / 2.6 49A31G (GCC→GGC) PP_1551replication protein O
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base G (2/2);  total (24/25)
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.68e-01

AAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTC  >  NC_002947/1750588‑1750879
                                                                                                                                                 |                                                                                                                                                  
aaaaaaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccg                                                                                                                                                 >  6:181114/1‑149 (MQ=255)
    aaaCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTcccggcccgg                                                                                                                                                >  6:294928/1‑146 (MQ=255)
      aCCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGGCCGGGAGAt                                                                                                                                           <  6:145126/149‑1 (MQ=255)
       ccGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGACTGCTGGCTATCGATCTCCCGGCCCGGGAGAt                                                                                                                                           >  3:74503/1‑148 (MQ=255)
         gCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCGACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGaa                                                                                                                                        <  7:139627/149‑1 (MQ=255)
              gCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGAtt                                                                                                                                    <  6:198056/148‑1 (MQ=255)
               cAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATtgt                                                                                                                                  <  1:96240/149‑1 (MQ=255)
                aGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGga                                                                                                                                              >  2:156269/1‑136 (MQ=255)
                aGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGga                                                                                                                                              <  1:156269/136‑1 (MQ=255)
                 gggCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGc                                                                                                                                <  4:36458/149‑1 (MQ=255)
                    cggcTTTCTGTACAGCATTAGAACGGGATAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGt                                                                                                                             <  6:123525/149‑1 (MQ=255)
                        gttCTCTAACGCATTTGAACCGGTTAAAGCAAGGCAAACATCCTCCCACATTACCAGTGCAGGGGGGTAACCCCGGTTGCCAACCAGCGAATTGAAGGCCGGCTGGCTATGGAACTCCCCGGCCGGGAGATGAAGATTGTGCTTTaact                                                                                                                         >  4:130493/2‑146 (MQ=255)
                              tACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGcc                                                                                                                     >  5:84298/1‑147 (MQ=255)
                               aCAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGa                                                                                                                                         <  6:22700/126‑1 (MQ=255)
                               aCAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGa                                                                                                                                         >  5:22700/1‑126 (MQ=255)
                                cAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAg                                                                                                                  <  3:170983/148‑1 (MQ=255)
                                    aTTATAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCAACCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTACCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                              >  5:130982/1‑148 (MQ=255)
                                      tAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                              <  8:197322/146‑1 (MQ=255)
                                      tAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTAACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGcca                                                                                                              >  7:197322/1‑146 (MQ=255)
                                         aaCGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCTTGGccacgtaca                                                                                                        <  5:294928/149‑6 (MQ=255)
                                                                     ccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGGTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCAGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATccc                                                                            <  7:284581/149‑1 (MQ=255)
                                                                     ccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                   <  5:210447/142‑1 (MQ=255)
                                                                     ccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAAc                                                                                   >  6:210447/1‑142 (MQ=255)
                                                                     ccACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATccc                                                                            <  4:230663/149‑1 (MQ=255)
                                                                              aaGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATAGATCTCACGGCCCGGGAGATGAAGATAGTGCTGTAAGTGGCCAACGCCCCCATCAACTACGGTGAGGGCGCTCAACGCATCCCGGCTACCGa                                                                   >  3:5981/1‑149 (MQ=255)
                                                                                gTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGc                                                                                                                <  6:169945/102‑1 (MQ=255)
                                                                                gTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGc                                                                                                                >  5:169945/1‑102 (MQ=255)
                                                                                   cAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcg                                                              <  6:19252/149‑1 (MQ=255)
                                                                                    aGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATcgc                                                             >  8:104951/1‑149 (MQ=255)
                                                                                        ggTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGaaa                                                         <  4:74503/149‑1 (MQ=255)
                                                                                             aCCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCa                                                     >  7:165368/1‑148 (MQ=255)
                                                                                               ccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATcc                                                  >  3:16836/1‑149 (MQ=255)
                                                                                               ccGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATcc                                                  <  6:84298/149‑1 (MQ=255)
                                                                                                     ggACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCt                                             >  7:320517/1‑148 (MQ=255)
                                                                                                     ggACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCt                                             >  6:208116/1‑148 (MQ=255)
                                                                                                       aCAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAc                                          >  4:220299/1‑149 (MQ=255)
                                                                                                        cAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCAcc                                         <  1:172387/149‑1 (MQ=255)
                                                                                                         aaCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGt                                                                                                >  5:312721/1‑93 (MQ=255)
                                                                                                         aaCCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGt                                                                                                <  6:312721/93‑1 (MQ=255)
                                                                                                           ccAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCt                                       >  5:148844/1‑148 (MQ=255)
                                                                                                                  atggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACAcc                                 <  3:130493/147‑1 (MQ=255)
                                                                                                                   tggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGc                                                                                 <  2:147327/98‑1 (MQ=255)
                                                                                                                   tggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGc                                                                                 >  1:147327/1‑98 (MQ=255)
                                                                                                                    ggatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCgtgt                             >  3:264388/1‑149 (MQ=255)
                                                                                                                     gatggCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTc                            <  8:255339/149‑1 (MQ=255)
                                                                                                                                 gCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTc                <  5:208116/149‑1 (MQ=255)
                                                                                                                                 gCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTc                >  1:222807/1‑149 (MQ=255)
                                                                                                                                 aatcttcaTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCAGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTc                >  2:108536/8‑149 (MQ=255)
                                                                                                                                       caTCTCCCGGGCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCa              <  1:108536/144‑1 (MQ=255)
                                                                                                                                        aTCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctg         >  6:126079/1‑149 (MQ=255)
                                                                                                                                         tctcCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgc        >  2:214711/1‑149 (MQ=255)
                                                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAATGCGGTTTCCAGCCtgctgc     <  4:16836/148‑1 (MQ=255)
                                                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCctgctg      <  8:165368/147‑1 (MQ=255)
                                                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCg    >  6:276623/1‑149 (MQ=255)
                                                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCg    >  5:86967/1‑149 (MQ=255)
                                                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCg    <  5:126079/149‑1 (MQ=255)
                                                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCg    >  3:74528/1‑149 (MQ=255)
                                                                                                                                             ccggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCg    >  3:243768/1‑149 (MQ=255)
                                                                                                                                               ggcccggGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTc  <  4:243768/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAAAAAACCGCCTGGCAGGGCGGCTTTCTCTACAGCATTAGAACGGGTTAAAGCATGACAAACATCGTCCCACTTGACAAGTCCAGGGGGTTCACCCGGATGGACAACCAGCTCATGGATGGCCTGCTGGCTATCGATCTCCCGGCCCGGGAGATGAAGATTGTGCTGTACGTGGCCAAGGCCACCATCAACTTCGGTGCGGGCGCTCAACGCATCCCGGCTACCGACATCGCGAAAGCCATCCACGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTC  >  NC_002947/1750588‑1750879

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