Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA3,657,5990TA25.0% 41.4 / 14.8 24intergenic (+32/+387)hdeD/arrSacid‑resistance membrane protein/antisense sRNA ArrS, function unknown
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base A (3/3);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAT  >  NC_000913_3_hsa_tpiA/3657455‑3657747
                                                                                                                                                |                                                                                                                                                    
gTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAt                                                                                                                                                  <  1:394538/149‑1 (MQ=255)
                 tCTGTAACGCTGCTTTCCACACTGGTGGGAATTGAACTGGTATTGAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                                                                 <  2:406126/149‑1 (MQ=255)
                       aCGCTGGTTTCCACACTGGTGGGAATTGAACTGATATGTCGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAct                                                                                                                           <  1:249055/149‑1 (MQ=255)
                        cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                                                                        <  1:193511/135‑1 (MQ=255)
                        cGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                                                                        >  2:193511/1‑135 (MQ=255)
                           tGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGCTCATAGTCACTTGATATAACTATGCTCCGATTAATActatct                                                                                                                       >  2:337764/1‑149 (MQ=255)
                                              aaTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTTATGTGACTATGATCCGATTCATACTCTCACCCCTACCAACTGTTGTaa                                                                                                    >  2:6246/1‑148 (MQ=255)
                                                 cgAACTTATATTGAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                                                                 <  1:417883/148‑1 (MQ=255)
                                                       gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                                                                       >  1:422683/1‑105 (MQ=255)
                                                       gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGCATCATAGTCACTTGATGTGACTATGAt                                                                                                                                       <  2:422683/105‑1 (MQ=255)
                                                                                       ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt                                                                       >  1:54893/1‑137 (MQ=255)
                                                                                       ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt                                                                       <  2:54893/137‑1 (MQ=255)
                                                                                                 cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCAGAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAc                                                 <  2:263498/149‑1 (MQ=255)
                                                                                                                atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATc                                  <  2:331682/149‑1 (MQ=255)
                                                                                                                 ttACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCt                                 >  1:56749/1‑149 (MQ=255)
                                                                                                                     cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTTAACCTACTGAATTTCATCTAAtt                             >  1:112875/1‑149 (MQ=255)
                                                                                                                       cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAtttt                           >  1:260939/1‑149 (MQ=255)
                                                                                                                                 ccGGACCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGCATTTCATCTAATTTTTTTCACTCTa                 >  2:300019/1‑149 (MQ=255)
                                                                                                                                     aTCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                     >  2:48950/1‑61 (MQ=21)
                                                                                                                                     aTCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                                     <  1:48950/61‑1 (MQ=21)
                                                                                                                                       cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTaa                               >  1:188953/1‑129 (MQ=255)
                                                                                                                                       cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTaa                               <  2:188953/129‑1 (MQ=255)
                                                                                                                                        atcatCACTTGATGTGACGTTTATCTGATTAGCACCCTCTCCGTTACGAAGCGTTGTAGATCAATTTCCCACTATCATTGAAATAATTACCTTCTAGTGATTATTTCACTCTACTGAATTTCATCTAATTTTTTTCTCTCTATGGCaaa          <  2:34860/145‑1 (MQ=255)
                                                                                                                                         tatgcACTTGGTGTGACTATGGTCCGATTAATACTTTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAAt         <  2:225778/145‑1 (MQ=255)
                                                                                                                                                ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTAACAATTTCAACCCACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAt  <  2:265123/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTGCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGGCAAATTAGCCAT  >  NC_000913_3_hsa_tpiA/3657455‑3657747

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: