Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,250,523 A→G 7.3% Y65Y (TAT→TAC)  nth ← DNA glycosylase/apyrimidinic (AP) lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,250,5230AG 7.3% 158.0 / 2.5 56Y65Y (TAT→TAC) nthDNA glycosylase/apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base A (21/29);  new base G (2/2);  total (23/33)
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.88e-01

GCCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGG  >  NC_002947/1250380‑1250665
                                                                                                                                               |                                                                                                                                              
gcCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGAt                                                                                                                                           >  8:99202/1‑149 (MQ=255)
 ccTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATg                                                                                                                                          <  7:99202/149‑1 (MQ=255)
  cTCCAGCGCTTCGCGGGTTTGCGGTACATGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATgg                                                                                                                                         <  6:292217/149‑1 (MQ=255)
   tCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGc                                                                                                                                        >  7:101807/1‑149 (MQ=255)
      aGCGGTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCTGTCGATCAGCAGCCGTCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTc                                                                                                                                     <  7:181111/149‑1 (MQ=255)
      aGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCATCTATGCCCTGGGCGTTGAAGGCCTc                                                                                                                                     <  1:146557/149‑1 (MQ=255)
               cgGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGtt                                                                                                                            >  5:297623/1‑149 (MQ=255)
                       cGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACgg                                                                                                                    <  1:188832/149‑1 (MQ=255)
                        ggTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACggg                                                                                                                   >  2:263431/1‑149 (MQ=255)
                             cTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACa                                                                                                              <  5:47598/149‑1 (MQ=255)
                                gcggTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGc                                                                                                           <  6:131382/146‑1 (MQ=255)
                                    tCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCgg                                                                                                         <  2:167322/147‑1 (MQ=255)
                                         gcgcTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTg                                                                                                   <  2:169857/148‑1 (MQ=255)
                                               gATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCa                                                                                                      >  1:321775/1‑139 (MQ=255)
                                               gATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCa                                                                                                      <  2:321775/139‑1 (MQ=255)
                                                aTCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCttgtt                                                                                           <  6:297623/149‑1 (MQ=255)
                                                aTCAGCAGCCGGCACGCTTATATGACATTCTTGGCCTTTCTGTTGTAAAGCCCGATTGTGTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCgttgtt                                                                                           <  6:15121/149‑1 (MQ=255)
                                                   agcagcCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGAc                                                                                        <  3:111164/149‑1 (MQ=255)
                                                   agcagcCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATAGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGTCCTTGTTGAc                                                                                        <  1:263431/149‑1 (MQ=255)
                                                    gcagcCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACg                                                                                       >  1:73668/1‑149 (MQ=255)
                                                     cagcCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGc                                                                                      >  6:136321/1‑149 (MQ=255)
                                                          ggCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACg                                                                                 >  6:50526/1‑149 (MQ=255)
                                                             aCGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGAc                                                                                  >  8:62704/1‑145 (MQ=255)
                                                             aCGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGAc                                                                                  <  7:62704/145‑1 (MQ=255)
                                                                 ttCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCGTAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTg                                                                           >  3:24752/1‑148 (MQ=255)
                                                                    tATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGAtt                                                                       <  2:73668/149‑1 (MQ=255)
                                                                     aTGACATTATTTGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGAtt                                                                       <  8:259293/148‑1 (MQ=255)
                                                                        aCATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAAc                                                                                                               >  6:86338/1‑105 (MQ=255)
                                                                        aCATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAAc                                                                                                               <  5:86338/105‑1 (MQ=255)
                                                                          aTTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGcgcg                                                                  >  3:226675/1‑148 (MQ=255)
                                                                               tGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCg                                                                              <  2:299854/131‑1 (MQ=255)
                                                                               tGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCg                                                                              >  1:299854/1‑131 (MQ=255)
                                                                                 gCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAgtgt                                                          <  5:50526/149‑1 (MQ=255)
                                                                                    ttGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCAc                                                       >  8:175842/1‑149 (MQ=255)
                                                                                        tgttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGTTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGc                                                    >  3:203154/1‑148 (MQ=255)
                                                                                        tgttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCTACGTGGAACAGGCGGGCAGTGGCCTTGTTGAATCCGCCGTAGGTGGATTGCGCGCACAGTGTCACGGc                                                    >  3:60955/1‑148 (MQ=255)
                                                                                        tgttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAgt                                                            <  5:19797/140‑1 (MQ=255)
                                                                                        tgttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAgt                                                            >  6:19797/1‑140 (MQ=255)
                                                                                        tgttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGc                                                    >  1:171132/1‑148 (MQ=255)
                                                                                          ttgtAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCGTAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGAc                                                 <  4:24752/149‑1 (MQ=255)
                                                                                             tAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTAGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACcag                                              >  7:128537/1‑149 (MQ=255)
                                                                                               aaGCCGGATTGTCTTGAGGTATTCGCTGTGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACcagca                                            <  4:60955/149‑1 (MQ=255)
                                                                                               aaGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCgg                                                                                                                                   <  5:29385/62‑1 (MQ=255)
                                                                                               aaGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCgg                                                                                                                                   >  6:29385/1‑62 (MQ=255)
                                                                                                 gCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGTTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGt                                          <  4:203154/149‑1 (MQ=255)
                                                                                                 gCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGt                                          <  2:171132/149‑1 (MQ=255)
                                                                                                          tCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCTGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACcagca                                            >  4:36771/1‑138 (MQ=255)
                                                                                                          tCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCTGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACcagca                                            <  3:36771/138‑1 (MQ=255)
                                                                                                              gATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACgggg                                <  1:15831/146‑1 (MQ=255)
                                                                                                                  tATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGAGAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTa                         <  4:226675/149‑1 (MQ=255)
                                                                                                                        cTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAACCGGGGTCGTGTAGGGCAg                   >  4:184867/1‑149 (MQ=255)
                                                                                                                          gAGGCCATCTATGCACTGGGCGTTGAAGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAttt                 >  7:273117/1‑149 (MQ=255)
                                                                                                                           aGGCCTTCAACGCCCAGTGCTTAGATGGCGTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGtt                 <  8:128537/148‑1 (MQ=255)
                                                                                                                              ccTGCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGt             <  6:324175/149‑1 (MQ=255)
                                                                                                                                    aagCCCAGGGCCAAGATGGCCGCGGGGGTGTTGGCGACGCTGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTgg       <  8:236328/147‑1 (MQ=255)
                                                                                                                                         ccgggCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTGAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCgg  <  3:184867/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GCCTCCAGCGCTTCGCGGGTTTGCGGTACCTGGCTGTCGTGGCGCTCGATCAGCAGCCGGCACGCTTCTATGACATTCTTGGCCTTGCTGTTGTAAAGCCCGATTGTCTTGATGTATTCGCTGAGGCCTTCAACGCCCAGGGCATAGATGGCCTCGGGGGTGTTGGCGACGGGGAACAGGCGGGCAGTGGCCTTGTTGACGCCGACGTCGGTGGATTGCGCGGACAGTGTCACGGCGACCAGCAGTTCAAACGGGGTCGTGTAGGCCAGTTCGGTCTTTGGGTCGG  >  NC_002947/1250380‑1250665

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