Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 2,679,525 G→T A1074D (GCT→GAT)  yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,679,5250GT97.1% 98.2 / ‑5.9 35A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (20/14);  total (21/14)
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.41e-01

GGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGT  >  NC_000913_3_pae_tpiA/2679382‑2679663
                                                                                                                                               |                                                                                                                                          
ggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGa                                                                                                                                       <  1:194624/149‑1 (MQ=255)
ggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGa                                                                                                                                       <  1:231457/149‑1 (MQ=255)
  tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGGCAATCCCGAta                                                                                                                                     >  2:591391/1‑148 (MQ=255)
  tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGAtt                                                                                                                                     >  1:3517/1‑149 (MQ=255)
       attaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCg                                                                                                                                >  1:227391/1‑149 (MQ=255)
                       tCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGcgcg                                                                                                                >  1:378689/1‑149 (MQ=255)
                          aaTTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAtt                                                                                                                       >  2:474511/1‑139 (MQ=255)
                          aaTTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAtt                                                                                                                       <  1:474511/139‑1 (MQ=255)
                           attattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATc                                                                                                                                 <  1:325063/128‑1 (MQ=255)
                           attattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATc                                                                                                                                 >  2:325063/1‑128 (MQ=255)
                           attattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTAGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAGCGCCGTCAATTCAGCGCGAGCa                                                                                                            >  2:531012/1‑149 (MQ=255)
                                tttGTGATGAGCGGCAGCGCGTGCCGCTCAAATATCACAGTTAATAGGGCTTAGTGAACGCGGTTAAAAAT‑AAAGGCATCACGGCGGTGCATAATGCCGCTTTCGACCCATCCCGCTTCATCGCCGTCAATCCAGCGCCACCACACTcc                                                                                                      >  2:564610/1‑148 (MQ=255)
                                       tgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTCATCGCGGTCAATTCAGCGCGAGCAGACTCAGAGTgc                                                                                                >  2:45934/1‑149 (MQ=255)
                                            ggCAGCGCGTGCTGCTCAAATATTACAGGTAATAGGGTTTAGTGAACGCTGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                                                                                           <  2:499737/149‑1 (MQ=255)
                                            ggCAGCGCGTGCCGCTCAAATATTACAGTAAATAGGGTTTAGTGGACGCGGTTAAAAATAAGCGGCATCACGGTGGTGAATAATACACCTTTCCGCCAATCAGGCTTACTCGCCGTCCATTAAGCGCCACCACACTCATATTacgccct                                                                                           >  2:202272/1‑146 (MQ=255)
                                                        cGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTcc                                                                               >  1:180261/1‑149 (MQ=255)
                                                        cGCTCAAATATTACAGTTAATAGGGTGTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTcc                                                                               >  2:565002/1‑149 (MQ=255)
                                                         gCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTccc                                                                              <  1:139548/149‑1 (MQ=255)
                                                           tCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCcgc                                                                            <  1:24154/149‑1 (MQ=255)
                                                                tatTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGAAACCTCCCCGAGCCCCGCGGACg                                                                       >  2:478192/1‑149 (MQ=255)
                                                                    aCAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                               >  1:626503/1‑105 (MQ=255)
                                                                    aCAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                               <  2:626503/105‑1 (MQ=255)
                                                                    aCAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTAATATTGCGCAACCTCCCCGAGCCACGCGGACGcccg                                                                   >  1:416133/1‑147 (MQ=255)
                                                                                gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCCAGTCCCGCGGACGCCAGCCTCGCCACCAt                                                       >  2:373333/1‑149 (MQ=255)
                                                                                            cgGTTAAAAATAAAAGGCATCACGGTGGTGAATTATGCCGCTTTCGGCTAATCCGGATTAATCGCCGTCAATTCAGCGCGCGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTc                                           <  1:535902/149‑1 (MQ=255)
                                                                                                tAAAAATAAAATGCATCTCGGTGGTGAATAATGCCGCTTTCGGCTAATCCGAATTAATCGCCGTCAATTCAGCGCGATTAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTa                                       <  2:386112/149‑1 (MQ=255)
                                                                                                   aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCTGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGGCACCATCAGGCATTTTTCTTTAtgc                                    >  1:575947/1‑149 (MQ=255)
                                                                                                        aaaGGCATCACGGTGGTGAATAATGCCGCATTAGTTAAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGctgctg                               >  1:156243/1‑149 (MQ=255)
                                                                                                            gCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgttgt                           <  1:629291/149‑1 (MQ=255)
                                                                                                                 aCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTATAAATc                      >  2:385756/1‑149 (MQ=255)
                                                                                                                      ggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCAt                 <  1:607052/149‑1 (MQ=255)
                                                                                                                         gAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGccacca                                                        <  1:583491/107‑1 (MQ=255)
                                                                                                                         gAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGccacca                                                        >  2:583491/1‑107 (MQ=255)
                                                                                                                         gAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACa              >  1:183353/1‑149 (MQ=255)
                                                                                                                                     tttCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGt  <  1:66247/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                          
GGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGT  >  NC_000913_3_pae_tpiA/2679382‑2679663

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