Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,682 A→C 10.9% 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,6820AC10.9% 138.3 / 8.6 55intergenic (+71/‑776)mngB/cydAalpha‑mannosidase/cytochrome bd‑I ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (19/30);  new base C (3/3);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

ATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTT  >  NC_000913/770541‑770824
                                                                                                                                              |                                                                                                                                              
aTCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATAc                                                                                                                                               <  1:1189004/143‑1 (MQ=255)
 tCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAGAGGGCGCTTTGTTCATGCCGGATACg                                                                                                                                              >  1:892586/1‑143 (MQ=255)
  cAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGGAGCCACATATGCGCGGGCAGATTGTTGACAAAGGGGGCTTTGTTCATGCCGGATACgg                                                                                                                                             >  2:240277/1‑143 (MQ=255)
  cAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTGATGCCGGATACgg                                                                                                                                             <  2:382715/143‑1 (MQ=255)
    agcttCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCa                                                                                                                                           <  1:676289/140‑1 (MQ=255)
    aGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                                                                           <  1:10762/143‑1 (MQ=255)
    aGGTTCAAATCTATCGGAGCCTTTTTTTCCCGTCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCa                                                                                                                                           <  1:173354/143‑1 (MQ=255)
      gTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCacg                                                                                                                                         >  2:31284/1‑141 (MQ=255)
            ggttATCGGGGCCTTTTTTACCCGGCCGGTCACGGACGTTCAGTTACCGGCTTGCCGGAATAGCAATCAAACCCACGCCACATATGCGCGGCCAGATCGTTGACAAAGGGCGCTTCGTTCATGCCGTATCCGGCATGAAcaaa                                                                                                                                   <  2:981399/140‑4 (MQ=255)
               ttgcGGGGCCTTTTCTGCCCGGGCAGTCGCGGACGGTCAGTTACCGGCTTGCCTGAACAGCGATCAAGCCGCTGCCACATATGCGCGGCCAGATTGTTGACCAAGGGCGCTTTGTTCATGCCGGATACGGCATGCCCGCtttt                                                                                                                                <  2:853344/140‑2 (MQ=255)
                    gcgCCTTTTTTACCCGGCCAGTCACCGACGTTCAGTTACCGGCTTGCCTGAATATCAATCAAACCTAAGCCACATATGCTCTGCCATATTGTTGACAAAGGGCTCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCgg                                                                                                                           <  1:569090/141‑1 (MQ=255)
                     gggcTTGTTTACCCGGCCGGTCACGGACGTTCAGTTACCGGCTTGCCGGAATAGCAATCAAACCGACGCCACATGTGCGCGGCCAGATTCTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGt                                                                                                                          <  2:888806/140‑1 (MQ=255)
                            tttCCCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGCCGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaa                                                                                                                   <  1:310125/143‑1 (MQ=255)
                             ttACCCGGCTATTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACaaaa                                                                                                                  <  1:553329/143‑1 (MQ=255)
                                cccGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGACGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCa                                                                                                               <  2:621434/143‑1 (MQ=255)
                                 ccGGCCCGTCACGGTCGTTTAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAGGCCACATATGCGCGGCCAGATTGTTGACATAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                                                              <  2:546468/143‑1 (MQ=38)
                                 ccGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGACGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAg                                                                                                              <  1:252581/143‑1 (MQ=255)
                                           aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCaa                                                                                                    <  1:674969/143‑1 (MQ=255)
                                           aCGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCaa                                                                                                    <  2:929738/143‑1 (MQ=255)
                                                 gTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGAATGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCATTATTCGGTCTACAAAAGCAGGCAAATTCAATATAtt                                                                                              >  2:20612/1‑143 (MQ=255)
                                                   tCAGTTACCAGCTCTTCCTAATAGCGATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCGTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGc                                                                                            <  1:1134183/143‑1 (MQ=38)
                                                    ccgTTACCGGCTTGCCTGAACAGCAATCAACCCGAAGCCACATATGCGCGGCCAGATTGTGGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCa                                                                                           <  1:740560/141‑1 (MQ=255)
                                                    cAGTTACCTGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCa                                                                                           <  2:844684/143‑1 (MQ=255)
                                                    cAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCa                                                                                           <  2:317931/143‑1 (MQ=255)
                                                      gTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCaga                                                                                         <  1:368933/143‑1 (MQ=255)
                                                         atcGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGGCCGATTGTTGACAAAGGGCGCGTGGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                                                                      <  2:529646/141‑1 (MQ=255)
                                                           cGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCTGGCAAATTCAAGAGATTGCAGAgatga                                                                                    >  2:838650/1‑143 (MQ=255)
                                                              cTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGt                                                                                 <  1:1103997/143‑1 (MQ=255)
                                                                  gCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGa                                                                              <  1:512286/143‑1 (MQ=255)
                                                                    cTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCAt                                                                                                                                          >  1:439107/1‑81 (MQ=255)
                                                                    cTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGTATCCGGCAt                                                                                                                                          <  2:439107/81‑1 (MQ=255)
                                                                      gAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACt                                                                          >  2:885167/1‑143 (MQ=255)
                                                                              aTCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATCCAATATATTGCAGAGAGGATGTAGACACTGGCAAGCg                                                                  >  1:1158381/1‑143 (MQ=255)
                                                                               tCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGCTACGGCATGAACGCTTTATTCGGCCTACAAAAGCAGGCAAACTCAATATATTGCAGAGCTGATGTAGACACTGGCAACCGt                                                                 >  2:389187/1‑143 (MQ=38)
                                                                                  aaCCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATTTATTGCAGAGATGATGTAGACACTGGCAAGCGTAgc                                                              <  1:824500/143‑1 (MQ=255)
                                                                                       aaGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                                                                      >  2:1140439/1‑114 (MQ=255)
                                                                                       aaGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAgat                                                                                      <  1:1140439/114‑1 (MQ=255)
                                                                                          ccACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAgg                                                      >  1:978519/1‑143 (MQ=255)
                                                                                            acaTATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACACAAGCAGGCACATTCAATCTATTGCAGCGATGATGTAGACACTGGCAACCCTACCGCCTCAGGTa                                                    >  2:426940/1‑143 (MQ=25)
                                                                                                 tGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAAtttt                                               <  2:935789/143‑1 (MQ=255)
                                                                                                  gcgcGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTg                                              >  2:944471/1‑143 (MQ=255)
                                                                                                             ttgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTATGCAAAAACACGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTc                                   <  1:933630/143‑1 (MQ=255)
                                                                                                              tgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCa                                  <  1:1257974/143‑1 (MQ=255)
                                                                                                              tgttgACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGAAGAGAGGATGTAGACACTGGCAAGCGTAGGGCATGAGGTAATTTAGCGTTTATCTTCa                                  >  2:631661/1‑143 (MQ=25)
                                                                                                                  gACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCACATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTCTCTTCActct                              >  2:502817/1‑143 (MQ=255)
                                                                                                                   gcAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTAGTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACtctc                             <  2:108895/142‑1 (MQ=255)
                                                                                                                   aCAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCACTATATTGCAGCGCTCATGGAGACACTGGCAAGCGTAGCGCATCACGCAATTTTGCGTTTATCTTCACtctc                             >  2:147226/1‑143 (MQ=255)
                                                                                                                   aCAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACtctc                             >  2:1010841/1‑143 (MQ=255)
                                                                                                                     aaaGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCaa                           <  1:178851/143‑1 (MQ=255)
                                                                                                                          gcgcTTTGTTCATGCCGGATACGGGATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCAc                      <  2:1062999/143‑1 (MQ=255)
                                                                                                                                  ttCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTTGCGCATCACGTGATTTTGGGTTTATCTTCACTCTCAAGCCACGTTTAtgt              >  1:650517/1‑143 (MQ=37)
                                                                                                                                  ttCATGCCGTATCCGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATAtgt              >  1:810048/1‑143 (MQ=255)
                                                                                                                                       gCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCtt         >  2:34356/1‑143 (MQ=255)
                                                                                                                                        ccGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCttt        >  1:260309/1‑143 (MQ=255)
                                                                                                                                          ggATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATTTTCACTCTCAAGCCACGTATATGTGGCTTTAt      >  1:498284/1‑143 (MQ=255)
                                                                                                                                              aCGTCATGGTCGCTTTATTCGGTCTACAAGAGCAGGCAAAGTCAATATCTTGCGGTGATGATGTAGACGCTGGCAAGCGTAGTGCATCATGTACTTTTGCGTTTATCTTCACTCTCTAGCCACGTATATGTGGCTTTAttttt  <  2:1007944/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
ATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTATTTTT  >  NC_000913/770541‑770824

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