Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,682 A→C 12.0% intergenic (+71/‑776) mngB → / → cydA alpha‑mannosidase/cytochrome bd‑I ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913770,6820AC12.0% 80.1 / 4.6 39intergenic (+71/‑776)mngB/cydAalpha‑mannosidase/cytochrome bd‑I ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base C (2/3);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

AATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTT  >  NC_000913/770540‑770822
                                                                                                                                              |                                                                                                                                             
aaTCAAGGTTCAAATCTACCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTCACCGGGTGGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATa                                                                                                                                               <  1:241892/143‑1 (MQ=255)
  tCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACg                                                                                                                                             >  1:265675/1‑143 (MQ=255)
    aaGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGc                                                                                                                                           <  2:491046/143‑1 (MQ=255)
       gTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATg                                                                                                                                        >  2:919271/1‑143 (MQ=255)
        ttCAAATCTATCGGGGCCTTTTCTACCCTGCCGGTGCTGTGTGCGCACTGACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGa                                                                                                                                       <  1:235991/143‑1 (MQ=255)
             cccTATCCGGGACTTTTTGGCCCGGCGTGTCTCGGACGTTCGGTCACCGGCTTGCCTGAAGAGCAATCAAACCGAAGCCACGTATGCGCGGCCAGATTGGTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAcaaa                                                                                                                                  <  2:410183/141‑4 (MQ=255)
              tCTATCGGGACCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAATCCGAAGCCACATATGCGCCGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCtt                                                                                                                                 <  1:637098/143‑1 (MQ=255)
                      ggCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                                                                         <  2:934064/143‑1 (MQ=255)
                                  ccGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                                                            <  1:852267/143‑1 (MQ=255)
                                   cgcccAGTCACGGCCTTTCATTTACCCGCTTGCCTGAATAGCAAGCAAACCGCAGCCACAGAGGCGCGGCCAGATTGTTGACAGAGGGCGCTTGGTTCATGCCGGAAACGGCATGAACGCTTTATTCGGTCTACAAAAGCAgg                                                                                                           <  2:871902/140‑1 (MQ=255)
                                            aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCaa                                                                                                  <  2:802624/143‑1 (MQ=255)
                                             cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCCGCCCAATTCact                                                                                                 >  1:882237/1‑141 (MQ=38)
                                             cGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAat                                                                                                 >  1:359183/1‑143 (MQ=255)
                                              ggACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGCTACGGCCTGAACGCTTTTTTTGGTCTTCAATAGCAGGCAAATTCAata                                                                                                >  1:331945/1‑143 (MQ=255)
                                                   ttCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTg                                                                                           <  1:316655/143‑1 (MQ=255)
                                                         tACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAACTTCTATATATTGCAgaga                                                                                     >  1:36884/1‑143 (MQ=255)
                                                                ttGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTa                                                                              <  2:72976/143‑1 (MQ=255)
                                                                       aaTAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAACTCAATATATTGCAGAGATGATGTAGACACTg                                                                       >  2:494596/1‑143 (MQ=255)
                                                                         tAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGc                                                                     >  1:913692/1‑143 (MQ=255)
                                                                           gCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCaa                                                                   <  1:919271/143‑1 (MQ=255)
                                                                            cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAg                                                                  >  1:595336/1‑143 (MQ=255)
                                                                            cAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAg                                                                  >  1:242432/1‑143 (MQ=255)
                                                                             aaTCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTATATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGc                                                                 >  2:587536/1‑143 (MQ=255)
                                                                               tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAcaaa                                                                                                                                  >  1:477848/1‑74 (MQ=37)
                                                                               tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCATGAAcaaa                                                                                                                                  <  2:477848/77‑4 (MQ=37)
                                                                                     cGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCCGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAAGCCAATACATTGCAGAGGTGATGTAGACACTGGCAACCCTAGCGCa                                                         >  2:330837/1‑143 (MQ=255)
                                                                                         gCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCCTTATTCGGTCTACACAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAg                                                     <  1:164812/143‑1 (MQ=255)
                                                                                              ttatGCGCGGCCAGATTGTTGACAGCGGGCGCTTGGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAGCATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAt                                                <  2:38888/142‑1 (MQ=255)
                                                                                                              tgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCa                                <  2:514568/143‑1 (MQ=255)
                                                                                                               gttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCAc                               <  2:265675/143‑1 (MQ=255)
                                                                                                                   aCAAAGGGCGCTTTGTTCATGCCGTATACGGCATGAACGCTTAATTGAACAATGGAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACtctc                           <  1:996725/143‑1 (MQ=37)
                                                                                                                        gggCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGcc                      <  1:45518/143‑1 (MQ=255)
                                                                                                                          gcgcTTGGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCAc                    <  1:364965/143‑1 (MQ=255)
                                                                                                                            gcggtGTTCATGCCGGAGACGGCATGAACGCGTTATGCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGt                  <  2:1039290/139‑1 (MQ=255)
                                                                                                                                    cATGCCGGATACGGCATGAACGCTCTATTCGGCCTACAAAAGCAGGCAAACTCAATATAGCGCAGAGAAGATGTAGACACTCGCAAACGTCGCGCATCAGGTAATTTTTCGTTTCTCTTCACTCTCCAGCCACGTATAtgtgc          >  1:1028415/1‑142 (MQ=255)
                                                                                                                                      tGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCt        >  2:928349/1‑143 (MQ=255)
                                                                                                                                        ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGGGCATCAGGTAATGTTGGGTTTATGTTCACTCTCAAGCCACGGATATGTGGCttt      >  1:8892/1‑143 (MQ=38)
                                                                                                                                         cgtaTACGGCATGAACGCTTTACTCGGTCTACAAAAGCAGGCAAGTTCAGTAGATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACATATACGTGGCttga     <  2:646683/140‑3 (MQ=255)
                                                                                                                                           gATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATGCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCATCACTCTCAAGGCACATATACGTGGCttgaga   <  2:550172/143‑5 (MQ=255)
                                                                                                                                            aTACGGGAGGAACGGTGTATTCGGTCTACAATAGCAGGCAAATTCAATATTTTGCAGAGATGATGTAGACACTGGTAAGCGTAGCGCATCTGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAttt  <  2:488287/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
AATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTT  >  NC_000913/770540‑770822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: