Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,751,022 C→G 8.9% D113E (GAC→GAG)  PP_1551 → replication protein O

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,751,0220CG 8.9% 163.0 / 5.2 55D113E (GAC→GAG) PP_1551replication protein O
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base G (3/2);  total (28/27)
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.86e-01

CGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAA  >  NC_002947/1750875‑1751171
                                                                                                                                                   |                                                                                                                                                      
cGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAATGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTc                                                                                                                                                      >  1:1362484/1‑150 (MQ=255)
 gCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCg                                                                                                                                                     >  2:771477/1‑150 (MQ=255)
       ccTGACACCGTGTCCAAGGCGGTTTCCGGCGTGCTGGGTCGCCGTGTCCTGTTCCGCGACAGTTTTAGTTCGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                               <  1:386209/150‑1 (MQ=255)
       ccTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                               <  2:1328637/150‑1 (MQ=255)
       ccTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGaa                                                                                                                                               <  2:1271586/150‑1 (MQ=255)
         tGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCGGCTgtttt                                                                                                                                             <  1:771477/150‑5 (MQ=255)
            cacCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGt                                                                                                                                          <  2:1168746/150‑1 (MQ=255)
             acCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGGGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTc                                                                                                                                         >  1:223202/1‑150 (MQ=255)
                     cAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCgg                                                                                                                                 >  2:70881/1‑150 (MQ=255)
                      aaGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGc                                                                                                                                <  1:377394/150‑1 (MQ=255)
                       aGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCGGCTgttttg                                                                                                                                            <  2:375603/137‑6 (MQ=255)
                       aGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGAGTCGGCTgttttg                                                                                                                                            >  1:375603/1‑132 (MQ=255)
                         gCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCa                                                                                                                             <  2:313872/150‑1 (MQ=255)
                                 cAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATGGGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaa                                                                                                                     >  1:728669/1‑150 (MQ=255)
                                 cAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGaa                                                                                                                     <  1:13111/150‑1 (MQ=255)
                                      tgctgcGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGa                                                                                                                <  2:334162/150‑1 (MQ=255)
                                               gCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAAccg                                                                                                       >  1:1184451/1‑150 (MQ=255)
                                                  gtgtCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAg                                                                                                                 <  2:1242897/137‑1 (MQ=255)
                                                  gtgtCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAg                                                                                                                 >  1:1242897/1‑137 (MQ=255)
                                                     tCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTc                                                                                                                                                      <  1:805811/97‑1 (MQ=255)
                                                     tCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTc                                                                                                                                                      >  2:805811/1‑97 (MQ=255)
                                                     tCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACt                                                                                                 >  2:1101508/1‑150 (MQ=255)
                                                          ttCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTg                                                                                                                                                 <  2:1387397/97‑1 (MQ=255)
                                                          ttCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTg                                                                                                                                                 >  1:1387397/1‑97 (MQ=255)
                                                                     tggtgCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaa                                                                                 >  2:688868/1‑150 (MQ=255)
                                                                         gCCCGGGGTGCCACCGGCGTGAACGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAATCCGCCGACTCCCTTCTTTATTCTaagaaa                                                                             <  2:678081/150‑1 (MQ=255)
                                                                            cGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAagaaagaa                                                                          >  1:893408/1‑150 (MQ=255)
                                                                             ggggTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAgaaagaaa                                                                         >  1:754811/1‑150 (MQ=255)
                                                                              gggTGACATCGGCGTGAATGACCCGAATGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAAc                                                                        <  2:1362484/150‑1 (MQ=255)
                                                                                     aTCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTAt                                                                 <  1:70881/150‑1 (MQ=255)
                                                                                      tCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTAtt                                                                                     >  1:1050769/1‑129 (MQ=255)
                                                                                      tCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTAtt                                                                                     <  2:1050769/129‑1 (MQ=255)
                                                                                      tCGGCGTGAATCACCCGAAGGACTGGCTTTACGTTGATCCGAAACAGACCAAAACAACCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATg                                                                >  2:1093746/1‑150 (MQ=255)
                                                                                              aaagACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATctt                                                       >  1:1425931/4‑150 (MQ=255)
                                                                                              aaTGACCCGAAGGACTGGATTTACGTTGATCCTAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATctt                                                       >  2:218358/1‑150 (MQ=255)
                                                                                                 gACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGa                                                                                                   <  2:1093733/104‑1 (MQ=255)
                                                                                                 gACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGa                                                                                                   >  1:1093733/1‑104 (MQ=255)
                                                                                                   cccGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATcttccttc                                                  <  1:688868/150‑1 (MQ=255)
                                                                                                   cccGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAAAAGACCAAAACCGCCGATTCCCTTCTTTATTCTAAGTAAGAAACCCCCTATGTAA‑ATcttccttc                                                  <  2:163013/150‑1 (MQ=255)
                                                                                                    ccGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCg                                                 >  2:1331862/1‑150 (MQ=255)
                                                                                                          gACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGCAGCCCCCTATGTAA‑ATCTTCCTTCGGAAGAt                                           <  1:37569/150‑1 (MQ=255)
                                                                                                            cTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGt                                         <  2:1184451/150‑1 (MQ=255)
                                                                                                            cTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGt                                         >  1:812993/1‑150 (MQ=255)
                                                                                                               gATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCTAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTAc                                      <  2:223202/150‑1 (MQ=255)
                                                                                                                  ttACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATg                                   <  2:893408/150‑1 (MQ=255)
                                                                                                                          aTCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAACATCTTCCGAAGGAAGAT‑TTACAtagg                                 >  2:66768/1‑141 (MQ=255)
                                                                                                                                  cAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccg                   <  1:1331862/150‑1 (MQ=255)
                                                                                                                                  cAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCcagccg                   >  2:805693/1‑150 (MQ=255)
                                                                                                                                    gACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCa                           <  2:659038/140‑1 (MQ=255)
                                                                                                                                    gACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCa                           >  1:659038/1‑140 (MQ=255)
                                                                                                                                     aCCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccggag                >  2:593255/1‑150 (MQ=255)
                                                                                                                                       cAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCgagccggagcc              >  2:185762/1‑150 (MQ=255)
                                                                                                                                          ttcAGCCGAGTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGGAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATACCCCCCCAGCGAGCCGGAGCCGAc           >  1:1078114/3‑150 (MQ=255)
                                                                                                                                              gCCGACTCGGCTGAAGTGTTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCg       <  2:812993/150‑1 (MQ=255)
                                                                                                                                                 gACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCCTCGGAAAATGTTACATGCCCCCCCAGCCAGCCGGAGCCGACCCCGGcc    >  2:1324681/1‑150 (MQ=255)
                                                                                                                                                   cTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGCTGTTACATGCCCCCCCGGCGAGCCGGAGCCGACTCCGGCCaa  <  1:1101508/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGCTCACCCTGACACCGTGTCCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAA‑ATCTTCCTTCGGAAGATGTTACATGCCCCCCCAGCGAGCCGGAGCCGACTCCGGCCAA  >  NC_002947/1750875‑1751171

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