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,5250GT95.5% 54.4 / NA 22A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (7/14);  minor base A (1/0);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAA  >  NC_000913_3_pae_tpiA/2679383‑2679661
                                                                                                                                              |                                                                                                                                        
gTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTAAGGGTAATAGGGGTTAGTGAACGCGGTTAAAAAGAAACGGCATCACGGTGGTGACTAATGCGGCTTTCGGCCAATcccggt                                                                                                                                    >  2:28328/1‑145 (MQ=255)
   tttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCTGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTACTGAACGCGGTTAAAAATAAACGGGATCAGGGTGGTGAATAATGCCGCTTTAGGCCAATCCGGCTTaa                                                                                                                                 >  2:476227/1‑149 (MQ=255)
      attaGCAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGTGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCg                                                                                                                              <  1:445039/149‑1 (MQ=255)
               ggCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCACACCGCATTTTAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAtt                                                                                                                     <  1:406315/149‑1 (MQ=255)
                     ttCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgcgc                                                                                                               <  1:298809/149‑1 (MQ=255)
                     ttCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCACCGTCAATTCAgcgc                                                                                                               >  1:42238/1‑149 (MQ=255)
                       cTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                             <  1:155326/149‑1 (MQ=255)
                           ttattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAg                                                                                                         >  2:145899/1‑149 (MQ=255)
                            tattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                                                                                                        <  1:129167/149‑1 (MQ=255)
                                 tgtgATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCat                                                                                                   <  1:178946/149‑1 (MQ=255)
                                               gcgcGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTTCGCGCGAGCAGACTCATATTGCGCAGCCTcc                                                                                     >  2:505320/1‑149 (MQ=255)
                                                 gcgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                                                                                   >  1:517912/1‑149 (MQ=255)
                                                   gTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGa                                                                                 <  1:145899/149‑1 (MQ=255)
                                                                           tAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAACCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCACCCCTAGTCCCGCGGACGCCAGCCACGcca                                                         >  2:237779/1‑149 (MQ=255)
                                                                                 ttggggAAGGCGGTTAAAAATAAAGGGCGTCACGGTGGTGAATGATGGCGTTTTCGGCCAATCCGGATTAATCGCTGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGTCGCCAGCATCGCCACCATCa                                                   <  2:487421/144‑1 (MQ=255)
                                                                                 ttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCa                                                   <  1:290699/149‑1 (MQ=255)
                                                                                            ggTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCt                                        <  2:502198/149‑1 (MQ=255)
                                                                                            ggTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCACTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCt                                        <  1:148722/149‑1 (MQ=255)
                                                                                               tAAAAAGAGAAGGCATGACGGTGGTGAATGATGCCGCTTTCGGGCAATTCGGATTAATCGCCGTCATTTCAGAGCGATCCGACTCATATTGCGCAACGTCCCCGGGTCCCGTGTACGCCAGCATCGCCACCATCAGGCATTTTTCTTTa                                     <  2:324271/149‑1 (MQ=255)
                                                                                                                          ttaatTCCGCTCTCGACCACTCGGGAGTAATCGCCCTCCAATCATTGCGAGCAGACTAACACTGCGTAACCTCCCCGAGTCCCGCGGGCGCCAGCATCGCCACCATCAGGCATTTTTCTTTATTCTGTTGTTGTAAATCGCCATACAGc          <  2:164620/148‑1 (MQ=255)
                                                                                                                                 cGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCCTCAGGCATTTTTCTTTATGCTGCTGTTGTAATTCGCCATACAGCGacaaca   >  1:203573/1‑149 (MQ=255)
                                                                                                                                  gCTTTCCGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTCTAAATCGCCCTACAGCGAcaacaa  <  1:476227/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                        
GTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAA  >  NC_000913_3_pae_tpiA/2679383‑2679661

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