Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE707,9100AG74.2% 98.5 / 40.5 62intergenic (‑1/‑121)znuC_1/ideRHigh‑affinity zinc uptake system ATP‑binding protein ZnuC/Iron‑dependent repressor IdeR
Reads supporting (aligned to +/- strand):  ref base A (9/7);  new base G (21/25);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTGAAGAA  >  USA300TCH1516_ALE/707773‑708045
                                                                                                                                         |                                                                                                                                       
ttGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGt                                                                                                                                      <  2:132098/141‑1 (MQ=255)
 tGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGtt                                                                                                                                     >  2:312713/1‑141 (MQ=255)
    tAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTtat                                                                                                                                  >  2:155368/1‑141 (MQ=255)
    tAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTtat                                                                                                                                  >  1:390100/1‑141 (MQ=255)
          aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATAcctcc                                                                                                                            <  1:731775/141‑1 (MQ=255)
          aaGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcctcc                                                                                                                            <  2:341785/141‑1 (MQ=255)
           aGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATAcctcct                                                                                                                           >  1:679755/1‑141 (MQ=255)
           aGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATAcctcct                                                                                                                           >  1:622114/1‑141 (MQ=255)
           aGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTGTCTAACAAGCGTTTATAcctcct                                                                                                                           >  2:24179/1‑141 (MQ=255)
               ttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAAtt                                                                                                                       <  1:176372/141‑1 (MQ=255)
               ttACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAAtt                                                                                                                       <  1:229812/141‑1 (MQ=255)
                     gCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGt                                                                                                                 <  2:253284/141‑1 (MQ=255)
                          aTTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAgg                                                                                                            >  2:342064/1‑141 (MQ=255)
                              gggCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACGTGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTaa                                                                                                        >  1:554427/1‑141 (MQ=255)
                                    aTGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCtaat                                                                                                  >  1:185317/1‑141 (MQ=255)
                                         aCCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATATTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTCACCTaattaatt                                                                                             >  2:187557/1‑141 (MQ=255)
                                               tattTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTataa                                                                                       <  1:713141/141‑1 (MQ=255)
                                                  ttCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTAtaataa                                                                                    >  2:172450/1‑141 (MQ=255)
                                                  ttCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTAtaataa                                                                                    >  1:722016/1‑141 (MQ=255)
                                                    cgccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAAc                                                                                  <  1:464742/141‑1 (MQ=255)
                                                     gccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACATACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACt                                                                                 >  1:471038/1‑141 (MQ=255)
                                                      ccgcGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTg                                                                                <  1:549836/141‑1 (MQ=255)
                                                        ccgTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAANAGTTTAGGCTAACCTAATTAATTGTATAATAAACTggg                                                                              <  1:161443/140‑3 (MQ=255)
                                                             cTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAatat                                                                         <  1:88934/141‑1 (MQ=255)
                                                              tGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAACAACAGATTTACATATTTTGTTTCTAACAAGAGCACACACCTCCTCATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATAtt                                                                        <  1:492558/141‑1 (MQ=255)
                                                               ggTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATAttt                                                                       >  1:100354/1‑141 (MQ=255)
                                                                    cGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCAt                                                                  >  2:140417/1‑141 (MQ=255)
                                                                    cGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCAt                                                                  >  2:140418/1‑141 (MQ=255)
                                                                     gATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATg                                                                 <  1:342539/141‑1 (MQ=255)
                                                                     gATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATg                                                                 <  1:563962/141‑1 (MQ=255)
                                                                         aGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTAGTTTCTAACAAGCGTTTATACCGCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCaa                                                             >  1:193288/1‑141 (MQ=255)
                                                                          ggAAATGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                                                            <  1:90405/141‑1 (MQ=255)
                                                                          ggAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAg                                                            <  1:456284/141‑1 (MQ=255)
                                                                            aaaTGTTTTTAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTa                                                          <  1:782276/141‑1 (MQ=255)
                                                                                   tttAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                                                   >  2:324858/1‑141 (MQ=255)
                                                                                   tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                                                   >  2:244345/1‑141 (MQ=255)
                                                                                   tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                                                   >  1:318094/1‑141 (MQ=255)
                                                                                   tttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGt                                                   >  1:23538/1‑141 (MQ=255)
                                                                                    ttAAGTACATGCTTATTACCTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                                  <  1:236369/141‑1 (MQ=255)
                                                                                    ttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                                  >  1:147266/1‑141 (MQ=255)
                                                                                    ttAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTATAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATAGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTg                                                  >  1:432073/1‑141 (MQ=255)
                                                                                          aCATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGatataa                                            >  1:679083/1‑141 (MQ=255)
                                                                                             tGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTTCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCACGTCAAGTAAATTCGTGatataatat                                         <  1:110527/141‑1 (MQ=255)
                                                                                               cTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAg                                       >  1:597677/1‑141 (MQ=255)
                                                                                                  attaCCTAAAAACAGATTTAAATCGTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACa                                    <  1:88666/141‑1 (MQ=255)
                                                                                                   ttaCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACaa                                   <  2:294902/141‑1 (MQ=255)
                                                                                                    taCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAat                                  <  1:703930/141‑1 (MQ=255)
                                                                                                     aCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatg                                 >  2:225696/1‑141 (MQ=255)
                                                                                                     aCCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATATATCATGTCAAGTAAATTCGTGATATAATATAGACAatg                                 >  1:257240/1‑141 (MQ=255)
                                                                                                      ccTATAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgt                                >  2:319593/1‑141 (MQ=255)
                                                                                                      ccTAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgt                                <  1:431421/141‑1 (MQ=255)
                                                                                                        tAAAAACAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAatgtat                              <  2:257982/141‑1 (MQ=255)
                                                                                                           aaaCAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTAtgtg                           <  1:596553/141‑1 (MQ=255)
                                                                                                            aaCAGATTTAAATCTTTTTTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGa                          <  1:977851/141‑1 (MQ=255)
                                                                                                              cAGATTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAgg                        <  1:164088/141‑1 (MQ=255)
                                                                                                                 aTTTAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATGTAGACAATGTATGTGAGGTGa                     <  2:271294/141‑1 (MQ=255)
                                                                                                                   ttAAATCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGaaa                   <  2:149612/141‑1 (MQ=255)
                                                                                                                        tCTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATg              <  1:329878/141‑1 (MQ=255)
                                                                                                                         cTTTTGTTTCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATTTTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGt             <  1:874087/141‑1 (MQ=255)
                                                                                                                               tttCTAACAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTg       >  1:700860/1‑141 (MQ=255)
                                                                                                                                ttCTAACTAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTga      <  2:381108/141‑1 (MQ=255)
                                                                                                                                    aaCAAGCGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTgaagaa  >  1:558124/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TTGATAAGGGAAGATTTACCAGCACCATTCGGGCCCATGATACCAATTATTTCGCCGCGTACTGGTATCGATAAGGAAATGTTTTTAAGTACATGCTTATTACCTAAAAACAGATTTAAATCTTTTGTTTCTAACAAACGTTTATACCTCCTAATTAAAAGTTTAGGCTAACCTAATTAATTGTATAATAAACTGAGAATATTTATCATGTCAAGTAAATTCGTGATATAATATAGACAATGTATGTGAGGTGAAAGTATGTTAACTGAAGAA  >  USA300TCH1516_ALE/707773‑708045

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑