Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,796,099 T→G 11.3% K317N (AAA→AAC U712_19580 ← Cytochrome d ubiquinol oxidase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,796,0990TG11.3% 124.7 / 2.0 56K317N (AAA→AACU712_19580Cytochrome d ubiquinol oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (18/30);  new base G (5/2);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATTTGCTGTTTTCAAGCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTCATTTGAGCTTTTTTCATTTTTTGTATCGAT  >  CP006881/3795964‑3796236
                                                                                                                                       |                                                                                                                                         
agacaagCTGTTTTCAAGCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCtttt                                                                                                                                        <  1:27348/133‑1 (MQ=255)
  aTTTGCTGTTTTCAAGCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACCCCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCgtttgc                                                                                                                                      >  3:330510/1‑137 (MQ=255)
  aTTTGCTGTTTTCAAGCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCGTTTcc                                                                                                                                      >  1:156647/1‑139 (MQ=255)
   tttGCTGTTTTCAAGCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCAGGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTGACTGGCGGGTAGTAGTCGCCtgttggg                                                                                                                                     >  1:344352/1‑136 (MQ=255)
            ttCAAGCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCGCTGGCGGTATGTAGTCGCCGTTTCGGTGTAttttt                                                                                                                            >  2:319543/1‑139 (MQ=255)
                 gCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTa                                                                                                                       <  5:87561/139‑1 (MQ=255)
                     tttACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACCCCTGCGCCTACCATGCTGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCGTTTCCGTATATTTTTTCGTACTc                                                                                                                    <  6:120414/138‑1 (MQ=255)
                     tttACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACt                                                                                                                     <  1:109717/137‑1 (MQ=255)
                     tttACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACt                                                                                                                     <  5:152235/137‑1 (MQ=255)
                     tttACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTc                                                                                                                    <  8:270262/138‑1 (MQ=255)
                     tttACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTc                                                                                                                    <  6:160105/138‑1 (MQ=255)
                     tttACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTc                                                                                                                    <  2:156647/138‑1 (MQ=255)
                     tttACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTc                                                                                                                    <  4:94573/138‑1 (MQ=255)
                           gcggTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGa                                                                                                             <  8:64860/139‑1 (MQ=255)
                            cggTTTAACCATAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGaa                                                                                                            <  5:139623/139‑1 (MQ=255)
                            cggTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGaa                                                                                                            <  2:292067/139‑1 (MQ=255)
                                    ccAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTc                                                                                                    <  3:201812/139‑1 (MQ=255)
                                     cAAACGCCGCCTAACGCAGCAAGACTCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCGTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCa                                                                                                   <  8:298009/139‑1 (MQ=255)
                                      aaaGGCCTCCTAACGCCGCAAGAATCCTGAGAACTCCTGCTCCTAGCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCAt                                                                                                  <  5:242818/139‑1 (MQ=255)
                                        aGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCtt                                                                                                                                          >  3:257216/1‑97 (MQ=255)
                                        aGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCtt                                                                                                                                          <  4:257216/97‑1 (MQ=255)
                                                        caAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCt                                                                                >  1:49292/1‑139 (MQ=255)
                                                        caAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCt                                                                                >  8:161686/1‑139 (MQ=255)
                                                                aTGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTACGTTTATTTTTTGGTACGCAGCCTGAAGGGTTTTGATCCCTTTGACACTGACGCGGAAtttgt                                                                        >  2:316729/1‑137 (MQ=255)
                                                                  gACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGaa                                                                                                            >  6:357336/1‑101 (MQ=255)
                                                                  gACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGaa                                                                                                            <  5:357336/101‑1 (MQ=255)
                                                                  gACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACAGGCGGAATGTAGTCGCCATTTCCGGATATTTTTTCGTAATCAGCCATCAGGGGTGTAATCCCTTTGACACGTCCGATGAAATTCTg                                                                       >  5:169415/1‑138 (MQ=255)
                                                                   acaacaCCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGAt                                                                     <  7:67603/139‑1 (MQ=255)
                                                                   acaacaCCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGAt                                                                     <  3:100121/139‑1 (MQ=255)
                                                                   acaacaCCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCGTTTACGGATATTTTTTAGTGCTAAGCCGTAAGGGTTTTCATCCCTTTTAAAAGTCCGCAGAATTTCTGa                                                                      >  4:84917/1‑138 (MQ=255)
                                                                    caacaCCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGAt                                                                     <  1:164034/138‑1 (MQ=255)
                                                                       cacCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGc                                                                 >  5:138003/1‑139 (MQ=255)
                                                                       cacCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTACGTATATTTTTTAGTACTCAGCCTGAAGGGTTTTAATCCCTTTGACACTTGCGCTGAATTTCTGATAAGc                                                                 >  8:167844/1‑139 (MQ=255)
                                                                        acCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGcc                                                                <  8:7009/139‑1 (MQ=255)
                                                                         ccTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCa                                                               <  6:138003/139‑1 (MQ=255)
                                                                            gCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAgg                                                                                                          >  8:12370/1‑93 (MQ=255)
                                                                            gCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAgg                                                                                                          <  7:12370/93‑1 (MQ=255)
                                                                            gCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCGTTTGCGTATATTTTTTGGGACTCAGCCTGGAGGGTTTTGAGACCTTTGACACTTCCGCTGAATTTGTGATAAGGCAAGt                                                            >  6:109778/1‑139 (MQ=255)
                                                                            gCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATCTAGTCGCCGTTTGCGTATATTTTTTAGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGt                                                            >  7:186170/1‑139 (MQ=255)
                                                                             cTCCTACCCTGATGCGGAAGCGCCAGAATGTCGTTTTCACTGGCGGAATGTAGGCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAg                                                             <  4:245213/137‑1 (MQ=255)
                                                                             cTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAg                                                             <  3:84917/137‑1 (MQ=255)
                                                                                  aCCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTc                                                                                                    <  2:233906/93‑1 (MQ=255)
                                                                                  aCCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTc                                                                                                    >  1:233906/1‑93 (MQ=255)
                                                                                  aCCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCGTTTACCCATCCGCGGAATTTCTGATAAGCCAAGTAGCTc                                                       >  1:97076/1‑138 (MQ=255)
                                                                                  aCCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTAAGCCTGAAGGGGTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTc                                                       >  3:332542/1‑138 (MQ=255)
                                                                                  aCCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTAGTAAACAACCTGAAGGGGTTTCATCCCTTTTACACATACGCTGAAATTGTGATAAGACAAATAGCTc                                                       >  2:680/1‑138 (MQ=255)
                                                                                  aCCATGATGCGGAAGCTCCAGAATGTCGTGTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTTCTCAGCCTGAAGGGTTTTCATCCCTGTGACACTTCCGCTGCATTTCTGCTAAGCCAAGTAtctc                                                       >  1:96246/1‑138 (MQ=255)
                                                                                  aCCATGATGAGGAAGCTCCAGAATGTCTTTTTCAGTGGCGGAATGTAGTCGCCTTTTGCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCa                                                      <  1:233280/139‑1 (MQ=255)
                                                                                          gCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTACGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATaa                                              >  6:58924/1‑139 (MQ=255)
                                                                                                    cAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACttt                                      <  5:58924/137‑1 (MQ=255)
                                                                                                            cGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCa                                                      >  1:43303/1‑113 (MQ=255)
                                                                                                            cGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCa                                                      <  2:43303/113‑1 (MQ=255)
                                                                                                            cGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTAGTACTCAGCCTGAAGGGTTTTAATCCCTTTGACAATTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTAAttt                            >  6:94480/1‑139 (MQ=255)
                                                                                                                 tCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTGATCCCTTTGCCCCTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTGATTTGCGCt                       >  3:158634/1‑139 (MQ=255)
                                                                                                                 tCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTAAGCCTGAAGGGTTTTCATGCCTTTGTCACTTCCGCTGCATTTCTGATAAGCCAAGTAGCTGAAGGCATTAGGAACTTTGATTTGATTTGAGCt                       >  3:153254/1‑139 (MQ=255)
                                                                                                                                  tagCCTTTTCCGTATATTTTTTCGTACTGAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAATGAACTTTGATTTCATTTGAGCTTTTTTCATTTTTTGTat      <  4:170251/137‑1 (MQ=255)
                                                                                                                                     ccTTTTCCGTATATTTTTTCGTACTAAGCCTGAATGGTTTTAATCCCTTTTACACTTCCGCTGAATTTCTGATAAGCCAAGTTGCTAAAGGCATAAGGAACTTTGATTTGATTTGTGCTTTTTTCATTTTTTGtttata   >  2:70712/1‑134 (MQ=255)
                                                                                                                                      cTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTCATTTGAGCTTTTTTCATTTTTTGTatcgat  <  5:92524/139‑1 (MQ=255)
                                                                                                                                      cTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGATGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTCATTTGAGCTTTTTTCATTTTTTGTatcgat  <  8:341983/139‑1 (MQ=255)
                                                                                                                                      cTTTTCCGTATATTTTTTCGTACACAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTCATTTGAGCTTTTTTCATTTTTTGTatcgat  <  4:330510/139‑1 (MQ=255)
                                                                                                                                      cTTTTCCGTATAATTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTCATTTGATCTTTTTTCATTTTTTGTatcgat  <  1:64809/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CCATTTGCTGTTTTCAAGCTTTTTACGGCGGTTTAACCAAAGGCCGCCTAACGCAGCAAGAATCATGACAACACCTGCTCCTACCATGATGCGGAAGCTCCAGAATGTCGTTTTCACTGGCGGAATGTAGTCGCCTTTTCCGTATATTTTTTCGTACTCAGCCTGAAGGGTTTTCATCCCTTTGACACTTCCGCTGAATTTCTGATAAGCCAAGTAGCTCAAGGCATAAGGAACTTTGATTTCATTTGAGCTTTTTTCATTTTTTGTATCGAT  >  CP006881/3795964‑3796236

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