Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 4,111,662 G→A L179L (CTG→TTG)  yiiQ ← DUF1454 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA4,111,6620GA97.5% 126.8 / NA 40L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/16);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTG  >  NC_000913_3_hsa_tpiA/4111543‑4111801
                                                                                                                       |                                                                                                                                           
agcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATa                                                                                                               <  2:624363/150‑1 (MQ=255)
agcagGGCTGTGAAACAGTATCGAGAATCAGCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATa                                                                                                               <  2:110158/150‑1 (MQ=255)
     gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATaccccc                                                                                                          >  2:266402/1‑145 (MQ=255)
        tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                                                                                                         >  2:648696/1‑148 (MQ=255)
        tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                                                                                                         <  1:648696/148‑1 (MQ=255)
          tgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGcc                                                                                                     >  2:200040/1‑149 (MQ=255)
            aaaCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAcc                                                                                                   >  1:688051/1‑150 (MQ=255)
              aCAGTAGAGAGATTCAGCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCAGAGACAACATAACGCAGTGCACCtt                                                                                                 <  2:92304/150‑1 (MQ=255)
                     agagATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATaa                                                                                                              <  1:3824/130‑1 (MQ=255)
                     agagATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATaa                                                                                                              >  2:3824/1‑130 (MQ=255)
                        gaTTCATCGGCACAAAGCATATGCTGATTGCCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCATCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTCCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                       <  2:688051/150‑1 (MQ=255)
                                     aacggttttgcttgTTGTCATTTCTTTAAATCTTCAAGCGATTCAGATAGCGCCAGCTTGATCGGTTCAACGGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACGTTCTGTCTCGGTGTAATAACgtttt                                                                          >  2:291487/7‑149 (MQ=255)
                                     aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCCCCTTCGGTCTCGGTGTAATAACgtttt                                                                          >  2:671121/1‑149 (MQ=255)
                                      aaGGCATTGCTTTTTGTCATTTATGCAAACCTTCCAGTGATTTAGATAGCGCCTGCTTACTTGGTTCAACAGCGAAGGGCAACCCCTTTTCGCCGTTGTCCGCGACAAAATAACGCAGTGCGCCTTCTGTCTCGGTGTAATAACgtttgt                                                                         <  2:252423/150‑1 (MQ=255)
                                          ctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGttttt                                                                     >  1:321518/1‑150 (MQ=255)
                                            ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCGAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGGCTCGGTGTAATAACGTTTGTTTTTcc                                                                   >  1:619917/1‑150 (MQ=255)
                                                        aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTtagta                                                       >  1:206132/1‑150 (MQ=255)
                                                            aTTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGAAAACATAACGCCGTGCACCTTCTGTATCGGTGTAATAACGTTTTTTTTTACCCGCGGTtagttgtga                                                   >  2:544183/1‑145 (MQ=255)
                                                               cAAACCCTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTTTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                <  1:230642/150‑1 (MQ=255)
                                                               cAAACATTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCCTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                <  1:582828/150‑1 (MQ=255)
                                                                aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGc                                               >  2:9907/1‑150 (MQ=255)
                                                                     ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTACTAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttt                                          >  1:589507/1‑150 (MQ=255)
                                                                         aGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                      >  2:180910/1‑150 (MQ=255)
                                                                           cGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATACCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGGTTTTTCTg                                    >  2:89106/1‑150 (MQ=255)
                                                                            gATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTgg                                   >  1:571459/1‑150 (MQ=255)
                                                                              ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCg                                                                                                                                <  2:469451/55‑1 (MQ=255)
                                                                              ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCg                                                                                                                                >  1:469451/1‑55 (MQ=255)
                                                                                   aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACCTTTGTTTTTCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                            <  1:444051/150‑1 (MQ=255)
                                                                                      gcgcCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                         <  1:375922/150‑1 (MQ=255)
                                                                                       cgcCAGCTGAATCGGTTCAACAGCGAAGGTCAACCCCTGTTCGCCGGTGTCCGCGACAACATAACGCACTGCACCTTCTGTCTCGGTGTAATAACATGTGTTTATCCCCGCCGTTACTAGCGACTGCGGGTTTTTCTGGCGTTGTGtttt                        >  1:596202/1‑150 (MQ=255)
                                                                                              ttAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATACAGTTTGTTTTTCCCCACCATTAGTAGCCACTGCATTTTTTTTTGGCTTTTTGTGTTGGCCAta                 >  2:303758/1‑149 (MQ=255)
                                                                                                 aTCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt              >  1:609225/1‑150 (MQ=255)
                                                                                                    ggTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg           <  1:180910/150‑1 (MQ=255)
                                                                                                    ggTTAAACAGCGTAGTGCAGCCTCTTTTCGCCGTCGTCGCCGACTACATAACGTTGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg           <  2:612285/150‑1 (MQ=255)
                                                                                                        cAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa       <  1:325412/150‑1 (MQ=255)
                                                                                                        agacaaagAAGGTCAACCCCCTTCCGCAGATGTCCGTGACAACATAACGTTGTGTACCTCCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTCTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa       <  2:20608/143‑1 (MQ=255)
                                                                                                          aCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt     >  1:75349/1‑150 (MQ=255)
                                                                                                             gCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGtgtg  >  2:241393/1‑150 (MQ=255)
                                                                                                                    tCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                       <  2:238165/122‑1 (MQ=255)
                                                                                                                    tCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                       >  1:238165/1‑122 (MQ=255)
                                                                                                                       |                                                                                                                                           
AGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTG  >  NC_000913_3_hsa_tpiA/4111543‑4111801

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