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,8990TG94.4% 51.9 / ‑3.7 18T88P (ACG→CCG) pntARe/Si‑specific NAD(P)(+) transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (3/14);  total (4/14)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTAT  >  NZ_CP009273/1671821‑1671979
                                                                              |                                                                                
ccATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGAt                                                                       <  2:615424/90‑1 (MQ=255)
ccATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGAt                                                                       <  2:152457/90‑1 (MQ=255)
ccATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGAt                                                                       <  2:330315/90‑1 (MQ=255)
        gTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAAc                                                               <  2:163173/90‑1 (MQ=255)
                acgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtc                                                       <  2:600840/90‑1 (MQ=255)
                  gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                     <  2:599436/90‑1 (MQ=255)
                      cGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCatcatc                                                 >  2:203463/1‑90 (MQ=255)
                            tttttGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACgg                                           <  2:93967/90‑1 (MQ=255)
                                   aTTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTg                                    <  1:7161/90‑1 (MQ=255)
                                                   gcgcAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAtct                    >  2:134705/1‑90 (MQ=255)
                                                        ggCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc               >  1:122079/1‑90 (MQ=255)
                                                             gATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa          <  2:248327/90‑1 (MQ=255)
                                                              aTAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAg         <  2:215281/90‑1 (MQ=255)
                                                                aaaaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGAc       <  1:287249/90‑1 (MQ=255)
                                                                  aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc               <  1:708904/80‑1 (MQ=255)
                                                                  aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc               >  2:708904/1‑80 (MQ=255)
                                                                     tCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTAt  <  2:352516/90‑1 (MQ=255)
                                                                     tCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTAt  <  1:666937/90‑1 (MQ=255)
                                                                              |                                                                                
CCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTAT  >  NZ_CP009273/1671821‑1671979

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