Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,396,099 G→C A413G (GCC→GGC)  nuoC ← NADH‑quinone oxidoreductase subunit C/D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,396,0990GC79.0% 42.2 / 7.9 19A413G (GCC→GGC) nuoCNADH‑quinone oxidoreductase subunit C/D
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (7/8);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCA  >  NZ_CP009273/2396013‑2396181
                                                                                      |                                                                                  
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:507338/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:589345/1‑90 (MQ=255)
   ccTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCaa                                                                              <  1:780728/90‑1 (MQ=255)
    cTTACGCACGTCGAAGTCGATCCCGGTGGCACGCAGGCCCGCGCCAGTGGCGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAAc                                                                             <  1:637559/90‑1 (MQ=255)
    cTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAAc                                                                             <  1:107402/90‑1 (MQ=255)
        cGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGccc                                                                         <  1:490173/90‑1 (MQ=255)
           aCGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTgg                                                                      >  1:217146/1‑90 (MQ=255)
             gTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGa                                                                    >  2:294041/1‑90 (MQ=255)
             gTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGa                                                                    >  1:444221/1‑90 (MQ=255)
                 aaGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGa                                                                <  1:128444/90‑1 (MQ=255)
                    tCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCt                                                             <  1:201403/90‑1 (MQ=255)
                                          ccGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAgcg                                       >  2:206440/1‑90 (MQ=255)
                                            gcgcCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGcgcc                                     >  2:334571/1‑90 (MQ=255)
                                            gcgcCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGcgcc                                     >  2:565436/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCttt                                >  2:811962/1‑90 (MQ=255)
                                                        ccccACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTaa                         <  2:407417/90‑1 (MQ=255)
                                                                  gcgcCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGAc               >  1:671369/1‑90 (MQ=255)
                                                                         tttCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCgg        >  2:85372/1‑90 (MQ=255)
                                                                            cgcgCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCAt     <  2:671369/90‑1 (MQ=255)
                                                                               gcCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCa  <  2:133896/90‑1 (MQ=255)
                                                                               gcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCa  <  2:709632/90‑1 (MQ=255)
                                                                                      |                                                                                  
GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCA  >  NZ_CP009273/2396013‑2396181

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