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,5250GT96.0% 70.7 / NA 27A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (13/11);  total (16/11)
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

TTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGG  >  NC_000913_3_pae_tpiA/2679385‑2679667
                                                                                                                                            |                                                                                                                                               
ttttATTAACAACGGCTTATTCTAATTATATTGTGATGAGCGGCAGCGCATGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGAGAAAAATCCCGCCTTCGGCCAACCCCGATTa                                                                                                                                         >  2:23108/1‑149 (MQ=255)
          aacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCa                                                                                                                              >  2:353908/1‑149 (MQ=255)
                   ttCTTATTAGTTTGTTATGTTCTGCAGCGAGTGCCGCTCATATATTACAGTTAATAGGGTTTAGTGAATGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgcgc                                                                                                                     <  2:168597/149‑1 (MQ=255)
                     cTAATGATTGTGTTAGTAGCGTCAGCGCGTGCCGCACAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATCATGCAGTTAGAGGCCAAACCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                                   >  2:366544/1‑149 (MQ=255)
                      tAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAg                                                                                                                  >  2:396966/1‑149 (MQ=255)
                             tttGTGATGAGTGGTAGCGTGTTATTCGCTGATATTATATTTAAAATGGTTTAGTTAACGCGTTTAAAAATAAAAGGCTTCTCGGTGGTGAAGAATGCCGCTTTCGGCCAATCCGGATTATTCGCCGTCAATTCAGCGCGAGCAGAcgc                                                                                                           <  2:492032/149‑3 (MQ=255)
                                    tgaGCGGAAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTGAAAAATAAATGGCATCACGATGGTGAATAATTACGCTTTCGGCCAAGCGGGATACATCGCCTGCGATGTCGCGGGAGCGGCCTCATGTTgc                                                                                                    >  2:86923/1‑149 (MQ=255)
                                                      gCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTCAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTccc                                                                                  >  1:448791/1‑149 (MQ=255)
                                                                   aGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCACCCTCCCCGAGTCCCGCGGACGCCAgac                                                                     >  2:275689/1‑147 (MQ=255)
                                                                     ttAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATc                                                                   <  1:300510/149‑1 (MQ=255)
                                                                     ttAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATc                                                                   >  1:332003/1‑149 (MQ=255)
                                                                        aTAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGcc                                                                >  1:319809/1‑149 (MQ=255)
                                                                        aTAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGGGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAATCCCGCGGACGCCAGCATCgac                                                                >  2:393887/1‑147 (MQ=255)
                                                                            ggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGccacca                                                            >  2:198683/1‑149 (MQ=255)
                                                                                     aatccgGTTAAAAAAAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTAAATTCAGCGCGAGCAGACTCATATTGCGCACCCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAtt                                                   <  2:165180/145‑1 (MQ=255)
                                                                                                        ggCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgttg                                <  1:396966/149‑1 (MQ=255)
                                                                                                         gCATCACGGTGGTTAATAATGCCGCTGTCGGCCAATCCGGACTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCGGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgttgt                               <  1:493139/149‑1 (MQ=255)
                                                                                                          cATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTa                              >  2:203665/1‑149 (MQ=255)
                                                                                                           aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTaa                             <  1:251735/149‑1 (MQ=255)
                                                                                                              aCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCg                                                                                                                                     <  1:515948/43‑1 (MQ=255)
                                                                                                              aCGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCg                                                                                                                                     >  2:515948/1‑43 (MQ=255)
                                                                                                                 gtggtgAATCATGCCGGTTTCGGCCAATCCGGCTTACTCCATGCTCCATCAGCTCGAGTCGACTCAACTTTCGTTACTCCCCCAAAATGCCCCGACGTCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGcc                       <  1:122650/149‑1 (MQ=255)
                                                                                                                 gtggtgAATAATGCCGCTTTCGGTCAATCCGGATTAATCGCCGTCAATTCAGTGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGcc                       <  2:304783/149‑1 (MQ=255)
                                                                                                                           atGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTCTAAATCGCCATTCAGCGac             >  1:165320/1‑149 (MQ=255)
                                                                                                                            tGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGaca            >  1:214757/1‑149 (MQ=255)
                                                                                                                                gCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGAcaacaa        <  1:504380/149‑1 (MQ=255)
                                                                                                                                      ggCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGTAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAgg  >  1:359694/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                               
TTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGG  >  NC_000913_3_pae_tpiA/2679385‑2679667

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: