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,0990GC96.0% 158.7 / ‑0.7 50A413G (GCC→GGC) nuoCNADH‑quinone oxidoreductase subunit C/D
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (27/21);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCA  >  NZ_CP009273/2396013‑2396181
                                                                                      |                                                                                  
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:856268/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:1176102/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:95634/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:483688/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:227604/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:462457/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:580064/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:737050/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTCCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  2:1457116/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGGAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                 >  1:693177/1‑90 (MQ=255)
    cTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAAc                                                                             <  1:1121416/90‑1 (MQ=255)
           aCGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTgg                                                                      >  1:404859/1‑90 (MQ=255)
             gTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGa                                                                    >  2:272247/1‑90 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                  <  1:462457/90‑1 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                  <  1:201322/90‑1 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                  <  1:580064/90‑1 (MQ=255)
                  aGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGAc                                                               <  2:1237438/90‑1 (MQ=255)
                     cGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCtt                                                            <  2:1146324/90‑1 (MQ=255)
                      gATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCttt                                                           >  1:458236/1‑90 (MQ=255)
                       aTCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTc                                                          <  2:1360744/90‑1 (MQ=255)
                       aTCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTc                                                          <  2:759090/90‑1 (MQ=255)
                       aTCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTc                                                          <  2:1426144/90‑1 (MQ=255)
                        tCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCa                                                         >  1:1128734/1‑90 (MQ=255)
                        tCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCa                                                         >  2:1014865/1‑90 (MQ=255)
                            ggTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                     <  1:1417044/90‑1 (MQ=255)
                            ggTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                     <  2:762890/90‑1 (MQ=255)
                            ggTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                     <  1:62267/90‑1 (MQ=255)
                              tAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATgg                                                   <  1:834715/90‑1 (MQ=255)
                                  aCGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGtt                                               >  1:21256/1‑90 (MQ=255)
                                     cAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTg                                            >  1:76234/1‑90 (MQ=255)
                                     cAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTg                                            >  2:443696/1‑90 (MQ=255)
                                     cAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTg                                            >  1:1286818/1‑90 (MQ=255)
                                      aGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGc                                           <  1:869145/90‑1 (MQ=255)
                                      aGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGc                                           >  2:443736/1‑90 (MQ=255)
                                          ccGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAgcg                                       >  1:46687/1‑90 (MQ=255)
                                            gcgcCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGcgcc                                     >  1:564342/1‑90 (MQ=255)
                                            gcgcCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGcgcc                                     >  2:1366380/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCttt                                >  2:984500/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCttt                                >  1:1264867/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCttt                                >  1:1518363/1‑90 (MQ=255)
                                                 aGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTGTGCAGCGCCGCttt                                >  1:1217955/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  1:992755/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  2:87066/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          <  1:1095787/90‑1 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  1:1118628/1‑90 (MQ=255)
                                                       gCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTa                          >  1:315899/1‑90 (MQ=255)
                                                        ccccACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTaa                         <  1:242150/90‑1 (MQ=255)
                                                            aCTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACg                     >  2:610369/1‑90 (MQ=255)
                                                            aCTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACg                     >  1:18482/1‑90 (MQ=255)
                                                              tCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGcc                   >  2:1106037/1‑90 (MQ=255)
                                                               ccAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCa                  <  1:1094655/90‑1 (MQ=255)
                                                                cAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAg                 <  2:227604/90‑1 (MQ=255)
                                                                       tctTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTc          >  2:796142/1‑90 (MQ=255)
                                                                       tctTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTc          >  1:777521/1‑90 (MQ=255)
                                                                         tttCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCgg        >  2:384803/1‑90 (MQ=255)
                                                                         tttCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCgg        >  2:1011933/1‑90 (MQ=255)
                                                                             gcgcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATc    >  1:285156/1‑90 (MQ=255)
                                                                              cgcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATcc   <  2:21256/90‑1 (MQ=255)
                                                                              cgcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATcc   <  2:1256882/90‑1 (MQ=255)
                                                                               gcCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCa  <  1:507936/90‑1 (MQ=255)
                                                                                      |                                                                                  
GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCA  >  NZ_CP009273/2396013‑2396181

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