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.

ATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAT  >  NC_000913_3_hsa_tpiA/1211932‑1212227
                                                                                                                                                      |                                                                                                                                                        
atTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑a                                                                                                                                                           >  1:841694/1‑149 (MQ=255)
     cccATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑c‑ccc‑c                                                                                                                                                  >  2:717167/1‑149 (MQ=255)
       cATATGGCTGAAAAAGATATACCACATGTAGGTTGATTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGATTATTTT‑A‑‑c‑ctc‑cc‑c                                                                                                                                               >  2:433395/1‑149 (MQ=255)
        atatGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTAACCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑c‑ccc‑cc‑c‑c                                                                                                                                             >  2:1042024/1‑149 (MQ=255)
                   aaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAaaata                                                                                                                                  >  2:693183/1‑149 (MQ=255)
                   aaaGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAaaata                                                                                                                                  >  2:319444/1‑149 (MQ=255)
                            ccACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATGTGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGGGTATATT‑T‑‑C‑CCCACC‑C‑CTCAAAGAAATAaatttatt                                                                                                                          >  1:323325/1‑149 (MQ=255)
                              acaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTaaa                                                                                                                       >  1:488865/1‑149 (MQ=255)
                               caTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTaaa                                                                                                                       >  1:226331/1‑149 (MQ=255)
                                      gTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCACATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGc                                                                                                               >  2:658265/1‑149 (MQ=255)
                                       ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCt                                                                                                              <  1:974717/149‑1 (MQ=255)
                                       ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCt                                                                                                              >  1:772179/1‑149 (MQ=255)
                                         gAATTACCGTGTCCATTACTATCCATTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTGTATT‑A‑‑C‑CCC‑CCTT‑CTCAAAAAAATAAATTAATTAAAATGATGGCtt                                                                                                             <  2:752024/149‑1 (MQ=255)
                                           aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCaaaaaaa                                                                                                                                    >  2:201009/1‑124 (MQ=255)
                                           aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCaaaaaaa                                                                                                                                    <  1:201009/124‑1 (MQ=255)
                                           aTTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGATCATTGGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CGCAAAAAAATAAATTAATTAAAATTATGGTTtat                                                                                                           >  1:91329/1‑149 (MQ=255)
                                             tACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTtatata                                                                                                        <  1:319444/149‑1 (MQ=255)
                                                 gtgtAAATTACTATCCACTTCATTTTTTATCTCTAATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATAAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaa                                                                                                     <  2:919412/149‑1 (MQ=255)
                                                      aaTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                                                                                <  2:226331/149‑1 (MQ=255)
                                                      aaTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                                                                                <  1:873164/149‑1 (MQ=255)
                                                           cTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATAAAAATGATGGCTTTTATAAAATAAAATTTaaa                                                                                          >  2:290371/1‑149 (MQ=255)
                                                             aTCCACTTCATTTGTTATGTCTTATCCCACTGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTTAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAg                                                                                         >  2:9587/1‑149 (MQ=255)
                                                              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTCC‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataaa                                                                                                 <  1:573117/140‑1 (MQ=255)
                                                              tCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaataaa                                                                                                 >  2:573117/1‑140 (MQ=255)
                                                                 aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAaattaa                                                                                                                            >  1:807067/1‑110 (MQ=255)
                                                                 aCTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAaattaa                                                                                                                            <  2:807067/110‑1 (MQ=255)
                                                                   ttCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATAATGGCTTCTTTAAAATAAAATATAAAGCaataa                                                                                   >  2:1020365/1‑146 (MQ=255)
                                                                     cATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCCTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTAAAAGCAAGGAAt                                                                                 >  1:415983/1‑149 (MQ=255)
                                                                       tttGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtct                                                                               >  1:950829/1‑149 (MQ=255)
                                                                           ttATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAAt                                                                           <  1:279653/149‑1 (MQ=255)
                                                                                tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑ACCC‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAttt                                                                                             >  1:548830/1‑127 (MQ=255)
                                                                                tCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑ACCC‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAttt                                                                                             <  2:548830/127‑1 (MQ=255)
                                                                                  ttATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGCGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAaattaa                                                                                                                            >  2:324299/1‑93 (MQ=255)
                                                                                  ttATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAaattaa                                                                                                                            <  1:324299/93‑1 (MQ=255)
                                                                                     tCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTaaa                                                                 >  2:504288/1‑149 (MQ=255)
                                                                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaaa                                                                                          <  2:560341/120‑1 (MQ=255)
                                                                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCCAAC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTaaa                                                                                          >  1:560341/1‑120 (MQ=255)
                                                                                         aCGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACaaaa                                                            <  1:658265/149‑1 (MQ=255)
                                                                                             tATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAAAAAGGAATCTCAATGGATATTAAACAAAATga                                                         >  2:636521/1‑149 (MQ=255)
                                                                                               tgtattatGGTTTATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgaga                                                       <  2:35151/144‑1 (MQ=255)
                                                                                                         ttCCTTAGGACCTGTATTACTTGATGTTGCATATCAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAAGAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGaaa                                             <  2:323325/149‑1 (MQ=255)
                                                                                                            aTTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATCAACTTTAAAACAAGGAATCTCAAGGTGGGTTAAACAAAATGAGATTTTGTGaagaca                                          >  1:243280/1‑145 (MQ=255)
                                                                                                                gtaTGTTTATTTCTAGATTTTGCATATGAATATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATaaa                                     <  2:772179/147‑1 (MQ=255)
                                                                                                                  aTGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTACATGAATGAGAAACGAAATGAGATTTTGTGAAAGCAATAAatt                                   >  1:978047/1‑149 (MQ=255)
                                                                                                                    gTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAattct                                 >  2:738535/1‑147 (MQ=255)
                                                                                                                            cTTGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTc                          >  2:448129/1‑149 (MQ=255)
                                                                                                                             ttGATTTTGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCg                         >  1:632822/1‑149 (MQ=255)
                                                                                                                                  tttGCATATGAGTATATT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGGTGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTa                    <  1:1042024/149‑1 (MQ=255)
                                                                                                                                  tttGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑CACTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTa                    <  2:755935/149‑1 (MQ=255)
                                                                                                                                       atatGAGTATATT‑C‑‑CACAC‑CC‑C‑CTCAAACAAAAAACTTCCTAAAAATAATCACCCATAGATAGCAAAAATGGAAGTAAAGAAGGTCAATGGATGATAAACAAAAAGAGATTTTGTGAAAGCAATAAAATATTGACTTCGTTTTAGAttt               >  2:735199/1‑149 (MQ=255)
                                                                                                                                            aGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAg         >  1:178208/1‑149 (MQ=255)
                                                                                                                                              tataTT‑AC‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGATTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAAAGAGATTTTGTTAAAGCAATAAATTATTAACTTCGTTTTATATTTGTTTAgt        >  1:502977/1‑148 (MQ=255)
                                                                                                                                                   t‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAt  >  1:576821/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                        
ATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATT‑A‑‑C‑CCC‑CC‑C‑CTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAT  >  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: ‑