Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,807,8090GC35.0% 100.1 / 15.1 62E70Q (GAA→CAA) mcbCMicrocin B17‑processing protein McbC
Reads supporting (aligned to +/- strand):  ref base G (29/10);  new base C (0/21);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.09e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGGCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAACTAAAAGCCCT  >  ACB122/4807677‑4807940
                                                                                                                                     |                                                                                                                                   
aaaGGCAAATATACGAAACAAAAAGCAAATCAACACTTTAGCGCCCTGCACAAAATACCACCAAAAGACCTCTAAAAACTTACCGCAAACGAACTAAAACTAAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAAc                                                                                                                                  <  7:267924/136‑1 (MQ=255)
    gCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GACCCGAc                                                                                                                              >  7:240328/1‑136 (MQ=255)
     cAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTCGCCCCAAGCAAAAACAAA‑GAACCGAct                                                                                                                             >  8:76814/1‑136 (MQ=255)
     cAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGAct                                                                                                                             >  5:129022/1‑136 (MQ=255)
      aaaaataCCAAACACAAAGCAACTCAACACTTCCGCGCCCTCCACAAAATACAACCAAAACACCTCGAAAAACTTCCCGGAAACCAACAAAAACTCAAACCCATTGACCCCCCAACAAAAAAAAAC‑CACCCGActc                                                                                                                            <  6:71556/132‑1 (MQ=255)
                 acacAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTg                                                                                                                 >  6:245352/1‑136 (MQ=255)
                 acacAAACCAACTCACCACTTCAGCGCCCTCCACCAAAAACCACAAAAAGACCTCTAAAAACTAACCGTAAACCAACAAAAACTCAATACAATTGTACCCCCAACCAAAAACAAC‑CAACCGACTCTTACACCTCTg                                                                                                                 <  7:175242/136‑1 (MQ=255)
                   acaAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTgg                                                                                                                >  2:88046/1‑135 (MQ=255)
                     aaaGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGc                                                                                                               >  7:253322/1‑134 (MQ=255)
                     aaaGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGc                                                                                                               >  3:146401/1‑134 (MQ=255)
                      aaGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGa                                                                                                             >  1:118634/1‑135 (MQ=255)
                      aaGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAg                                                                                                            >  8:309178/1‑136 (MQ=255)
                          aaCTCAACACTTGAGCGCCCTGCACCAAATACAACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCca                                                                                                        >  8:175242/1‑136 (MQ=255)
                             tCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCcaccc                                                                                                     >  5:323547/1‑134 (MQ=255)
                               aaCACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACAccc                                                                                                   >  4:287870/1‑136 (MQ=255)
                                  acTTGCGCCCCCTCCACCAAAAACCACCAAACCACCTCTAAAAACTAACCGCAACCCAACAAAAACACAAAACCATTGTCCCCCCACCCAAAAACAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCt                                                                                                <  5:301454/136‑1 (MQ=255)
                                  acTTGACCGCCCTGCACCAAAAACCACCAAAACACCTCTAAAAAATTACCGCAAACCAAAAAAAACTCAAAACCATTGTAGCCCCAACAAAAAACAAA‑CACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCt                                                                                                <  5:119806/136‑1 (MQ=255)
                                  acTTCAGCGCCCTCCACAAAAACCCACCAAAACACCTCTAAAAACTTACCGTAAACCAACAAAAACACAATACCATTGTCGCCCCAACAAAAAACAAA‑CACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCt                                                                                                <  7:135464/136‑1 (MQ=255)
                                       aGCGCCGTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCAACCCCCGCTACg                                                                                             >  8:315219/1‑134 (MQ=255)
                                       aGCGCCCTCAACCAAAAACCACAAAAACACCTCAAAAAAATAACCGAAACCCAAAAAAAACTCAAAACCTTTG‑ACCCCCAAACAAAAACAAACCACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACga                                                                                            <  3:287870/135‑1 (MQ=255)
                                        gcgcACTACACAAAAAACAAACAAACGACCTCTAAAAAATTCAGTTAAACAAAATAAAAAACAATAAATTTCTAGCCCAACAAAAAAAAAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACgaag                                                                                          <  8:300396/136‑1 (MQ=255)
                                              tGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAg                                                                                    >  7:32758/1‑136 (MQ=255)
                                                 aCCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGtt                                                                                  >  7:53620/1‑135 (MQ=255)
                                                 aCCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGt                                                                                   >  1:337845/1‑134 (MQ=255)
                                                    aaaTACCACCAAACCACCTCAAAAAACTAACCGCAACCGAACAAAAACTCAATACCATTGTCCCCCCAACCAAAAACAAA‑CACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATg                                                                              <  5:161807/136‑1 (MQ=255)
                                                     aaTACCACCAAAACCCCTCCAAAAACTTACCGGAAACGAACTAAAACTCAAAACCATTGAACCCCCAAGCAAAAACAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGc                                                                             <  3:286512/136‑1 (MQ=255)
                                                     aaTACCACAAAACCACATCTAAAAAATAACCTTAACCCAAATAAAACTAATTACCTTTTTAGCCACAACAAAAAACAAA‑TAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGc                                                                             <  6:121687/136‑1 (MQ=255)
                                                        accaccAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAAAACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTAAACCTCTGGCGAGCACACACCCGCTACCAAGAAATAGTTCATGCCaa                                                                          >  5:121687/1‑136 (MQ=255)
                                                             caAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGCCTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGAtt                                                                     >  5:71360/1‑136 (MQ=255)
                                                             caAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGAtt                                                                     >  7:21431/1‑136 (MQ=255)
                                                             caAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCGACGAAGAAATAGTTCATGCGAACGAtt                                                                     >  4:60778/1‑136 (MQ=255)
                                                                  gACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAg                                                                >  6:104905/1‑136 (MQ=255)
                                                                   accccaaaaaaCAAACCCCAAACCAAAAAAAACAAAAAACCATTGACCCCCCAACCAAAAAAAAA‑CACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGc                                                               <  7:309178/130‑1 (MQ=255)
                                                                      tctaaaaaCTACCCGCAAACAAAATAAAACAAAATACCATTGTAGCCCAAACAAAAAACAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCAc                                                             <  1:88046/132‑1 (MQ=255)
                                                                       cGAAAAACTTACCGGAACCGAACTAACACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCCCACCCGCTACGAAGAAATAGTTCCTGCGAACGATTTGAAGCACGc                                                           >  7:199166/1‑136 (MQ=255)
                                                                       cGAAAAACTTACCGGAAACGAACGAAAACTCAATACCATTGTAGCCCCAAGCAACAACAAA‑GATCCGACGCTTACACCTCTGGCGAGCCCCCCCCCGCTACGAAGAAATAGTTCATGCGCACGCTTTGAAGCACGc                                                           >  7:204213/1‑136 (MQ=255)
                                                                          aacacatACCGCAAACGCACAAAAACTCCATACCATTTTCGCCCACACAAAAAAAAAC‑GAACCGACTCTTACACATCTGGCGAGCCCACCCCAGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct                                                        <  8:204213/130‑1 (MQ=255)
                                                                          aacaCTTACCGGAAACGAACAAAACCTCAATCCCTTTGACGCCCCAACCAAAAACAAC‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct                                                        <  8:199166/133‑1 (MQ=255)
                                                                          aaaaCTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct                                                        <  4:316286/136‑1 (MQ=255)
                                                                          aaaaCTTACCGCAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct                                                        <  4:119267/136‑1 (MQ=255)
                                                                          aaaaCTTACCCTAAAACAAAAAAAAATAAATACATTTCTCCCCCCAACAAAAAAAAAC‑CCCCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct                                                        <  6:323547/136‑1 (MQ=255)
                                                                          aaaaATAACCGTAAACAAAAAAAAAATAAATACAATTGTAGCCCCACCAAAAAACAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct                                                        <  8:128965/136‑1 (MQ=255)
                                                                           aaaCTAACCGCAAACCAACAAAACCAAAATCCCATTGTCGCCCCAACCAAAAAAAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCct                                                        <  1:111181/135‑1 (MQ=255)
                                                                              cTTACCTCAAACCAACAAAAACACAAAACAATTTTCCCCCACACAAAAAACAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTc                                                    <  5:117318/136‑1 (MQ=255)
                                                                               ttACCGGAAACGAAATAAAACTCAATACCATTGTAGCCCCAATCAAAAACAAA‑CAACCGACTCTTACACCTCCGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCa                                                   <  8:16499/136‑1 (MQ=255)
                                                                                   cGGAAACGAAATAAAACTCAATACCATTGTAGCCCCAAGGAAAAACAAA‑GAACCCACTCATACTCCTCTGGCGAGCCCCCACACGCTACGAAGAAAAAGTTCAGGCGAACTATTTGAAGCACGCCCTCTTCACtaa                                               >  5:107235/1‑134 (MQ=255)
                                                                                    ggAAACGAAATAAAACTCAATACCATTGTAGCCCCAAGCAAAACCAAC‑GACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg                                              <  8:240328/136‑1 (MQ=255)
                                                                                    gcaaaCCAACTAAAACTCAATACCATTGTAGCCCAAAGAAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg                                              <  8:21431/134‑1 (MQ=255)
                                                                                     gAAACTAACTAAAAATAAATACCATTGTAGCCCCAAGCAAAAACAAA‑TAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTg                                              <  2:98586/135‑1 (MQ=255)
                                                                                         cGAACTAAAACTCAATACAATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGAtata                                         <  4:146401/136‑1 (MQ=255)
                                                                                              tAAAACTCAAAAACATTGTCGCCCCAAGCAAAAACAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCa                                    <  8:323062/136‑1 (MQ=255)
                                                                                               aacaCTCAATCCCATTGTAGCCCCACGCACAAACAAC‑GACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCa                                    <  2:276814/132‑1 (MQ=255)
                                                                                               aaaaCCAAATACCATTGTACCCCCACCCAAAAACAAC‑GACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCa                                    <  1:309620/135‑1 (MQ=255)
                                                                                                aaaCTCAAAACCATGGTAGCCCCAACAAAAAACAAA‑CACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCaa                                   <  2:79693/135‑1 (MQ=255)
                                                                                                        tACCATTGTAGCCCCAAGCAAAAACAAA‑CAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCAttgctt                          <  6:61606/136‑1 (MQ=255)
                                                                                                         aCCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCAttgcttg                         >  3:302177/1‑136 (MQ=255)
                                                                                                                tAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCTCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGct                  >  4:79580/1‑136 (MQ=255)
                                                                                                                cacccccAACCAAAAACAAA‑CACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGct                  <  5:308046/133‑1 (MQ=255)
                                                                                                                   ccccAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGctct                >  8:47583/1‑135 (MQ=255)
                                                                                                                           aaaaaCAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACTCCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAACTa          >  1:36349/1‑133 (MQ=255)
                                                                                                                           aaaaaCAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTCATATAGTCCAAGCATTGCTTGCACCGCTCTTAACTa          >  1:281943/1‑133 (MQ=255)
                                                                                                                                cAAA‑CACCCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAACTAAAAGCCCt  <  7:142056/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
AAAGGCAACTATACGATACACAAAGCAACTCAACACTTGAGCGCCCTGCACCAAATACCACCAAAAGACCTCGAAAAACTTACCGGAAACGAACTAAAACTCAATACCATTGTAGCCCCAAGCAAAAACAAA‑GAACCGACTCTTACACCTCTGGCGAGCCCACACCCGCTACGAAGAAATAGTTCATGCGAACGATTTGAAGCACGCCCTCTTCACTTTGATATAGTCCAAGCATTGCTTGCACCGCTCTTAACTAAAAGCCCT  >  ACB122/4807677‑4807940

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