Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,091,006 +C coding (660/708 nt) endA → DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,091,0051.C92.9% 32.6 / ‑3.4 14coding (659/708 nt)endADNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (9/4);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATG  >  NC_000913/3090777‑3091203
                                                                                                                                                                                                                                      |                                                                                                                                                                                                      
cccgcccgtcgtacgcatcatttcgaactgcaatccgtgcgtcgcaacaccttcccctaccgctgcaagtgccaggagcatcagcttaccgtacgccgccataatcgcgtagttcgtggcgatgccgtctatcgctgtgttcactgcggtgaacagctggttgcgaaagaaccatctgaactagcattaactgtcctgatctggctgatttcataccaaaGGTGCAGGGCCAATCATAACCCGTatgtgn                                                                                                                                                                                    <  1:86846‑M2/30‑2 (MQ=255)
                                                             gctgcaagtgccaggagcatcagcttaccgtacgccgccataatcgcgtagttcgtggcgaggccgtctatcgctgtgttcactgcggtgaacagctggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCt                                                                                                                      >  1:337126‑M2/160‑251 (MQ=255)
                                                                       ccaggagcatcagcttaccgtacgccgccataatcgcgtagttcgtggcgaggccgtctatcgctgtgttcactgcggtgaacagctggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTCACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAcgcg                                                                                                            >  1:159705‑M2/150‑251 (MQ=255)
                                                                                tcagcttaccgtacgccgccataatcgcgtagttcgtggcgaggccgtctatcgatgtgttcactgcggtgaacagctggttgcgacataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACcgc                                                                                                   >  1:211959‑M2/141‑251 (MQ=255)
                                                                                                 gccataatcgcgtagttcgtggcgaggccgtctatcgctgtgttcactgcggtgaacagctggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCATGGCCAATCATAACCCTTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCttttt                                                                                                                                <  1:143904‑M2/82‑1 (MQ=255)
                                                                                                                   gtggcgaggccgtctatcgctgtgttcactgcggtgaacagctggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGc                                                                                                                                                                         <  1:20625‑M2/41‑1 (MQ=255)
                                                                                                                                    nnctgtgtnnactgcggtgaacagctggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGc                                                                                                                                                                                                       <  2:14520‑M2/10‑1 (MQ=255)
                                                                                                                                    cgctgtgttcactgcggtgaacagctggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGc                                                                                                                                                                                                       >  1:14520‑M2/89‑98 (MQ=255)
                                                                                                                                                      tnngagctgnntgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTAtgtg                                                                                                                                                                                     <  2:30035‑M2/29‑1 (MQ=255)
                                                                                                                                                      gaacagctggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcagaccaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACAGTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATTCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGcc                             <  1:359132‑M2/181‑1 (MQ=255)
                                                                                                                                                              ggttgcgaaataaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTAtgtg                                                                                                                                                                                     >  1:30035‑M2/63‑91 (MQ=255)
                                                                                                                                                                   cgaaatcaccatctgaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCAATCATCACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAAc                                                                                                                          >  1:96928‑M2/58‑145 (MQ=255)
                                                                                                                                                                    gaaataaccatctgaactatcaggaactttcctgctctggctgattgcatacaaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGCCCGGCGCGaa                                                                                                                                                               >  1:391874‑M2/57‑107 (MQ=255)
                                                                                                                                                                                 gaactatcaggaactttcctgatctggctgattgcataccaaaGGTGCAGGGCCCATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGACCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGCCGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGCACCACTGACCAGCCATTCTCACATCGCGCTTTGCGACGATg  >  1:377487‑M2/44‑251 (MQ=255)
                                                                                                                                                                                                          ggctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTCACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTCCCACCGCGTgggg                                                                                             >  1:432256‑M2/19‑135 (MQ=255)
                                                                                                                                                                                                            ctgattgcataccaaaGGTGCAGGGCCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGACAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAc                                                                                                               >  1:202982‑M2/17‑115 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                      
GTGAAGGCCAGTACGGTCAATGCGCCATGAAGGTCGATTTCAAAGAAAAAGCTGCCGAACCACCAGCGCGTGCACGCGGTGCCATTGCGCGCACCTACTTCTATATGCGCGACCAATACAACCTGACACTCTCTCGCCAGCAAACGCAGCTGTTCAACGCATGGAACAAGATGTATCCGGTTACCGACTGGGAGTGCGAGCGCGATGAACGCATCGCGAAGGTGCAGGGCAATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACCGCGTGGGGAGACGACGCGGATTTTTAACTATGCGTATCCCCCGCATTTATCATCCTGAACCACTGACCAGCCATTCTCACATCGCGCTTTGCGAAGATG  >  NC_000913/3090777‑3091203

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: