Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,682 A→C 17.5% intergenic (+71/‑776) 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,6820AC17.5% 95.5 / 19.5 40intergenic (+71/‑776)mngB/cydAalpha‑mannosidase/cytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (15/18);  new base C (3/4);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

AATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTT  >  NC_000913/770540‑770824
                                                                                                                                              |                                                                                                                                              
aaTCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCAt                                                                                                                                          <  4:172772/149‑1 (MQ=255)
aaTCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCAt                                                                                                                                          <  2:305886/149‑1 (MQ=255)
  ctaaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGa                                                                                                                                        <  1:502678/147‑1 (MQ=255)
   cAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGaa                                                                                                                                       <  3:365264/149‑1 (MQ=255)
    aaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAAc                                                                                                                                      <  3:104396/149‑1 (MQ=255)
      ggTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCTACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                                                                    <  3:82952/149‑1 (MQ=255)
        ttCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCa                                                                                                                                           >  1:269823/1‑140 (MQ=255)
        ttCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCa                                                                                                                                           <  2:269823/140‑1 (MQ=255)
        ttCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                                                                                  >  2:127730/1‑149 (MQ=255)
                 aTCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                                                                         <  4:454450/149‑1 (MQ=255)
                      ggCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                                                                    <  1:547132/133‑1 (MQ=255)
                      ggCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGc                                                                                                                                    >  2:547132/1‑133 (MQ=255)
                          ttttttACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGc                                                                                                                >  1:586813/1‑149 (MQ=255)
                                   cGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAtt                                                                                                       <  4:347310/149‑1 (MQ=255)
                                      ccAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCaa                                                                                                    <  3:281616/149‑1 (MQ=255)
                                      ccAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCaa                                                                                                    <  1:390230/149‑1 (MQ=255)
                                              ggACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                                                                         <  4:191129/120‑1 (MQ=255)
                                              ggACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                                                                         >  3:191129/1‑120 (MQ=255)
                                                 cGTTCAGTTACCGGCTTGCCTGAATAGCAATCCAACCGAACACACATATGCGCGGACAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCCATATATTGCaga                                                                                         <  4:295987/149‑1 (MQ=37)
                                                  gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCagag                                                                                        >  4:205802/1‑149 (MQ=255)
                                                        ttACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCttt                                                                                                                                 <  2:356354/102‑1 (MQ=255)
                                                        ttACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCttt                                                                                                                                 >  1:356354/1‑102 (MQ=255)
                                                          aCCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTa                                                                                <  2:92935/149‑1 (MQ=255)
                                                                               tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                                                                         <  2:384245/87‑1 (MQ=255)
                                                                               tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTc                                                                                                                         >  1:384245/1‑87 (MQ=255)
                                                                                   aCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAg                                                       >  2:157484/1‑149 (MQ=255)
                                                                                                  gcgcGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTa                                        >  2:252927/1‑149 (MQ=255)
                                                                                                    gcgGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATc                                      >  3:293190/1‑149 (MQ=255)
                                                                                                           gATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCActc                               >  4:579710/1‑149 (MQ=255)
                                                                                                               gttgACAAAGGGCGCTTTGTTCATGCTGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCAc                                 >  4:434743/1‑143 (MQ=255)
                                                                                                               gttgACAAAGGGCGCTTTGTTCATGCTGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCAc                                 <  3:434743/143‑1 (MQ=255)
                                                                                                                ttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAg                          >  2:526668/1‑149 (MQ=255)
                                                                                                                                  ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCAc                      >  1:137962/1‑135 (MQ=255)
                                                                                                                                  ttCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCAc                      <  2:137962/135‑1 (MQ=255)
                                                                                                                                     aTGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGc                                             >  2:520822/1‑109 (MQ=255)
                                                                                                                                     aTGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGc                                             <  1:520822/109‑1 (MQ=255)
                                                                                                                                     aTGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtt     <  3:205802/149‑1 (MQ=255)
                                                                                                                                       gCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAtttt   >  2:398248/1‑149 (MQ=255)
                                                                                                                                        ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttttt  <  4:394536/149‑1 (MQ=255)
                                                                                                                                        ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttttt  >  1:197812/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                              
AATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTT  >  NC_000913/770540‑770824

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