Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,011,800 A→G 11.1% E199G (GAG→GGG)  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,8000AG11.1% 92.0 / 3.3 36E199G (GAG→GGG) repsingle‑stranded DNA dependent ATPase
Reads supporting (aligned to +/- strand):  ref base A (17/15);  new base G (2/2);  total (19/17)
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.39e-01

CCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCGTGGTGGGCGACGACGACCAGTCGAT  >  NC_002947/6011652‑6011944
                                                                                                                                                    |                                                                                                                                                
ccAGAAGCCCTGGTAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAg                                                                                                                                                 >  2:164904/1‑150 (MQ=255)
       cccTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCTCCATCGTCTACTCCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGTACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTgg                                                                                                                                          <  2:96466/150‑1 (MQ=255)
        ccTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCTCGGGGTGCtcct                                                                                                                                         >  2:385038/1‑147 (MQ=255)
        ccTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGa                                                                                                                                         >  1:106454/1‑150 (MQ=255)
         cTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCTCGGGGTGCtcctg                                                                                                                                        <  1:385038/150‑5 (MQ=255)
          tGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAgc                                                                                                                                       >  2:81960/1‑150 (MQ=255)
          tGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAgc                                                                                                                                       >  1:22610/1‑150 (MQ=255)
                                             gccATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGAc                                                                                                    <  1:555284/150‑1 (MQ=255)
                                                    tCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACctg                                                                                                             <  2:430477/134‑1 (MQ=255)
                                                    tCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACctg                                                                                                             >  1:430477/1‑134 (MQ=255)
                                                     cTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCa                                                                                            <  1:528261/150‑1 (MQ=255)
                                                             aCTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGaca                                                                                       >  2:256419/1‑147 (MQ=255)
                                                             aCTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGaca                                                                                       <  1:256419/147‑1 (MQ=255)
                                                                   aGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACcgcg                                                                                                                        >  2:233384/1‑108 (MQ=255)
                                                                   aGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACcgcg                                                                                                                        <  1:233384/108‑1 (MQ=255)
                                                                      gcACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACgccagcc                                                                           >  2:361859/1‑150 (MQ=255)
                                                                            tCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTgg                                                                                                                                  <  1:257130/89‑1 (MQ=255)
                                                                            tCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTgg                                                                                                                                  >  2:257130/1‑89 (MQ=255)
                                                                               aGGCGTTCAACGCGGTGGACTACGACGACCTGACCCAGACGCCGTTCAAGCTGTTCCAGGACAACCCCGAGGTGCTTGGGATGGAGCAGAACCGCGTGAGCTACCAGAGGCGGGGCCAATCCGAGGACACCCACGCCAGCCTGTatcagc                                                                  >  2:39399/1‑145 (MQ=255)
                                                                                ggCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGaa                                                                                                                             <  1:262871/90‑1 (MQ=255)
                                                                                ggCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGaa                                                                                                                             >  2:262871/1‑90 (MQ=255)
                                                                                                      gacgacCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATg                                           >  2:602482/1‑150 (MQ=255)
                                                                                                             tGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGc                                                      >  1:231090/1‑132 (MQ=255)
                                                                                                             tGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGc                                                      <  2:231090/132‑1 (MQ=255)
                                                                                                                  cTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAg                                                                                           <  1:121916/90‑1 (MQ=255)
                                                                                                                  cTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAg                                                                                           >  2:121916/1‑90 (MQ=255)
                                                                                                                         cGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCgtgg                        <  1:316255/150‑1 (MQ=255)
                                                                                                                          ggTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCgtggt                       <  1:464512/150‑1 (MQ=255)
                                                                                                                           gTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCgtggtg                      >  1:531694/1‑150 (MQ=255)
                                                                                                                              aaGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACctgctg                                                                <  2:132580/105‑1 (MQ=255)
                                                                                                                              aaGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACctgctg                                                                >  1:132580/1‑105 (MQ=255)
                                                                                                                                gCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCGTGGTGGGcga                 >  1:372500/1‑150 (MQ=255)
                                                                                                                                        tccaGCACCTCGGGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACctgctg                                                                >  2:154518/4‑95 (MQ=255)
                                                                                                                                        tccaGCACCTCGGGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACctgctg                                                                <  1:154518/92‑1 (MQ=255)
                                                                                                                                              aCCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCGTGGTGGGCGACGACGACCAGTCGa   <  1:602482/150‑1 (MQ=255)
                                                                                                                                               ccccGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCGTGGTGGGCGACGACGACCAGTCGAt  <  1:418785/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCAGAAGCCCTGGAAAAGGCGCGCAACCCGCGCGAGCAGACCGCCGCCATCGTCTACACCCACTATCAGCGCACGCTCAAGGCGTTCAACGCGGTGGACTTCGACGACCTGATCCTGCAGCCGGTCAAGCTGTTCCAGGAGCACCCCGAGGTGCTGGAGCGCTGGCAGAACCGCGTGCGTTACCTGCTGGTGGACGAATACCAGGACACCAACGCCAGCCAGTACCTGCTGGTGAAAATGCTGATCGGCATGCGCAACCAGTTTACCGTGGTGGGCGACGACGACCAGTCGAT  >  NC_002947/6011652‑6011944

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