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,6620GA95.7% 141.2 / NA 47L179L (CTG→TTG) yiiQDUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/16);  minor base T (0/2);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGA  >  NC_000913_3_hsa_tpiA/4111519‑4111805
                                                                                                                                               |                                                                                                                                               
aGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCt                                                                                                                                            <  1:216115/149‑1 (MQ=255)
  gAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCttt                                                                                                                                          >  2:433064/1‑149 (MQ=255)
    agagGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCACACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTc                                                                                                                                        <  2:381845/149‑1 (MQ=255)
              cccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAAGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCCCCTTTGTcc                                                                                                                              >  2:232714/1‑149 (MQ=255)
               ccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTAATTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACACCGAAGGTCAACCCCTTTTCTCCGTTATCcg                                                                                                                             >  2:516670/1‑149 (MQ=255)
                     ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaa                                                                                                                       >  1:248148/1‑149 (MQ=255)
                               cTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                                                                                                             <  1:28849/149‑1 (MQ=255)
                                tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGt                                                                                                            >  1:38354/1‑149 (MQ=255)
                                   gAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                                                                                                         >  2:537147/1‑149 (MQ=255)
                                        aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAGTCGGGTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc                                                                                                    <  1:415299/149‑1 (MQ=255)
                                         gTATGGAGAGTCAGCGGCACAAAGGCCTTGCTATCTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCACCTTAATCGGTTCAACAGCTAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                   <  2:85706/149‑1 (MQ=255)
                                         gTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                   >  1:498361/1‑149 (MQ=255)
                                         gTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt                                                                                                   <  1:147317/149‑1 (MQ=255)
                                             agagATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtct                                                                                               >  1:12530/1‑149 (MQ=255)
                                                     aTCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtac                                                                                       >  2:167408/1‑148 (MQ=255)
                                                      tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCACGCAATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGATCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGtaat                                                                                      >  2:100378/1‑149 (MQ=255)
                                                        ggCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAGCTTCTGTCTCGGTGtaataa                                                                                    >  1:437169/1‑149 (MQ=255)
                                                               aGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACttttgt                                                                             >  2:296480/1‑149 (MQ=255)
                                                               aGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttgt                                                                             <  1:239206/149‑1 (MQ=255)
                                                                                   tATTCAACCCGTCAAGCGATTCAGATAGCGCCAGCGTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCTTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGc                                                         <  1:275548/149‑1 (MQ=255)
                                                                                        aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGCAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                    >  2:250830/1‑149 (MQ=255)
                                                                                             ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATCACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAggt                                               >  1:479228/1‑147 (MQ=255)
                                                                                                 ctcGATTGAGATAGCGCCAGTTTAATCGGTTCAACAGTGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCCGCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttttt                                           <  2:522481/147‑1 (MQ=255)
                                                                                                    gATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTGGCGACTGCAGTTTTTTCTg                                        >  2:355364/1‑149 (MQ=255)
                                                                                                    gATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTcag                                        >  1:522481/1‑147 (MQ=255)
                                                                                                      ttCAGATAGCGCCAGCTTAATCGGTTAAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGc                                      >  1:447111/1‑149 (MQ=255)
                                                                                                       tCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACCACATAACGCAGAGCACCTTCTGTCTCTGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTACAGTTTTGTCTTGCt                                     >  2:300014/1‑149 (MQ=255)
                                                                                                       tCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCGCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                     >  1:167091/1‑149 (MQ=255)
                                                                                                         aGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCttt                                   >  2:422728/1‑149 (MQ=255)
                                                                                                           aTAACGCAAGCTTAAGCGGTGCAACAACGAAGGTCATCCCCTTTTTGACGTATTCCGAGACAAAATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                                 <  2:13131/149‑1 (MQ=255)
                                                                                                                 ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAc                                                                                                        <  1:159062/72‑1 (MQ=255)
                                                                                                                 ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAc                                                                                                        >  2:159062/1‑72 (MQ=255)
                                                                                                                  cAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                          <  1:510285/149‑1 (MQ=255)
                                                                                                                     ctccctCGGTTCACCAGCGAAGGTCATCCCCTTTTCGCCGTTGTCCGCGAACACATAACGCAGTGCACCTTCTGGCTCGGTGTAATAACTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                       <  1:300014/144‑1 (MQ=255)
                                                                                                                       tAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                     <  1:172876/149‑1 (MQ=255)
                                                                                                                          tCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAGTAACGTTTGTTTTTGCACGCCGTTAGTTGGGACTGCCGTTTTTTTTGGCTTTGTGTTTTGGTCCTTact                  >  2:468716/1‑147 (MQ=255)
                                                                                                                          tCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                             <  1:419424/122‑1 (MQ=255)
                                                                                                                          tCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt                                             >  2:419424/1‑122 (MQ=255)
                                                                                                                          gcGGTACAACAGTGAAGGTCAACCCCTTTTCGCCGTTGCCCGCTGTAACTTAATCAAGTGCACCTTCTGTCTCGGTGGAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                  <  2:166423/148‑1 (MQ=255)
                                                                                                                             gTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTGATggg               >  1:373863/1‑149 (MQ=255)
                                                                                                                             gTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGAACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCGGTTTTTTCTGGCTTTGTGTTTTGGTCATTTATggg               >  1:203492/1‑149 (MQ=255)
                                                                                                                                  aCAGAGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATACCGAAGTGCACCTTCTGTCTCGGTGTCATAAAGGTCGTTTATACCCGCCTCTAGTACCGACTGCAGTTTTTTCTTCCTTAATGGCTTGGTCCTCAATAGGGCGat          >  2:158224/1‑148 (MQ=255)
                                                                                                                                   cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt         <  1:471132/149‑1 (MQ=255)
                                                                                                                                   cAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTGTTGGTCATTAATGGGGTGAgt         >  1:10643/1‑149 (MQ=255)
                                                                                                                                   cAGCGAAGGCCAACCCCTTTATGCCGTTGTCCGCGAAAACATAACGCCTTGCTCCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt         <  2:167091/149‑1 (MQ=255)
                                                                                                                                      cGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGtgtg      >  1:173641/1‑149 (MQ=255)
                                                                                                                                          ggTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGGGCGTGTGCGGa  >  1:195379/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
AGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGA  >  NC_000913_3_hsa_tpiA/4111519‑4111805

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