Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,524,942 T→G 12.2% G112G (GGT→GGG)  U712_07850 → Cytochrome c oxidase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,524,9420TG12.2% 115.3 / 5.5 60G112G (GGT→GGG) U712_07850Cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base G (5/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTAATGACTATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGCA  >  CP006881/1524810‑1525055
                                                                                                                                    |                                                                                                                 
aGTAATGACTATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTGCATTTTTAAACGCGCTTGGTTTTTgg                                                                                                             >  3:242331/1‑139 (MQ=255)
aGTAATGACTATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTGCAGTTTTTAACGCGCTTGGTTTTTgg                                                                                                             >  3:208385/1‑139 (MQ=255)
         tATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTACACGCGCTGGGTTTTTGGCGATTCttt                                                                                                    >  2:116515/1‑139 (MQ=255)
         tATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCttt                                                                                                    >  6:137849/1‑139 (MQ=255)
         tATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCttt                                                                                                    >  3:186511/1‑139 (MQ=255)
         tATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTAGTTTCCATTTTTTAACGCGCTTGGGTTTTGGCTCTTCttt                                                                                                    >  3:303234/1‑139 (MQ=255)
           tGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGATGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCttttt                                                                                                  >  7:160519/1‑139 (MQ=255)
           tGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCttttt                                                                                                  >  8:19439/1‑139 (MQ=255)
           tGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGGTTTTGGCTCTTCttttt                                                                                                  >  4:231889/1‑139 (MQ=255)
           tGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCGTGGGTTTTGGCTCTTCttttt                                                                                                  >  7:348458/1‑139 (MQ=255)
            gCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTTGTTTGCATTTTTTAACGCGCTTGGTTTTTGGCtcttggttt                                                                                                   >  3:48541/1‑132 (MQ=255)
            gCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAAGGCGCTTGGTTTTTGGCTCTTCTTTTTc                                                                                                 >  2:38385/1‑139 (MQ=255)
                ggAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGctc                                                                                                          <  6:292871/126‑1 (MQ=255)
                ggAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAAAGCGCTTGGTTTTTGGctc                                                                                                          >  5:292871/1‑126 (MQ=255)
                 gAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTAGGTTGCATTTTTAAACGCGCTTGGTTTTTGGCTCTTATTTTTCtgggg                                                                                            >  4:90656/1‑139 (MQ=255)
                  aaCGCCGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCAGTTTTAAACGCGCGTGGTTTTTGGCTCTTCTTTTTCggggg                                                                                            >  4:281778/1‑138 (MQ=255)
                  aaCGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCggggg                                                                                            >  6:175641/1‑138 (MQ=255)
                  aaCGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCggggg                                                                                            >  3:88789/1‑138 (MQ=255)
                    cgacgaTGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCAt                                                                                         <  4:88789/139‑1 (MQ=255)
                       cgaTGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCAtttt                                                                                      >  2:220876/1‑139 (MQ=255)
                         atgatCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCCTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCt                                                                                    <  7:232915/139‑1 (MQ=255)
                            atCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTGTCCATTTTTAAACTCGCTTGGGTTTTGGCTGGTCTTTTTGGGGGGGAttttggtt                                                                                  >  6:217963/1‑134 (MQ=255)
                             tCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGTATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTaa                                                                                 <  4:331752/138‑1 (MQ=255)
                                 tCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCt                                                                             <  7:19439/138‑1 (MQ=255)
                                    ggCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGc                                                                         <  2:271206/139‑1 (MQ=255)
                                     gCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGGTTCCATTTTTAAACGCGCGTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCt                                                                        >  8:167370/1‑139 (MQ=255)
                                           aTGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTGGTTTGCATTTTTATACGGGGTTGGTTTTTGGGTGATGGTTTTGGGGGGGCTTTTGGTAATTATAAGGTGGgttt                                                                  >  2:219049/1‑139 (MQ=255)
                                                   ctTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAAAGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCg                                                            >  7:209684/1‑137 (MQ=255)
                                                       tttGCCTTGATGAATGCGGTTGTGCCTTTGCAGAGCGGGGCGCGGGACGCTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAgcg                                                       <  4:357069/138‑1 (MQ=255)
                                                       tttGCCTTGATGAATGCGGGTGGGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAgcgc                                                      <  5:111107/139‑1 (MQ=255)
                                                           ccTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTGTCGGTTCCCTTTTTAAACGCGCGTGGGGGTTGGGGGTTCGTTTTGGGGGGGGTGTTCTTAAATCTGAGCTGGTTTTTGGTGGgcgcgcgcgc                                                  >  1:347045/1‑135 (MQ=255)
                                                            cTTGATGAATGCGTTTGTGCCTTTGCAGATCGGGGCGCGGGAGGTTTCGTTTCCATTTTTAAAGGCGCTTGGTTTTTGTCTCTTCTTTTTCGGGGGCATTTTCTTAAATTTAAGCTGGTTTTTGGGCTGAGCGCCAGAt                                                 <  2:329258/139‑1 (MQ=255)
                                                            cTTGATGAATGCGGTGGTGCCTTTGCAGATCGGGTCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGGTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGAt                                                 <  3:281778/139‑1 (MQ=255)
                                                              tgatgaATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCGTAGTTTTCGGGGGCAGTTTCTGAAATCGAAGCTGGTTTTTGGGAGGAGCGCCAGATGc                                               >  5:306999/1‑139 (MQ=255)
                                                                 tgaATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCggg                                            <  5:137849/139‑1 (MQ=255)
                                                                 tgaATGCGGTTGTGCCTTTGCAGATCGGGGCGCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCggg                                            <  2:13977/139‑1 (MQ=255)
                                                                     tGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATg                                                <  8:60392/131‑1 (MQ=255)
                                                                     tGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATg                                                >  7:60392/1‑131 (MQ=255)
                                                                     tGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATATAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTgg                                        >  6:336393/1‑139 (MQ=255)
                                                                     tGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATATTTAAAAGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAAACTAAGCTGGTTTTTGGGCGGAGGGCAATATCCGGGCTgg                                        >  1:124704/1‑139 (MQ=255)
                                                                       cGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATCTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGa                                       <  7:292032/138‑1 (MQ=255)
                                                                          ttGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAAAGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGGATTTTCTTAATTCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGAcat                                    >  3:262952/1‑138 (MQ=255)
                                                                                   tGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCt                           <  5:336393/138‑1 (MQ=255)
                                                                                            gggCTCGGGATGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCt                 <  8:160519/139‑1 (MQ=255)
                                                                                            gggCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCt                 >  7:91188/1‑139 (MQ=255)
                                                                                            gggCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCt                 <  4:186511/139‑1 (MQ=255)
                                                                                            gggCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCt                 <  3:90656/139‑1 (MQ=255)
                                                                                                         ttcgtttcCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCt                 >  7:67704/1‑126 (MQ=255)
                                                                                                         ttcgtttcCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCt                 <  8:67704/126‑1 (MQ=255)
                                                                                                         ttcgtttcCATTTTTAAACGCGCGTGGTTTTTGGCTGTTCTTTTTCGGGGGCATTTTCTTTAATCTAAGCTGGGTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCGCCCTTTCGCTTCATTGTAAAtgg    >  4:222691/1‑139 (MQ=255)
                                                                                                           cgtttcCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGc   <  1:220876/138‑1 (MQ=255)
                                                                                                           cgtttcCATTTTTAAACGCGCTTGGGTTTTGGCTCTGCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGCa  <  8:41268/139‑1 (MQ=255)
                                                                                                           cgcttcCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGc   <  1:116515/135‑1 (MQ=255)
                                                                                                            gtttcCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGCa  <  3:231889/138‑1 (MQ=255)
                                                                                                            gtttcCATTTTTAAACGCGCTTGGTGTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTACGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGCa  <  4:242331/138‑1 (MQ=255)
                                                                                                            gtttcCATTTTTAAACGCGATTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTTGTGTTTGGGCTGAGCGCCTGTTGCGGGCTGGACATCATATGCCTCACTTTCGGTTCATTCTAAAGGGCa  <  5:296327/138‑1 (MQ=255)
                                                                                                            gtctcCATTTCTACACGAGCTTGGGTTATGGCTCTTCATATTCGGGGGCATTTTCTTAAATCTCCGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCGGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGCa  <  7:265236/135‑1 (MQ=255)
                                                                                                            ggttcCTTTTTTAACCGGGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGCa  <  8:209684/136‑1 (MQ=255)
                                                                                                                            gcgcTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCgg                                                           <  8:11005/65‑1 (MQ=255)
                                                                                                                            gcgcTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCgg                                                           >  7:11005/1‑65 (MQ=255)
                                                                                                                                    |                                                                                                                 
AGTAATGACTATGCACGGAACGACGATGATCTTTCTGGCTGCGATGCCGCTCTTATTTGCCTTGATGAATGCGGTTGTGCCTTTGCAGATCGGGGCTCGGGACGTTTCGTTTCCATTTTTAAACGCGCTTGGTTTTTGGCTCTTCTTTTTCGGGGGCATTTTCTTAAATCTAAGCTGGTTTTTGGGCGGAGCGCCAGATGCGGGCTGGACATCATATGCCTCACTTTCGCTTCATTCTAAAGGGCA  >  CP006881/1524810‑1525055

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: ‑