Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913445,5580AC26.2% 87.3 / 14.9 43V370G (GTA→GGA) yajRputative major facilitator superfamily transport protein YajR
Reads supporting (aligned to +/- strand):  ref base A (25/6);  new base C (0/11);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 2.89e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGG  >  NC_000913/445423‑445696
                                                                                                                                        |                                                                                                                                          
aaaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTTCGCCAGCGAGAAATAcccc                                                                                                                                        >  1:1194747/1‑140 (MQ=255)
aaaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAACCCCACTGCGGCCAGCATTTCGCCCGCGAGAAATAcccc                                                                                                                                        >  1:618914/1‑140 (MQ=255)
 aaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAAAAccccc                                                                                                                                       >  1:482957/1‑140 (MQ=255)
  aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCCGCATTTCGCCAGCGAGAAATAcccccc                                                                                                                                      >  2:1312968/1‑139 (MQ=255)
   cGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTg                                                                                                                                     <  2:112571/140‑1 (MQ=255)
    gCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGa                                                                                                                                    <  2:95397/140‑1 (MQ=255)
              cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                          <  1:890632/140‑1 (MQ=255)
               gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                         <  1:127921/140‑1 (MQ=255)
                ccTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGAGTGGCACCCACACTTCTGCCAGCATTTCGGCCGCGGGGAAAAACCCCTGCCCCTCAAACacg                                                                                                                        >  2:596920/1‑138 (MQ=255)
                  tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAAAATGcc                                                                                                                      >  1:188374/1‑140 (MQ=255)
                    gTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGt                                                                                                                    >  1:498423/1‑140 (MQ=255)
                           gCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGGGACCGTCAGCCACACTGCGGCCAGCATTGCGCCCGCGAGAAAAAACCCCTGACCGTCAAACATGCCGTTAATccc                                                                                                             >  2:936723/1‑139 (MQ=255)
                            cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATGCCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                            <  1:596920/140‑1 (MQ=255)
                                      cGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGGACTGGCGACTGTCAGCCACACTGCGGCCCGCATTGCGCCAGCGAGAAAAACCCCCTGACCCGCAAACATGCCGTTAAACCAGCCGCCCcgcg                                                                                                  >  2:263310/1‑140 (MQ=255)
                                                    gcAAACTGCTGACATACGGCGGTTCTTTTA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTTCGCCAGCGCGAAAAACCCCCTGGCCGTCAAAAAAGCCGTTAATCCAGCCCCCCCGCGAACCGCCCAAcgcc                                                                                    >  2:991619/1‑140 (MQ=255)
                                                               aCATCCGGCGGTTTTTTCA‑TGGTCCTGGCGACTGCCAGCCACCCTGCGGCCAGCATTGCGCAAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGaaa                                                                         <  2:498423/140‑1 (MQ=255)
                                                                    ggggggTTTTTTCA‑TGGTCCTGGCGCCTGTCAGCAACACTGGGGCCAGCATGGCGCCAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggc                                                                    <  2:1048822/136‑1 (MQ=255)
                                                                    cggcggTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggc                                                                    <  1:962045/140‑1 (MQ=255)
                                                                    cggcggTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCCTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACTTCCCATCGCCACGCCAAGAAActagc                                                                    >  2:1255964/1‑140 (MQ=255)
                                                                        ggTTCTTTCA‑TGGTACGGGCGACTGCCAGCCACCCTGGGGCCAGCATGGCCCCAGCGAGAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                                <  2:188374/140‑1 (MQ=255)
                                                                         gTTTTTTCATTGGTACGGGCGGCTTTCACCCACGCCGCGGCCGGCGTCGCGCCGGCGGGAAATACCCCCTGACGCACAATCATGCCGTTAATCCGGCCCCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                                <  2:166412/140‑1 (MQ=255)
                                                                         gTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                               >  1:587737/1‑140 (MQ=255)
                                                                          ttCTTTCA‑TGGTACTGGCGACTGTTAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAGCCGCCAATCGCCACGCCAAGAAACTGGCtggtgg                                                              <  2:1150738/140‑1 (MQ=255)
                                                                                    ggTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAacac                                                     >  2:653676/1‑140 (MQ=255)
                                                                                     gTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                    >  1:168280/1‑140 (MQ=255)
                                                                                     gTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                    >  2:598025/1‑140 (MQ=255)
                                                                                       aCTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                                  >  2:1299634/1‑140 (MQ=255)
                                                                                        cTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                                 >  1:838757/1‑140 (MQ=255)
                                                                                          ggCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                               >  1:962092/1‑140 (MQ=255)
                                                                                          ggCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                               <  1:263310/140‑1 (MQ=255)
                                                                                           gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaaa                                                        <  1:393458/130‑1 (MQ=255)
                                                                                           gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                              <  1:1342569/140‑1 (MQ=255)
                                                                                           gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaaa                                                        >  2:393458/1‑130 (MQ=255)
                                                                                           gCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                              >  2:972067/1‑140 (MQ=255)
                                                                                            cGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGcc                                             >  1:1105798/1‑140 (MQ=255)
                                                                                                tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                         >  1:444819/1‑140 (MQ=255)
                                                                                                tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                         >  1:444719/1‑140 (MQ=255)
                                                                                                tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                         >  1:793371/1‑140 (MQ=255)
                                                                                                tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                         >  1:95336/1‑140 (MQ=255)
                                                                                                tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                         >  2:1056251/1‑140 (MQ=255)
                                                                                                tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                         >  1:159555/1‑140 (MQ=255)
                                                                                                       cacaCTGCGGCCAGCATGGCGCCAGGGAGAAATCCCCCCTGACCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTa                                  <  2:618914/140‑1 (MQ=255)
                                                                                                            tGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCGGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                             <  1:1312968/140‑1 (MQ=255)
                                                                                                              cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                           >  1:1401095/1‑140 (MQ=255)
                                                                                                              cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                           >  2:714741/1‑140 (MQ=255)
                                                                                                                gCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTgg                         >  2:248419/1‑140 (MQ=255)
                                                                                                                   ggctttGCCCCAGGGAAAAATCCCCCCGCCCCGTAAACCATGCCGTTAATCCACCCCCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                      <  2:162797/136‑1 (MQ=255)
                                                                                                                   aGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                      <  1:631449/140‑1 (MQ=255)
                                                                                                                     cATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                    <  2:838757/140‑1 (MQ=255)
                                                                                                                         gcgcCAGCGAGAAACCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt                <  2:95336/140‑1 (MQ=255)
                                                                                                                            caggCGAAAAATCCCCCCGGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACt             <  2:962092/137‑1 (MQ=255)
                                                                                                                              aGCGAGAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTCACTGa           <  1:714741/140‑1 (MQ=255)
                                                                                                                              aGCGAAAAATCCCCCCTGACCGAAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGa           <  1:12313/140‑1 (MQ=255)
                                                                                                                                   gaAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGt      <  1:1299634/140‑1 (MQ=255)
                                                                                                                                   aaaaTCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGt      <  1:936723/139‑1 (MQ=255)
                                                                                                                                    aaaTACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTg     >  2:878985/1‑140 (MQ=255)
                                                                                                                                     aaTCCCCCCTGCCCGTAAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGa    <  2:1105798/140‑1 (MQ=255)
                                                                                                                                       tACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAgg  >  1:684143/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                          
AAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCA‑TGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGG  >  NC_000913/445423‑445696

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