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,8990TG97.6% 130.2 / ‑3.7 41T88P (ACG→CCG) pntARe/Si‑specific NAD(P)(+) transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (20/20);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTT  >  NZ_CP009273/1671816‑1671985
                                                                                   |                                                                                      
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                                  >  1:957016/1‑90 (MQ=255)
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                                  >  2:947248/1‑90 (MQ=255)
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                                  >  2:409268/1‑90 (MQ=255)
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                                  <  2:395663/90‑1 (MQ=255)
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                                  >  1:939623/1‑90 (MQ=255)
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                                  >  1:1506579/1‑90 (MQ=255)
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                                  >  1:1279930/1‑90 (MQ=255)
  tCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAg                                                                                >  1:1118445/1‑90 (MQ=255)
  tCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAg                                                                                <  1:1000235/90‑1 (MQ=255)
     ccATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGAt                                                                             >  2:335005/1‑90 (MQ=255)
     ccATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGAt                                                                             <  2:1498577/90‑1 (MQ=255)
     ccATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGAt                                                                             <  2:1084245/90‑1 (MQ=255)
       aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTc                                                                           >  1:1207976/1‑90 (MQ=255)
           cGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTa                                                                       <  1:1490519/90‑1 (MQ=255)
            ggTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTaa                                                                      <  2:1507811/90‑1 (MQ=255)
                     acgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTTACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtc                                                             <  1:28997/90‑1 (MQ=255)
                     acgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtc                                                             >  1:585692/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                           >  2:466720/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                           >  2:1228103/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                           >  2:1021948/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGCGTCCCAGGATTCAGTAACGCAATTtcat                                                           >  2:1031196/1‑90 (MQ=255)
                           cGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCTAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCatcatc                                                       >  2:518162/1‑90 (MQ=255)
                           cGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCatcatc                                                       >  1:1302204/1‑90 (MQ=255)
                           cGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCatcatc                                                       >  1:841596/1‑90 (MQ=255)
                                        aTTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTg                                          <  1:409268/90‑1 (MQ=255)
                                        aTTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTg                                          <  1:1031196/90‑1 (MQ=255)
                                              tCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTc                                    <  2:1506579/90‑1 (MQ=255)
                                               ccGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCa                                   <  1:379614/90‑1 (MQ=255)
                                               ccGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCa                                   <  1:1539427/90‑1 (MQ=255)
                                                 ggATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGa                                 >  1:99712/1‑90 (MQ=255)
                                                           cAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTg                       <  2:939623/90‑1 (MQ=255)
                                                             ggCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc                     <  1:753864/90‑1 (MQ=255)
                                                                  gATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa                <  2:1367622/90‑1 (MQ=255)
                                                                   aTAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAg               <  1:1196477/90‑1 (MQ=255)
                                                                       aaCTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc           >  2:1213418/1‑90 (MQ=255)
                                                                       aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc           <  2:1412425/90‑1 (MQ=255)
                                                                          tCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTAt        <  2:86851/90‑1 (MQ=255)
                                                                           cACCAGCGGTGTTCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTAtt       >  1:1197824/1‑90 (MQ=255)
                                                                           cACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTAtt       >  1:1213193/1‑90 (MQ=255)
                                                                            aCCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTc      <  2:338995/90‑1 (MQ=255)
                                                                                gCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCtt  <  1:518162/90‑1 (MQ=255)
                                                                                   |                                                                                      
CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTT  >  NZ_CP009273/1671816‑1671985

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