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.4% 109.1 / NA 38A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (20/17);  minor base C (1/0);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGG  >  NC_000913_3_pae_tpiA/2679379‑2679667
                                                                                                                                                   |                                                                                                                                              
tAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                                                                                                              <  1:640467/149‑1 (MQ=255)
   ggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGa                                                                                                                                           <  1:399853/149‑1 (MQ=255)
        ttattaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                                                                                                      >  1:330190/1‑149 (MQ=255)
            taACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGcc                                                                                                                                  <  1:304986/149‑1 (MQ=255)
                 acGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCaa                                                                                                                             >  1:170267/1‑149 (MQ=255)
                             aaTTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGc                                                                                                                 >  1:548765/1‑149 (MQ=255)
                                        tttgaGCGGCA‑CGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATAtt                                                                                                      >  1:39168/3‑149 (MQ=255)
                                        gatgaGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATCGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATACGGATTAATCCCCGTCAATTCATCGCGAGCAGACTCATAtt                                                                                                      >  2:31749/1‑149 (MQ=255)
                                        gatgaGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAACCCGGATTAATCGCCGTCAATTCAGCGCGAGCAAACTCATAtt                                                                                                      >  2:357999/1‑149 (MQ=255)
                                           gaGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATGATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcg                                                                                                   >  2:402489/1‑149 (MQ=255)
                                                           ccGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTCCGCAACCTCCCCGATTc                                                                                    >  2:480539/1‑149 (MQ=255)
                                                                 aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCATCAATTCAGCGCGAGCAGACTCATATCGCGCAACCGCCCCGCGTTCcgccg                                                                              >  2:46050/1‑147 (MQ=255)
                                                                       tACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGCATTAATCCCCGTCAATTCAGCGCGAGCAGACTCATATTCCGCAACCTCCCCCAGCCCCCCTCACGCCa                                                                        >  2:413575/1‑149 (MQ=255)
                                                                         cAGTTAATGGGGTTTAGTGAACGCGGTGAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGc                                                                      <  2:62452/149‑1 (MQ=255)
                                                                          aGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCa                                                                     >  1:13011/1‑149 (MQ=255)
                                                                             tAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCg                                                                  <  1:128013/149‑1 (MQ=255)
                                                                                    gTTTAGCGAACTCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCAt                                                           <  1:34494/149‑1 (MQ=255)
                                                                                     tttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCACCGCCACCATc                                                          >  1:442244/1‑149 (MQ=255)
                                                                                      ttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCa                                                         <  1:386778/149‑1 (MQ=255)
                                                                                       tAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAg                                                        <  2:170267/149‑1 (MQ=255)
                                                                                         gTGAACGAGATTAAAAATAAAACGTATCACGTTGGTAAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGc                                                      <  1:1234/149‑1 (MQ=255)
                                                                                          tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTAATATTGCGCAACCTCCCCGAGTCCCGCGGACGGCAGCATCGCCACAATCACGCa                                                     >  2:470906/1‑149 (MQ=255)
                                                                                                 ggTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCt                                              <  1:443680/149‑1 (MQ=255)
                                                                                                 agTTAAAAATAAAGGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCt                                              <  2:548284/148‑1 (MQ=255)
                                                                                                   ttAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCttt                                            >  1:223649/1‑149 (MQ=255)
                                                                                                                gCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGGGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCCGGCATTTTTCTTTATGCTGCtgttgt                               >  1:114320/1‑149 (MQ=255)
                                                                                                                   tCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTaaa                            >  1:386090/1‑149 (MQ=255)
                                                                                                                   tCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTaaa                            >  1:224227/1‑149 (MQ=255)
                                                                                                                      cGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGTGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCg                         <  2:39168/149‑1 (MQ=255)
                                                                                                                         tggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCa                      <  2:365308/149‑1 (MQ=255)
                                                                                                                               ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGc                <  2:114320/149‑1 (MQ=255)
                                                                                                                                 aatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGa              >  1:277625/1‑149 (MQ=255)
                                                                                                                                      cGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGcc                                                                         <  1:487632/85‑1 (MQ=255)
                                                                                                                                      cGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGcc                                                                         >  2:487632/1‑85 (MQ=255)
                                                                                                                                      cGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGacaaca         >  1:578673/1‑149 (MQ=255)
                                                                                                                                       gCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGAcaacaa        <  1:402489/149‑1 (MQ=255)
                                                                                                                                           tCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCa    <  1:216590/149‑1 (MQ=255)
                                                                                                                                             ggCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAgg  >  1:189635/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCG‑TGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGG  >  NC_000913_3_pae_tpiA/2679379‑2679667

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: ‑