Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA1,212,0801.C59.2% 83.6 / 78.4 49intergenic (‑77/+623)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base . (13/7);  new base C (17/12);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAT  >  NC_000913_3_hsa_tpiA/1211932‑1212227
                                                                                                                                                      |                                                                                                                                                        
atTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTa                                                                                                                                                           >  1:841694/1‑149 (MQ=255)
     cccATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccc                                                                                                                                                  >  2:717167/1‑149 (MQ=255)
       cATATGGCTGAAAAAGATATACCACATGTAGGTTGATTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGATTATTTTA‑‑cctcccc                                                                                                                                               >  2:433395/1‑149 (MQ=255)
        atatGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTAACCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑cccccccc                                                                                                                                             >  2:1042024/1‑149 (MQ=255)
                   aaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAaaata                                                                                                                                  >  2:693183/1‑149 (MQ=255)
                   aaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAaaata                                                                                                                                  >  2:319444/1‑149 (MQ=255)
                            ccACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATGTGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGGGTATATTT‑‑CCCCACCCCTCAAAGAAATAaatttatt                                                                                                                          >  1:323325/1‑149 (MQ=255)
                              acaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTaaa                                                                                                                       >  1:488865/1‑149 (MQ=255)
                               caTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTaaa                                                                                                                       >  1:226331/1‑149 (MQ=255)
                                      gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCACATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGc                                                                                                               >  2:658265/1‑149 (MQ=255)
                                       ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCt                                                                                                              <  1:974717/149‑1 (MQ=255)
                                       ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCt                                                                                                              >  1:772179/1‑149 (MQ=255)
                                         gAATTACCGTGTCCATTACTATCCATTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTGTATTA‑‑CCCCCCTTCTCAAAAAAATAAATTAATTAAAATGATGGCtt                                                                                                             <  2:752024/149‑1 (MQ=255)
                                           aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaaaaa                                                                                                                                    >  2:201009/1‑124 (MQ=255)
                                           aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaaaaa                                                                                                                                    <  1:201009/124‑1 (MQ=255)
                                           aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGATCATTGGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCGCAAAAAAATAAATTAATTAAAATTATGGTTtat                                                                                                           >  1:91329/1‑149 (MQ=255)
                                             tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                                                                        <  1:319444/149‑1 (MQ=255)
                                                 gtgtAAATTACTATCCACTTCATTTTTTATCTCTAATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATAAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaa                                                                                                     <  2:919412/149‑1 (MQ=255)
                                                      aaTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                                                                                <  2:226331/149‑1 (MQ=255)
                                                      aaTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                                                                                <  1:873164/149‑1 (MQ=255)
                                                           cTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATAAAAATGATGGCTTTTATAAAATAAAATTTaaa                                                                                          >  2:290371/1‑149 (MQ=255)
                                                             aTCCACTTCATTTGTTATGTCTTATCCCACTGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTTAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                                                         >  2:9587/1‑149 (MQ=255)
                                                              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTCC‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataaa                                                                                                 <  1:573117/140‑1 (MQ=255)
                                                              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataaa                                                                                                 >  2:573117/1‑140 (MQ=255)
                                                                 aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaattaa                                                                                                                            >  1:807067/1‑110 (MQ=255)
                                                                 aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaattaa                                                                                                                            <  2:807067/110‑1 (MQ=255)
                                                                   ttCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATAATGGCTTCTTTAAAATAAAATATAAAGCaataa                                                                                   >  2:1020365/1‑146 (MQ=255)
                                                                     cATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCCTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTAAAAGCAAGGAAt                                                                                 >  1:415983/1‑149 (MQ=255)
                                                                       tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtct                                                                               >  1:950829/1‑149 (MQ=255)
                                                                           ttATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAAt                                                                           <  1:279653/149‑1 (MQ=255)
                                                                                tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAttt                                                                                             >  1:548830/1‑127 (MQ=255)
                                                                                tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAttt                                                                                             <  2:548830/127‑1 (MQ=255)
                                                                                  ttATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGCGTATATTACCCCCCCCCTCAAAAAAATAaattaa                                                                                                                            >  2:324299/1‑93 (MQ=255)
                                                                                  ttATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaattaa                                                                                                                            <  1:324299/93‑1 (MQ=255)
                                                                                     tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                                                 >  2:504288/1‑149 (MQ=255)
                                                                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaaa                                                                                          <  2:560341/120‑1 (MQ=255)
                                                                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCAACCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaaa                                                                                          >  1:560341/1‑120 (MQ=255)
                                                                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaaa                                                            <  1:658265/149‑1 (MQ=255)
                                                                                             tATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAAAAAGGAATCTCAATGGATATTAAACAAAATga                                                         >  2:636521/1‑149 (MQ=255)
                                                                                               tgtattatGGTTTATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgaga                                                       <  2:35151/144‑1 (MQ=255)
                                                                                                         ttCCTTAGGACCTGTATTACTTGATGTTGCATATCAGTATATTACCCCCCCCCTCAAAAAAAGAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGaaa                                             <  2:323325/149‑1 (MQ=255)
                                                                                                            aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATCAACTTTAAAACAAGGAATCTCAAGGTGGGTTAAACAAAATGAGATTTTGTGaagaca                                          >  1:243280/1‑145 (MQ=255)
                                                                                                                gtaTGTTTATTTCTAGATTTTGCATATGAATATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATaaa                                     <  2:772179/147‑1 (MQ=255)
                                                                                                                  aTGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTACATGAATGAGAAACGAAATGAGATTTTGTGAAAGCAATAAatt                                   >  1:978047/1‑149 (MQ=255)
                                                                                                                    gTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAattct                                 >  2:738535/1‑147 (MQ=255)
                                                                                                                            cTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTc                          >  2:448129/1‑149 (MQ=255)
                                                                                                                             ttGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCg                         >  1:632822/1‑149 (MQ=255)
                                                                                                                                  tttGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGGTGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTa                    <  1:1042024/149‑1 (MQ=255)
                                                                                                                                  tttGCATATGAGTATATTA‑‑CCCCCCCACTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTa                    <  2:755935/149‑1 (MQ=255)
                                                                                                                                       atatGAGTATATTC‑‑CACACCCCCTCAAACAAAAAACTTCCTAAAAATAATCACCCATAGATAGCAAAAATGGAAGTAAAGAAGGTCAATGGATGATAAACAAAAAGAGATTTTGTGAAAGCAATAAAATATTGACTTCGTTTTAGAttt               >  2:735199/1‑149 (MQ=255)
                                                                                                                                            aGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAg         >  1:178208/1‑149 (MQ=255)
                                                                                                                                              tataTTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGATTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAAAGAGATTTTGTTAAAGCAATAAATTATTAACTTCGTTTTATATTTGTTTAgt        >  1:502977/1‑148 (MQ=255)
                                                                                                                                                   tA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAt  >  1:576821/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
ATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAT  >  NC_000913_3_hsa_tpiA/1211932‑1212227

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: