Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 770,679 G→T 11.0% intergenic (+68/‑779) 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,6790GT11.0% 131.1 / 9.9 55intergenic (+68/‑779)mngB/cydAalpha‑mannosidase/cytochrome bd‑I ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (19/30);  new base T (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 = 3.20e-01

GAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAT  >  NC_000913/770537‑770820
                                                                                                                                               |                                                                                                                                             
gAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCGTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCgg                                                                                                                                               >  2:824500/1‑143 (MQ=255)
 aaaaTCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGCTTGTTGACAAAGGGCGCTTTGTTCATGTCGGa                                                                                                                                              >  1:509398/1‑143 (MQ=255)
    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)
                                                                                                                                               |                                                                                                                                             
GAAAATCAAGGTTCAAATCTATCGGGGCCTTTTTTACCCGGCCAGTCACGGACGTTCAGTTACCGGCTTGCCTGAATAGCAATCAAACCGAAGCCACATATGCGCGGCCAGATTGTTGACAAAGGGCGCTTTGTTCATGCCGGATACGGCATGAACGCTTTATTCGGTCTACAAAAGCAGGCAAATTCAATATATTGCAGAGATGATGTAGACACTGGCAAGCGTAGCGCATCAGGTAATTTTGCGTTTATCTTCACTCTCAAGCCACGTATATGTGGCTTTAT  >  NC_000913/770537‑770820

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: