Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913445,5580AC27.8% 65.3 / 13.5 36V370G (GTA→GGA) yajRputative major facilitator superfamily transport protein YajR
Reads supporting (aligned to +/- strand):  ref base A (19/7);  new base C (0/10);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 7.65e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAA  >  NC_000913/445424‑445690
                                                                                                                                      |                                                                                                                                    
aaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGGACTGGCCACTGTCAGCCACCCAGCGGCCAGCCTTTCGCCAACGAGAAAAAccccc                                                                                                                                 >  2:1258088/1‑140 (MQ=255)
 aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                <  2:1056118/140‑1 (MQ=255)
 aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCt                                                                                                                                <  2:556059/140‑1 (MQ=255)
  cGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTg                                                                                                                               <  2:52262/140‑1 (MQ=255)
              cccTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                   <  1:1258088/139‑1 (MQ=255)
                 tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGCGACTGTCAGCCACACCGCGGCCCGCATTGCGCCCGCGAGAAAAAACCCCCGACCCTCAAACATGcc                                                                                                                >  2:321054/1‑140 (MQ=255)
                   gTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGCTCTTTCTTGGTACTGGCGACCGTCCGCCACACCGCGGCCCGCATTGCGCCAGCGAGAAATACCCCCCGACCGGCAAACATGCCgg                                                                                                              >  2:731557/1‑139 (MQ=255)
                    tttGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGtt                                                                                                             <  1:353254/140‑1 (MQ=255)
                           cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                      <  1:894555/140‑1 (MQ=255)
                            aaTGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTTATGGTACTGGGGACCGTCACCCACACTGCGGCCAGCATTGCGCCAGCGAGAAAAACCCCCTGACCGTCAAACATGCCGTTAATCCaac                                                                                                     >  2:1176075/1‑138 (MQ=255)
                                                  cgcAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgc                                                                               <  1:894520/140‑1 (MQ=255)
                                                    aaacctgctgACATAGGGCGGTTCTTTCAGGGTACTGGCGCCGGTCAGCCACACGGCGGCCAGCATGGCGCGAGCGAAAAATCCCCCCTGACCTTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                             <  2:636740/136‑1 (MQ=255)
                                                          gctgACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTTCGCCAGCGAGAAATACCCCCCGACCGTCAAACATGCCGGTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCaa                                                                       >  2:88718/1‑140 (MQ=255)
                                                              aCATCCGGCGTTTCTTTCATGGTCCTGGCGACTGTCAGCCCCACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGaaa                                                                   <  2:1268204/140‑1 (MQ=255)
                                                                 tACGGCGGTTTTTTTATGGGACTGGCGACTGGCAGCCACACTGCGGCCAGCATTTCGCCCAGGAGAAAAACCCCCCGACCGTCCAACATGCCCTTTATCCCGCCGCCCAGCGAACCGCCCATCCCCCCGCCCAGaaaaag                                                                >  2:607739/1‑137 (MQ=255)
                                                                   cggcggTTCTTTTATGGGACTGGCGACTGTCAGCCACACTGCGGCCAGCATTTCGCCAGCGAGAAAAAACCCCCGACCGTCCAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCCATCGCCACCCCCAGAAActggc                                                              >  2:1057162/1‑140 (MQ=255)
                                                                       ggTTCTTTCATGGTACGGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                          <  1:257465/140‑1 (MQ=255)
                                                                       gccattttCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCCCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggt                                                          <  2:426993/135‑1 (MQ=255)
                                                                        gTTCTTTCATGGTACGGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCGGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                         <  1:1176075/140‑1 (MQ=255)
                                                                        catttttCATGGTTCGGGCGCCTGTCGCCCACACTGCGGCCAGCTTGGCGCCGGGGAAAAATCCCCCCTGACCGTAAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                         <  2:1130955/136‑1 (MQ=255)
                                                                         ttCTTTTATGGGACTGGCGACTGTTCGCCACACTTCCGCCAGCCTTGCGCCAGCGAGGAAAACCCCCCCACCGTCAAACATGCCGTTTATCCAGCCCCCCCGCGAACCCCCCATCGCCACGCCCAGAAAACGGCtggggg                                                        >  2:47140/1‑140 (MQ=255)
                                                                         ttCTTTCATGGTACTGGCGCCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtgg                                                        <  1:321054/140‑1 (MQ=255)
                                                                               cATGGTACGGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTaaa                                                  <  1:1057162/140‑1 (MQ=255)
                                                                                   gTACGGGCGACTGTCCGCCACACGGCGGCCAGCTTTGCGCCAGGGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                              <  2:77808/140‑1 (MQ=255)
                                                                                     aCTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                            >  1:552854/1‑140 (MQ=255)
                                                                                     aCTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                            >  1:543708/1‑140 (MQ=255)
                                                                                      cTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAccctg                                           >  1:1130955/1‑138 (MQ=255)
                                                                                      cTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCAt                                           >  1:916174/1‑140 (MQ=255)
                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcctgtc                                          >  1:426993/1‑136 (MQ=255)
                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                          >  2:1129198/1‑140 (MQ=255)
                                                                                       tGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATc                                          >  1:1482554/1‑140 (MQ=255)
                                                                                           gccTGTCACCCACACTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                      <  2:199679/138‑1 (MQ=255)
                                                                                            ccTGCCACCCACCCTGCGGCCAGCTTTGCGCCGGGGAAAAATCCCCCCTGCCCGCAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                     <  2:1420727/139‑1 (MQ=255)
                                                                                            aCTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                     >  1:199679/1‑140 (MQ=255)
                                                                                             cTGTCACCCCCCCGGCGGCCAGCGTTGCGCCAGCGAAAAATACCCCCTGCCCGCAAAACATGCCGTTAACCCAGCCGCCCAGCGACCCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                    <  2:1461341/140‑1 (MQ=255)
                                                                                              tgccAGCCCCCCTGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCTGCCCGTCAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                   <  2:552854/137‑1 (MQ=255)
                                                                                              tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTTGAGTAAACACCCATCGCCGTAc                                   >  2:1391733/1‑140 (MQ=255)
                                                                                              tGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                   >  1:1420727/1‑140 (MQ=255)
                                                                                                tCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCt                                 >  2:762045/1‑140 (MQ=255)
                                                                                                   gCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTg                              >  1:1461341/1‑140 (MQ=255)
                                                                                                   gCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTg                              >  1:387442/1‑140 (MQ=255)
                                                                                                          tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                       >  1:1204094/1‑140 (MQ=255)
                                                                                                            cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                     >  2:278025/1‑140 (MQ=255)
                                                                                                            cGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                     >  2:351461/1‑140 (MQ=255)
                                                                                                                 agcttggcgcCGGCGAAAAACCCCCCCTGACCGCAAAACATGCCGTAAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                <  2:916174/134‑1 (MQ=255)
                                                                                                                 aGCTTTGCGCCGGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                <  2:144531/140‑1 (MQ=255)
                                                                                                                 aGCATGGCGCCAGCGAGAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                <  1:498091/140‑1 (MQ=255)
                                                                                                                  gCATTGCGCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAc               <  1:793596/140‑1 (MQ=255)
                                                                                                                   aattcggcCAGGGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAACCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct              <  2:1309162/134‑1 (MQ=255)
                                                                                                                     ttGCGCCGGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct            <  1:816716/140‑1 (MQ=255)
                                                                                                                     taggcccAGAGAAAAACACCCCCGAACCGTCAACCATGCCGTTAACCCAGCGGCCCGGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct            <  2:1214667/135‑1 (MQ=255)
                                                                                                                        cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                >  2:874844/1‑133 (MQ=255)
                                                                                                                        cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                <  1:874844/133‑1 (MQ=255)
                                                                                                                        cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa         >  1:393214/1‑140 (MQ=255)
                                                                                                                        cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa         >  1:752883/1‑140 (MQ=255)
                                                                                                                               gagaAATCCCCCCGGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATaa  <  1:808174/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
AACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAA  >  NC_000913/445424‑445690

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