Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0122.G22.4% 101.3 / 44.7 49intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (24/14);  new base G (7/4);  total (31/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 = 7.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACAC  >  NC_000913/379902‑380140
                                                                                                                 |                                                                                                                                
agacagataCTATAGGGGGGTATTATATATGTCAATGCAAACCCCCCTATAGTATAAAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                       <  1:988114/133‑1 (MQ=255)
                 gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGGAGGGAAGGCCGCCCCGGAGAgcggg                                                                                      >  2:13052/1‑139 (MQ=255)
                 gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAtctgtct                                                                                      >  2:988114/1‑133 (MQ=255)
                           tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGTCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAgg                                                                            >  2:72330/1‑139 (MQ=255)
                              gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt                                                                         >  1:622169/1‑139 (MQ=255)
                                 aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                                                       >  2:913201/1‑139 (MQ=255)
                                 aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAg                                                                      >  1:912643/1‑139 (MQ=255)
                                   tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATGGGCGGGGAAGGCCGCCCCGGATAGGGGGGGGCAGAAGGAATAAGaa                                                                    >  1:941987/1‑139 (MQ=255)
                                   tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGCGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                                                    >  2:1074413/1‑139 (MQ=255)
                                        aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                               >  1:661729/1‑139 (MQ=255)
                                                  aGTATAAAAAAATGGAGAAAAAAATTTATCCGGGGAAAGGGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACg                                                     <  2:941987/139‑1 (MQ=255)
                                                      tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGcg                                                 >  1:951687/1‑138 (MQ=255)
                                                      tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                 >  2:296072/1‑139 (MQ=255)
                                                      tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAAACGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt                                                 >  1:420242/1‑139 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACACAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAGGGCCGCCCCGGATAGCGGGCGGGGGAAGGAATCAGAATTCCCAGGTCAGACGGGCGGCAAg                                            >  1:464170/1‑139 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCAAATAGGCCGGGGAGGGCGCACCGGAAAGCGGGCGGCAGAAGGGAACGAAAGTTCCAGGTCAAACgggagggca                                             >  1:745953/1‑133 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGccc                                                                                                 >  1:205687/1‑86 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGGAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                            >  2:262478/1‑139 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                            >  1:719278/1‑139 (MQ=255)
                                                           cacTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGccc                                                                                                 <  2:205687/86‑1 (MQ=255)
                                                              tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                                                       <  2:919024/94‑1 (MQ=255)
                                                              tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                                                       >  1:919024/1‑94 (MQ=255)
                                                              tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                       >  1:589143/1‑77 (MQ=255)
                                                              tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                       <  2:589143/77‑1 (MQ=255)
                                                              tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                         >  2:860518/1‑139 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                        >  2:983505/1‑139 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                       >  2:1008518/1‑139 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                       <  2:661729/139‑1 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                       <  1:860518/139‑1 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGAAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTgaa                                       >  1:754932/1‑137 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg                                      >  2:370736/1‑139 (MQ=255)
                                                                gagaATAAAATTTATCCGGTGAATGGGGGGGCAACAAAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                       <  1:262478/139‑1 (MQ=255)
                                                                 ggaAAAAAATTTTTCCGGGGAATGTGGGGGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                       <  1:185165/138‑1 (MQ=255)
                                                                 agaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                     <  1:370736/139‑1 (MQ=255)
                                                                   aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AATGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGGTAGCGGACGGCAGAAGGGATGAGAATTTCCAGGTCAGACGGGCTGCAAGTTGAAGa                                     >  1:819603/1‑139 (MQ=255)
                                                                     tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGAcc                                   >  1:681499/1‑139 (MQ=255)
                                                                             aTCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCGGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaaac                           >  2:185165/1‑138 (MQ=255)
                                                                               ccGGTGAATGTGGTCGGAAAACAAAGGGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGaa                                                                              <  1:545354/85‑1 (MQ=255)
                                                                               ccGGTGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGaa                                                                              >  2:545354/1‑85 (MQ=255)
                                                                                  gggAAGGTGGTCGGAAAACAAAGAGGAGAG‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCgg                     <  2:464170/137‑1 (MQ=255)
                                                                                  gTGAATGTGGTCGGAAAACAAAGAGGA‑AA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcatc                       >  2:49681/1‑136 (MQ=255)
                                                                                  gTGAATGTGGTCGGAAAACAAAGAGGA‑AA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcatc                       <  1:49681/136‑1 (MQ=255)
                                                                                   tGAATGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATGGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGt                    >  1:27218/1‑139 (MQ=255)
                                                                                     aaTGTGGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTg                  <  1:1074413/139‑1 (MQ=255)
                                                                                       tgtgGTCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                <  2:622169/139‑1 (MQ=255)
                                                                                            tCGGAAATCAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                              >  1:772784/1‑104 (MQ=255)
                                                                                            tCGGAAAACAAAGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                              <  2:772784/104‑1 (MQ=255)
                                                                                               gAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGa                                                                           <  2:556873/73‑1 (MQ=255)
                                                                                               gAAAACAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGa                                                                           >  1:556873/1‑73 (MQ=255)
                                                                                                 aaaCAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc       <  2:681499/139‑1 (MQ=255)
                                                                                                 aaaCAAAGAGGA‑AAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc       <  1:913201/139‑1 (MQ=255)
                                                                                                     aaaGAGGA‑AAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCacac  <  2:719278/139‑1 (MQ=255)
                                                                                                            a‑aa‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                               >  2:674360/1‑86 (MQ=255)
                                                                                                            a‑aa‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                               <  1:674360/86‑1 (MQ=255)
                                                                                                                 |                                                                                                                                
TGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGA‑AA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACAC  >  NC_000913/379902‑380140

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑