Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,611,7230GA25.4% 110.5 / 12.9 57K196K (AAG→AAAnikEnickel transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (7/34);  new base A (13/1);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 6.32e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGGAGAGAAAT  >  CP009273/3611590‑3611853
                                                                                                                                     |                                                                                                                                   
tAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCt                                                                                                                                   <  1:4022363/136‑1 (MQ=255)
tAAGCGGCGGCCAGCGCCCGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTac                                                                                                                                 <  1:2848731/138‑1 (MQ=255)
  aGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTacaac                                                                                                                              <  1:1047689/139‑1 (MQ=255)
    cggcggCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCATTCGCCTGCTGGAAAAACTacaacaa                                                                                                                            >  1:97437/1‑139 (MQ=255)
          ccAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTacaacaa                                                                                                                            <  1:5032333/133‑1 (MQ=255)
               tCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTgg                                                                                                                    <  1:2680470/136‑1 (MQ=255)
                    cgcgTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCATTCGCCTGCTGAAAAAACTACAAAAAAAGTTTGGCAcccg                                                                                                              >  1:2453420/1‑135 (MQ=255)
                       gTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCgcctgc                                                                                                          <  1:4971520/138‑1 (MQ=255)
                         cTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCCTTCGCCTTCTGTAAAAACTGCCACAACAGTTTTGCACCgcctgcccg                                                                                                       >  1:2573866/1‑137 (MQ=255)
                           gCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGtt                                                                                                     <  1:2485418/139‑1 (MQ=255)
                                   cgcgcgCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCAccc                                                                                             <  1:1307395/139‑1 (MQ=255)
                                    tcgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCa                                                                                            <  1:4054854/138‑1 (MQ=255)
                                    gcgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCa                                                                                            <  1:1023599/139‑1 (MQ=255)
                                    gcccgttGGCCGGCGGACCGAAACTACTGATTGTGGATTAAGCCCGTTCTAACCTTGGTCCCTTTTTACAGGCTGGTGTCATTCGCCTCCTGAAAAAACTACAACAACAGTTTGGCATCGCCTGCCTGTTCATCACCCa                                                                                            <  1:796816/133‑1 (MQ=255)
                                     cgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTTGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCAc                                                                                           >  1:1199646/1‑139 (MQ=255)
                                     cgcgcTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCAATTGCCTGCTGAAAAAACTACAACAAAAGGTTGGCACCGCCTGCCTGTTCATTACCCAc                                                                                           >  1:318368/1‑139 (MQ=255)
                                          tGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGGGTTACAGGCGGGGGGCATTCGCCTGCTTAAAAAACTACAAAAAAAGTTTGGGACCGGCTGCCTGTTTATCACCCACGACtt                                                                                      >  1:1200616/1‑139 (MQ=255)
                                             cGGTTGAACCCCAACTTCTGATTTTGGATTAAGCCGTTTTTAACCTTGATCTCGTTTTACAGGCGGGTGTCTTTCGCCTCCTGAAAAATCTCCAACAACCGTTTGGCACCGCCTGCCTGTTCATCACCCCCGACTTAc                                                                                    <  1:2095254/138‑1 (MQ=255)
                                                tCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTacaacaa                                                                                                                            <  1:3844811/95‑1 (MQ=255)
                                                tCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGcc                                                                                 <  1:2901047/138‑1 (MQ=255)
                                                         aaCTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACg                                                                       <  1:962708/139‑1 (MQ=255)
                                                          actactGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCAATCGCCTGCTGAAAAAAATACAAAAAAAGGTTggga                                                                                                                  >  1:1263478/1‑93 (MQ=255)
                                                            tactGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCtt                                                                   <  1:1000920/139‑1 (MQ=255)
                                                               tGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGGTATTTGCCTTCTTAAAAAACCACAACAACAATTTGTGGCCGCCCTCCCTTTCCTCCCCCCCACCTTTCGCCTGGTAGAACGcgctt                                                                 >  1:3572178/1‑134 (MQ=255)
                                                                  ttCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTgg                                                                                                                    <  1:3931867/85‑1 (MQ=255)
                                                                  ttCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCgcc                                                                                                             <  1:4984886/92‑1 (MQ=255)
                                                                   tCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAg                                                            <  1:3232720/139‑1 (MQ=255)
                                                                     tGGATGAAGCCGTTTCTACCCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAATCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCg                                                          <  1:1642672/139‑1 (MQ=255)
                                                                          gAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:4692085/139‑1 (MQ=255)
                                                                          gAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:875372/139‑1 (MQ=255)
                                                                          gAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:990461/139‑1 (MQ=255)
                                                                          gAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:1695286/139‑1 (MQ=255)
                                                                          gAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:2467091/139‑1 (MQ=255)
                                                                          gAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:2639052/139‑1 (MQ=255)
                                                                           aaGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:3562084/138‑1 (MQ=255)
                                                                           aaGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTTATTCGCCTTCTTAAAAAAATACAACAAAAGTTTGGGCCCCCCTGCCCGTTCATCCCCCCCGACCTATGCCTCGTGAAAAGCTTTTGCCAGCGGGTaaa                                                    >  1:1306080/1‑138 (MQ=255)
                                                                            aGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGTGGTCATTCCCCTGCTGAAAAAGCTACAACAACAGGTTGGCCCCGCCCGCCTGTTTATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGGAATg                                                   >  1:1513488/1‑139 (MQ=255)
                                                                            aGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGGGTCAGTTGCCTTCTGAAAAAACTACAACAAAAGTTTTGCACCGCCCGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTTCCAGCGGGTAATg                                                   >  1:980345/1‑139 (MQ=255)
                                                                             gCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTTATTTGCCTTCTGAAAAAACTACCACAACAGGTTTGCACCGCCTGCCTGTTCCTCACCCACGACTTACGCCTGGTAGAACGCTTTTTGCCGCGGGTAATgg                                                  >  1:650485/1‑139 (MQ=255)
                                                                                   tCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCa                                                             <  1:3757763/122‑1 (MQ=255)
                                                                                     tAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGAc                                          <  1:1989834/139‑1 (MQ=255)
                                                                                      aaCCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACa                                         <  1:1563749/139‑1 (MQ=255)
                                                                                       aCCTTGATCTCGTGTTACAGGCGGAGGTGGTTTACCTGCGGAAAAAAATAAAACAAAACTGTGGGGCCGCCCTCCCGTTCATCACCAACACCTTACGCCAGGTGGAACGCCTTTTGCCGCGGGGAAGGGATATgttcaa                                        >  1:4590064/1‑134 (MQ=255)
                                                                                             aTCTCGTGTTACAGGCGGGTGGCAATTGCCTGCTGAAAAAACTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTAAGCCTGGTAGAACGCTTTTGCCAGCGGGTAATg                                                   >  1:1176283/1‑122 (MQ=255)
                                                                                                  gtgtTACAGGCGGGTGTCATTCGCCTGCTGAAAAAACTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGGAGAACGCTTTTTCCCGCGGGGAATGGTTTTGGACAACCGACAAAtcg                             >  1:607393/1‑139 (MQ=255)
                                                                                                   tgtTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTaa                                                     <  1:284776/114‑1 (MQ=255)
                                                                                                   tgtTACAGGCGGGTGTCATTCGCCCTCTGAAAAAGCTACCACAACAGTTTGGCACCGCCCTCCTGTTCCTCACCCACGACCTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGGTATGGACAACGGACAAAtcg                             >  1:2364183/1‑138 (MQ=255)
                                                                                                    gtTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCa                                                             <  1:34165/105‑1 (MQ=255)
                                                                                                    gtTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAAtcgtc                           <  1:1231328/139‑1 (MQ=255)
                                                                                                               gtgtCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACtt                                                                                      <  1:3609482/70‑1 (MQ=255)
                                                                                                               gtgtCATTCGCCTCCTGAAAAAGCTACAACACCAGTTTGGCGCCGCCTGCCTGTTCATCGCCCACGACGTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAgg                 <  1:469777/138‑1 (MQ=255)
                                                                                                                   cATCCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCCCCCACGACTTAc                                                                                    <  1:2278234/68‑1 (MQ=255)
                                                                                                                     ttCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTggg          >  1:2557405/1‑139 (MQ=255)
                                                                                                                         ccTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTa                                                                            >  1:1027773/1‑70 (MQ=255)
                                                                                                                         ccTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGGGGTGGgagag      >  1:2356297/1‑139 (MQ=255)
                                                                                                                             ctgAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGa                                                                          >  1:4857485/1‑68 (MQ=255)
                                                                                                                              tgAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGGAGAGAAAt  <  1:3781086/138‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
TAAGCGGCGGCCAGCTCCAGCGCGTCTGCCTGGCTCGCGCGCTGGCGGTCGAACCGAAACTACTGATTCTGGATGAAGCCGTTTCTAACCTTGATCTCGTGTTACAGGCGGGTGTCATTCGCCTGCTGAAAAAGCTACAACAACAGTTTGGCACCGCCTGCCTGTTCATCACCCACGACTTACGCCTGGTAGAACGCTTTTGCCAGCGGGTAATGGTTATGGACAACGGACAAATCGTCGAAACCCAGGTGGTGGGAGAGAAAT  >  CP009273/3611590‑3611853

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: