Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,671,899 T→G T88P (ACG→CCG)  pntA ← Re/Si‑specific NAD(P)(+) transhydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,671,8990TG90.0% 86.7 / 2.7 30T88P (ACG→CCG) pntARe/Si‑specific NAD(P)(+) transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (13/14);  total (14/16)
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.36e-01

CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC  >  NZ_CP009273/1671816‑1671981
                                                                                   |                                                                                  
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                              >  2:672183/1‑90 (MQ=255)
       aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAGTTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCACAGGATTc                                                                       >  2:198096/1‑90 (MQ=255)
       aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGTGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTc                                                                       <  1:758312/90‑1 (MQ=255)
       aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTc                                                                       >  1:28665/1‑90 (MQ=255)
                 cgttacgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCaa                                                             <  2:66373/90‑1 (MQ=255)
                   ttacgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAAtt                                                           >  1:314509/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTtcat                                                       <  1:437078/90‑1 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                       <  2:611882/90‑1 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                       >  2:14790/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                       <  2:13108/90‑1 (MQ=255)
                        ttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcatc                                                      <  1:14790/90‑1 (MQ=255)
                           cGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCatcatc                                                   >  2:488264/1‑90 (MQ=255)
                                 tttttGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACgg                                             <  1:729047/90‑1 (MQ=255)
                                 tttttGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACgg                                             <  1:443148/90‑1 (MQ=255)
                                       cATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCAtt                                       <  1:138315/90‑1 (MQ=255)
                                        aTTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTg                                      <  1:89974/90‑1 (MQ=255)
                                         ttAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGa                                     >  2:570175/1‑90 (MQ=255)
                                          tAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAc                                    >  1:573631/1‑90 (MQ=255)
                                           aaTTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAcc                                   >  2:814096/1‑90 (MQ=255)
                                             ttCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCtt                                 <  2:60170/90‑1 (MQ=255)
                                                      cTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAt                        >  1:755313/1‑90 (MQ=255)
                                                      cTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAt                        >  2:780980/1‑90 (MQ=255)
                                                             ggCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc                 <  1:563366/90‑1 (MQ=255)
                                                             ggCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc                 <  2:755313/90‑1 (MQ=255)
                                                                  gagaaaaCTCACCAGCGGTGTCCCAGGATTCAGTACCGCAATTTCATCATCTACCGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa            <  2:522458/87‑1 (MQ=255)
                                                                  gATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa            <  2:682072/90‑1 (MQ=255)
                                                                   aTAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAg           <  1:420090/90‑1 (MQ=255)
                                                                       aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc       >  2:556363/1‑90 (MQ=255)
                                                                       aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc       >  1:731588/1‑90 (MQ=255)
                                                                            aCCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTc  >  1:297183/1‑90 (MQ=255)
                                                                                   |                                                                                  
CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC  >  NZ_CP009273/1671816‑1671981

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