Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,011,556 C→A 12.2% L118M (CTG→ATG)  rep → single‑stranded DNA dependent ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,011,5560CA12.2% 160.9 / 9.3 57L118M (CTG→ATG) repsingle‑stranded DNA dependent ATPase
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base A (3/4);  total (28/29)
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.24e-01

CACCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATC  >  NC_002947/6011411‑6011702
                                                                                                                                                 |                                                                                                                                                  
cACCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTc                                                                                                                                                >  2:733059/1‑150 (MQ=255)
 aCCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACAGCAACGCCTTGATGTCg                                                                                                                                               <  1:587380/150‑1 (MQ=255)
   cTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACAGCAACGCCTTGATGTCGGa                                                                                                                                             <  2:1448859/150‑1 (MQ=255)
        tgcGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGAcatca                                                                                                                                        >  1:5391/1‑150 (MQ=255)
          cgcCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATg                                                                                                                                      <  2:712431/150‑1 (MQ=255)
           gcCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGc                                                                                                                                     <  2:5391/150‑1 (MQ=255)
             ccGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACAGCAACGCCTTGATGTCGGActc                                                                                                                                          <  1:1132434/143‑3 (MQ=255)
             ccGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACAGCAACGCCTTGATGTCGGActc                                                                                                                                          >  2:1132434/1‑141 (MQ=255)
                          gTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCg                                                                                                                      >  1:99268/1‑150 (MQ=255)
                            ggCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGc                                                                                                                    <  2:585465/150‑1 (MQ=255)
                                  aCCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGTACATCATGCAGAAGGAATACTCCGGCgacgac                                                                                                              >  1:465574/1‑150 (MQ=255)
                                   ccGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGacgacg                                                                                                             <  2:99268/150‑1 (MQ=255)
                                          tACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAgg                                                                                                                               >  1:387618/1‑125 (MQ=255)
                                          tACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAgg                                                                                                                               <  2:387618/125‑1 (MQ=255)
                                           aCGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgac                                                                                                     <  1:733059/150‑1 (MQ=255)
                                            cGTGCCACAACCTGGGCCTGAACAGCATCCGCAAGGAACACGAGCGCCTTGGCTATAAGCCGGGTTTCTCGATCTTCCACGAGTCAGACATCAATGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGGTACGACGGCATcgacg                                                                                                    <  2:810055/150‑1 (MQ=255)
                                             gatCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGc                                                                                                                    <  1:472164/131‑1 (MQ=255)
                                             gatCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGc                                                                                                                    >  2:472164/3‑133 (MQ=255)
                                             gTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgacga                                                                                                   <  2:1445201/150‑1 (MQ=255)
                                                  acaACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgacga                                                                                                   <  2:79447/145‑1 (MQ=255)
                                                  acaACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATcgacga                                                                                                   >  1:79447/1‑145 (MQ=255)
                                                  acaACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCa                                                                                              >  2:288326/1‑150 (MQ=255)
                                                   caACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGAAGGCATCGACGAGATCaa                                                                                             >  1:640809/1‑150 (MQ=255)
                                                      ccTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGAcat                                                                                                                                          <  1:581365/102‑1 (MQ=255)
                                                      ccTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGAcat                                                                                                                                          >  2:581365/1‑102 (MQ=255)
                                                        tGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACa                                                                                        >  2:645954/1‑150 (MQ=255)
                                                             cTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATc                                                                                   <  2:829812/150‑1 (MQ=255)
                                                                aaCATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACCTCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGACCAAGAACATCATCggc                                                                                >  1:701886/1‑149 (MQ=255)
                                                                aaCATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGTACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGt                                                                                <  2:465574/150‑1 (MQ=255)
                                                                aaCATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGt                                                                                <  1:421368/150‑1 (MQ=255)
                                                                  catcatCCGCAAGGAACACGAGCGCCTGGGCTATAAGCGGGGTTTCTCGATCTTCGACGAGTCCGACACCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAAGATGATCGGGGc                                                                              >  1:658423/1‑150 (MQ=255)
                                                                                  cacGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTg                                                              <  1:645954/150‑1 (MQ=255)
                                                                                         gcCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTgc                                                       >  1:317896/1‑150 (MQ=255)
                                                                                          ccTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTgcc                                                      <  2:317896/150‑1 (MQ=255)
                                                                                            tGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCg                                                                                  >  2:1588573/1‑120 (MQ=255)
                                                                                            tGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCg                                                                                  <  1:1588573/120‑1 (MQ=255)
                                                                                              ggCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCa                                                  >  1:1489122/1‑150 (MQ=255)
                                                                                                 tataAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGaa                                               >  2:453265/1‑150 (MQ=255)
                                                                                                      gCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTc                                                                                                                                                <  1:804041/48‑1 (MQ=255)
                                                                                                      gCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTc                                                                                                                                                >  2:804041/1‑48 (MQ=255)
                                                                                                      gCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCt                                          <  1:453265/150‑1 (MQ=255)
                                                                                                       ccGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTg                                         <  2:640809/150‑1 (MQ=255)
                                                                                                                tcGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGgcg                                >  2:1553404/1‑150 (MQ=255)
                                                                                                                     cTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAg                                                                                                                                   <  2:1627115/46‑1 (MQ=255)
                                                                                                                     cTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAg                                                                                                                                   >  1:1627115/1‑46 (MQ=255)
                                                                                                                        cgacgaGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCACTGGAAAAGGCGCGCAAccc                        >  2:914030/1‑150 (MQ=255)
                                                                                                                            gaGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCcgcgc                    <  2:1489122/150‑1 (MQ=255)
                                                                                                                             tgTCCGACAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGgaaga                                                                      <  2:973639/98‑1 (MQ=38)
                                                                                                                             tgTCCGACAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGgaaga                                                                      >  1:973639/2‑99 (MQ=38)
                                                                                                                              ccgatctaTCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACCGTCTCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGa                  <  1:514563/143‑1 (MQ=255)
                                                                                                                                ccGACAGCAACGCCTTGATGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGCACATGATCGGTGCCCGGAAGAACGACCCGCTCCTGCCGCCAGACCACCTGGGAAAACGGCGGAGCCCGCGCGAGc                >  1:2250/1‑150 (MQ=255)
                                                                                                                                 cGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGTATCGACGAGATCACGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCa               >  2:388079/1‑150 (MQ=255)
                                                                                                                                 cGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCa               <  1:914030/150‑1 (MQ=255)
                                                                                                                                 cGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCa               >  2:1402653/1‑150 (MQ=255)
                                                                                                                                   aCATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGCGCAGa             >  2:1235620/1‑150 (MQ=255)
                                                                                                                                      tCAAGGCGTTGCTGTCGAACATCAGGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGGCATGCCGCCATAAGCCCTGGCAAAGGCGCGCAACCCGCGCGAgagatcgg          >  2:514563/1‑143 (MQ=255)
                                                                                                                                              ttGCTGTCGGAAATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATc  <  1:463364/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CACCTTGCTGCGCCCGGGGGAAGGCCGTGGCCTGACCGTGTGTACGTTCCACAACCTGGGCCTGAACATCATCCGCAAGGAACACGAGCGCCTGGGCTATAAGCCGGGTTTCTCGATCTTCGACGAGTCCGACATCAAGGCGTTGCTGTCGGACATCATGCAGAAGGAATACTCCGGCGACGACGGCATCGACGAGATCAAGAACATGATCGGTGCCTGGAAGAACGACCTGGTCCTGCCGCCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATC  >  NC_002947/6011411‑6011702

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: