Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,033,739 (T)8→7 noncoding (91/91 nt) rseX → sRNA antisense regulator of ompA and ompC translation, Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,033,7320T.97.3% 163.1 / ‑3.2 37noncoding (84/91 nt)rseXsRNA antisense regulator of ompA and ompC translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (18/18);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATCTTTGATATCGGCTTTATTAATTCATTTAATCAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATACATTTATAATTTTCATGAATATTTATATTTAGAATTCATA  >  NC_000913_3_hsa_tpiA/2033586‑2033876
                                                                                                                                                  |                                                                                                                                                
atCTTTGATATCGGCTTTATTAATTCATTTAATCAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGAATGATGCCGGCttt                                                                                                                                                >  1:370120/1‑149 (MQ=255)
  cTTTGATATCGGCTTTATTAATTCATTTAATCAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCttttt                                                                                                                                              <  1:225388/149‑1 (MQ=255)
         atCGGCTTTATTAATTCATTTAATCAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGcc                                                                                                                                      >  1:126970/1‑149 (MQ=255)
         atCGGCTTTATTAATTCATTTAATCAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGcc                                                                                                                                      <  1:555506/149‑1 (MQ=255)
           cgccTTTATTAATTTATTTAATAAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTTTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCtt                                                                                                                                    <  2:370120/146‑1 (MQ=255)
                    taATTCATTTAATCAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAAtgt                                                                                                                           <  1:322707/149‑1 (MQ=255)
                                    tatatTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCa                                                                                                           <  1:260480/149‑1 (MQ=255)
                                             aCTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTgg                                                                                                  >  1:481393/1‑149 (MQ=255)
                                              cTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGgc                                                                                                 <  1:463642/149‑1 (MQ=255)
                                                 aTTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTTTTTTTGTTTGCTGCCGGCTTTTTCTTTGCCTTCTTTAATGTGCGCCTGATGACCCCACCCGTTTGgcgcc                                                                                              >  2:292774/1‑148 (MQ=255)
                                                          ataCCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTCTGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCAtt                                                                                     <  1:269788/149‑1 (MQ=255)
                                                              cTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGAt                                                                                                                  >  2:267015/1‑116 (MQ=255)
                                                              cTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGAt                                                                                                                  <  1:267015/116‑1 (MQ=255)
                                                               tttttATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATTTTTTTTTGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTTATCCCACCAGCCGTTTGTCGCAACAATCTTTGTTAc                                                                                >  2:100172/1‑149 (MQ=255)
                                                                    attattTTGTTGTATGATGTTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCGTTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCt                                                                           <  2:229700/149‑1 (MQ=255)
                                                                    attattCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCt                                                                           >  1:356233/1‑149 (MQ=255)
                                                                       attCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATg                                                                        >  2:157417/1‑149 (MQ=255)
                                                                        ttCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAAt                                                                                                                             >  2:136745/1‑95 (MQ=255)
                                                                        ttCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAAt                                                                                                                             <  1:136745/95‑1 (MQ=255)
                                                                                 atgatgCTTCCGTTATTAGCCTTTTATCGGCTTGTTTATATTGTTTGGGCCGGCATGGTGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATGATTGATACCCCCGATCTTTCCGCTcc                                                              >  1:126572/1‑148 (MQ=255)
                                                                                   gatgCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAAc                                                            >  1:28850/1‑149 (MQ=255)
                                                                                       cTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCAt                                                        <  2:356233/149‑1 (MQ=255)
                                                                                         tCCGTTATTAGCTTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATaa                                                      <  2:126572/149‑1 (MQ=255)
                                                                                         tCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATaa                                                      >  2:390984/1‑149 (MQ=255)
                                                                                             ttattaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACata                                                  >  1:434176/1‑149 (MQ=255)
                                                                                              tattaGCTTTTTACCGTCTGGTTTACATTTGTTAGGCCGGCATGATTCCGGCTTTTTTTATCCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTCGCGCAACAATCCTTGATACCCCCTAAGTTTCCGCTCAACTCATAACatat                                                 <  2:116382/149‑1 (MQ=255)
                                                                                              tattaGCCTTTTATCGTTTGGTTTATATTTTTGGGGCGGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACatat                                                 <  1:203360/149‑1 (MQ=255)
                                                                                               attaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCa                                                                                                           <  2:226923/90‑1 (MQ=255)
                                                                                               attaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCa                                                                                                           >  1:226923/1‑90 (MQ=255)
                                                                                                        ttATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATACAt                                       <  1:27606/149‑1 (MQ=255)
                                                                                                             gTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGc                                                                 >  2:187917/1‑118 (MQ=255)
                                                                                                             gTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGc                                                                 <  1:187917/118‑1 (MQ=255)
                                                                                                              tCTTGTTTATATTTTTTGTGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACCCCAGCCGTTTGGCGCAACAATCATTGATATTCCCTATGTTTCCGCTCAACTAATAACATATTGATATACATTTATaa                                 >  2:484795/1‑149 (MQ=255)
                                                                                                                   tttATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATACATTTATAATTTTc                            <  2:550358/149‑1 (MQ=255)
                                                                                                                    atataTTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTAAACTCATAACATATTGCTATACATTTATAATTTTCa                           >  1:567262/2‑149 (MQ=255)
                                                                                                                       tatTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATACATTTATAATTTTCATGa                        >  1:356929/1‑149 (MQ=255)
                                                                                                                                ggCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATACATTTATAATTTTCATGAatatttata               >  1:320033/1‑149 (MQ=255)
                                                                                                                                    ggCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCTCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATACATTTATAATTTTCATGAAtatttatattta           <  2:126970/149‑1 (MQ=255)
                                                                                                                                             ccGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATTCATTTATAATTTTCATGAATATTTATATTTAGAATTCATa  >  2:123391/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCTTTGATATCGGCTTTATTAATTCATTTAATCAATATATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCATAACATATTGATATACATTTATAATTTTCATGAATATTTATATTTAGAATTCATA  >  NC_000913_3_hsa_tpiA/2033586‑2033876

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: