Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,810 T→C 19.0% intergenic (+199/‑648) mngB → / → cydA alpha‑mannosidase/cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,8100TC19.0% 46.3 / 6.0 21intergenic (+199/‑648)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base C (2/2);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGA  >  NC_000913/770689‑770950
                                                                                                                         |                                                                                                                                            
gAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttgaga                                                                                                                                   >  4:1371754/1‑129 (MQ=255)
gAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttgaga                                                                                                                                   <  3:1371754/133‑5 (MQ=255)
    gCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                                                                                                                   <  4:1857539/129‑1 (MQ=255)
    gCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                                                                                                                   >  3:1857539/1‑129 (MQ=255)
                     aaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                                                                                                              <  1:1655190/101‑1 (MQ=255)
                     aaaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                                                                                                              >  2:1655190/1‑101 (MQ=255)
                      aaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAAt                                                                                                   >  4:624468/1‑143 (MQ=255)
                      aaGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAAt                                                                                                   <  3:624468/143‑1 (MQ=255)
                              aaaTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAt                                                                                                                                    <  4:365955/102‑1 (MQ=255)
                              aaaTTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAt                                                                                                                                    >  3:365955/1‑102 (MQ=255)
                                             cAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt                                                                                                                                 <  4:598042/90‑1 (MQ=255)
                                             cAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt                                                                                                                                 >  3:598042/1‑90 (MQ=255)
                                                agaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt                                                                                                                                            <  4:980208/76‑1 (MQ=255)
                                                agaTGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt                                                                                                                                            >  3:980208/1‑76 (MQ=255)
                                                    gatgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                                                                                                              <  3:1821835/70‑1 (MQ=255)
                                                    gatgTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTatat                                                                                                                                              >  4:1821835/1‑70 (MQ=255)
                                                                 gCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                                                                                                                   <  2:1282236/68‑1 (MQ=255)
                                                                 gCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt                                                                                                                                   >  1:1282236/1‑68 (MQ=255)
                                                                                        tttGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGGGTTAATtat                                                                                                <  2:664858/80‑1 (MQ=255)
                                                                                        tttGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGGGTTAATtat                                                                                                >  1:664858/1‑80 (MQ=255)
                                                                                                           taaaGCCACATATACGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtataa                                                                                              <  2:716274/61‑1 (MQ=38)
                                                                                                           taaaGCCACATATACGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATtataa                                                                                              >  1:716274/3‑63 (MQ=37)
                                                                                                                ccACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAg                                    <  2:1425111/116‑1 (MQ=255)
                                                                                                                ccACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAg                                    >  1:1425111/1‑116 (MQ=255)
                                                                                                                 cACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTAtt           <  2:354546/140‑1 (MQ=255)
                                                                                                                 cACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTAtt           >  1:354546/1‑140 (MQ=255)
                                                                                                                 cACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGa  >  2:770982/1‑149 (MQ=255)
                                                                                                                         |                                                                                                                                            
GAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTTAACAAAATAATAACCCTGGGTGAGTTAATTATAATATAATTATAAGTTAACTAAATGTTAATATTGGCGGGGTGGATTTATGCCTTTATTAGTAATCCTGAAACTCTGCGTCGTATTAGCCAGTGA  >  NC_000913/770689‑770950

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