Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,693,848 C→G 9.6% D774H (GAC→CAC)  ligD ← DNA ligase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,693,8480CG 9.6% 149.3 / 5.1 50D774H (GAC→CAC) ligDDNA ligase D
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base G (2/3);  total (25/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.31e-01

GCGTACCTTGCCCATGTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGAT  >  NC_002947/3693711‑3693992
                                                                                                                                          |                                                                                                                                                
gCGTACCTTGCCCATGTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGAAAGGTACCGATACTGGCAGCCCGTCACGGGCTCg                                                                                                                                         >  6:214168/1‑147 (MQ=255)
ccGTACCTCGCCCCTGGCCCCCCAAGGGCCGTCGACCTCCCCCCGCTGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGCTTTGCGCCAGTTCCTCACGCCAGACAGGTACCGATCCTGCCAGCCCGTGACGGGCtgcca                                                                                                                                       <  8:255329/148‑5 (MQ=255)
  gTACCTTGCCCATGTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTGACGGGCTGCCAGt                                                                                                                                     <  2:121724/149‑1 (MQ=255)
          cccATGTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGgc                                                                                                                             <  8:165565/149‑1 (MQ=255)
           ccATGTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGgcg                                                                                                                            >  6:66614/1‑149 (MQ=255)
              tGTCTGCCCTAGGGTCGTGGACCTGCGCCAGATGGGCGTGCA‑TGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTGGCAGTGAATAGGCGCAc                                                                                                                         <  2:275520/149‑1 (MQ=255)
               gTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTGACGGGCtgccag                                                                                                                                      <  5:243521/135‑6 (MQ=255)
               gTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTGACGGGCtgccag                                                                                                                                      >  6:243521/1‑130 (MQ=255)
                      cAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCAGGGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTgg                                                                                                                  >  1:48574/1‑149 (MQ=255)
                               tcgACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCtt                                                                                                        <  7:226460/149‑1 (MQ=255)
                               tcgACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCtt                                                                                                        >  4:277964/1‑149 (MQ=255)
                                       cGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGtt                                                                                                  <  4:219283/147‑1 (MQ=255)
                                        gCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTgcg                                                                                               >  5:269564/1‑149 (MQ=255)
                                         ccAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTgcgc                                                                                              <  2:243174/149‑1 (MQ=255)
                                            gATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCTGTTGCGCAgg                                                                                           <  2:170002/149‑1 (MQ=255)
                                               gCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAg                                                                                        >  7:20861/1‑149 (MQ=255)
                                                   gTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGa                                                                                    >  2:255359/1‑149 (MQ=255)
                                                   gTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGa                                                                                    <  8:20861/149‑1 (MQ=255)
                                                   gTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGa                                                                                    >  2:95373/1‑149 (MQ=255)
                                                    tGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCTGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGAt                                                                                   <  5:254009/149‑1 (MQ=255)
                                                             ggTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGt                                                                                          >  6:272168/1‑134 (MQ=255)
                                                             ggTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGt                                                                                          <  5:272168/134‑1 (MQ=255)
                                                                      ccATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCAcgc                                                                  >  6:166568/1‑149 (MQ=255)
                                                                       cATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCAcgc                                                                  >  7:21669/1‑148 (MQ=255)
                                                                            gttggCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGAtttt                                                            >  2:327379/1‑149 (MQ=255)
                                                                             ttggCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTg                                                                                                 <  8:258765/111‑1 (MQ=255)
                                                                             ttggCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTg                                                                                                 >  7:258765/1‑111 (MQ=255)
                                                                                gcggcCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCg                                                         <  5:90754/148‑1 (MQ=255)
                                                                                gcggcCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCg                                                         >  6:90754/1‑148 (MQ=255)
                                                                                 cggcCTTGATTTGCGCCTGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGc                                                       <  6:269564/149‑1 (MQ=255)
                                                                                 cggcCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGc                                                       >  3:144642/1‑149 (MQ=255)
                                                                                  ggcCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGcc                                                      <  5:166568/149‑1 (MQ=255)
                                                                                    ccTTGATTTGCGCCAGTTCCTCACGCCAGAAAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGcccc                                                    <  5:214168/149‑1 (MQ=255)
                                                                                    ccTTGATTTGCGCCAGTTCCGCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGcccc                                                    >  7:178051/1‑149 (MQ=255)
                                                                                        gATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGaca                                                >  4:94748/1‑149 (MQ=255)
                                                                                                  aGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAg                                      <  5:58239/149‑1 (MQ=255)
                                                                                                   gTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAgg                                     >  6:153676/1‑149 (MQ=255)
                                                                                                   gTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAgg                                     <  8:21669/149‑1 (MQ=255)
                                                                                                        tCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTc                                >  6:276147/1‑149 (MQ=255)
                                                                                                           cGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCgg                             <  8:178051/149‑1 (MQ=255)
                                                                                                                  tAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGc                      <  1:327379/149‑1 (MQ=255)
                                                                                                                    ggTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCtt                    <  3:277964/149‑1 (MQ=255)
                                                                                                                     gTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTg                   <  5:276147/149‑1 (MQ=255)
                                                                                                                        ccGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTggt                 <  2:48574/148‑2 (MQ=255)
                                                                                                                            tACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACTGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAgg             >  6:27126/1‑148 (MQ=255)
                                                                                                                               tgcgaGCCCGTGACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAAt         >  3:113154/5‑149 (MQ=255)
                                                                                                                                 gCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAAtt        <  5:66614/148‑1 (MQ=255)
                                                                                                                                   aGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGTGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGAc     >  1:214349/1‑149 (MQ=255)
                                                                                                                                      ccGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGAt  >  6:327232/1‑149 (MQ=255)
                                                                                                                                      ccGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGAt  >  3:201199/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                
GCGTACCTTGCCCATGTCGGCCCAAGGGTCGTCGACCTCCGCCAGATGCTCGTGCA‑GGTTGGTGATGTTCCATTGGTTGGCGGCCTTGATTTGCGCCAGTTCCTCACGCCAGATAGGTACCGATACTGGCAGCCCGTCACGGGCTCGCAGTGAATAGGCGCACGCCGTGGTAGCGCCCTTGCCGTTGCGCAGGTAGTCGATGAAAATGCGTCCCACGCGATTTTTCGGCCCCGACACCGCGCTCAGGCGCTCCGGAAACAGCTTGGCCAGGTAATTGACGAT  >  NC_002947/3693711‑3693992

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