Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,6010CA42.0% 40.3 / 23.2 50L267F (TTG→TTTcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (17/12);  new base A (0/21);  total (17/33)
Fisher's exact test for biased strand distribution p-value = 7.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCG  >  NZ_CP009273/353464‑353727
                                                                                                                                          |                                                                                                                               
cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                 >  1:239843/1‑138 (MQ=255)
cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                 >  1:732452/1‑138 (MQ=255)
cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                 >  2:451805/1‑138 (MQ=255)
    tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                            >  1:226833/1‑139 (MQ=255)
     cgggCAATCCCGGCCGCAGCCCTGTGTCAGCGCCTCCCGTGATTCATTTCCCCCAACACCCGCGCATTTACCCACCCCCTTGTCCAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAAcc                                                                                                                           <  2:338724/137‑1 (MQ=255)
     cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                           >  1:606963/1‑139 (MQ=255)
     cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                           >  1:586260/1‑139 (MQ=255)
           ttcccggCCGCAGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCGCAAACACCCGCCCTTTTACCCAACCCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCGGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGtt                                                                                                                     <  1:660913/133‑1 (MQ=18)
             tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                   >  1:293012/1‑139 (MQ=255)
             tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                   >  1:364501/1‑139 (MQ=255)
                 gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                >  1:74918/1‑138 (MQ=255)
                 gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                >  1:812740/1‑138 (MQ=255)
                       gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                         >  2:90603/1‑139 (MQ=255)
                         ccggtgcCAGCGTCTCCCGTGATTCTTTCCCCCCAACACCCGCGCTTTTACCCAACCCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                       <  1:917033/136‑1 (MQ=255)
                               tCAGCGCCTCCCGTGATTAATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTAAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                 <  1:497272/139‑1 (MQ=21)
                                cggcgcCTCCCTTGATTCTTTCCCCCAAACACCCGCTCTTTTACCCACCGCCTTGTGAAAAAAGCCCTTCGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                <  1:554527/137‑1 (MQ=255)
                                cAGCGGCTCCCGTAATTCATTCCCCCCAACACCCGCGCATTTACCAACCGCCTGGCAAAAAAATGCCTCCGGGGCGGCTGTCGCCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                <  2:606963/139‑1 (MQ=255)
                                  gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcg                                                                                                >  1:774909/1‑136 (MQ=255)
                                  gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAag                                                                                              >  1:565404/1‑139 (MQ=255)
                                       tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAaggg                                                                                            >  1:139546/1‑134 (MQ=255)
                                       tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGaa                                                                                         >  1:94814/1‑138 (MQ=255)
                                       ccccGTGATTCTTTCCCCCCAACACCCGCGCTTTTATCAACCCCCATGGAAAAAAACGCCTCCGGGGCGCCTGTCTGCAGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                        <  2:732452/137‑1 (MQ=255)
                                         ccGTGTTCCATTCCCCCCAACACCCGCGCTTTTACCAAACCCAATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                       <  2:475195/138‑1 (MQ=255)
                                                  aTTCCCGCCAAAACCCGCCATTTTATCCAACGCCATGTAAAAAAATGCCTCCGCGGGGGCTGTCTGCCAGTTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                             <  1:382589/139‑1 (MQ=255)
                                                    tccccccAACACCCGCCCATTTACCCACCCCCATGTGAAAAATTCCCTCCGCGGGGGCTTTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                             <  1:451805/132‑1 (MQ=255)
                                                       caccAAAACCCGCGCTTTTATCCAACGCCAGGTGAAAAAATGCCGGCGGGGCGCCTGTCGGCCTGTTATTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                        <  2:35215/137‑1 (MQ=255)
                                                       cGCCAACACCCGCGCATTTATCCACCCCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCTCCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGAGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                        <  1:824931/139‑1 (MQ=255)
                                                       cGCCAACAACCGCGCTTTTACCCAACGCCATGCAAAAAAATGCCTCCGGGCCGGCGGTCTGCAGTCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                        <  2:289586/139‑1 (MQ=255)
                                                          caacacCCGCGCTTTTACCCAACGCAATGGCAAAAATTGCTGTCGGGGCGCTTGTCTGCCAGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                     <  2:586260/139‑1 (MQ=255)
                                                          caacacCCGCGCATTTACCCACCTCCATGTAAAAAAATGCTATCGCGGCGTCTGTCTCCCTCCTTTTTTACCGCCGCAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                     <  1:604302/139‑1 (MQ=255)
                                                                      cattACCCACCGCCGTGTTAAAAAATCCCTTCGCGCCGGCTGTCTGCCACCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                         <  1:96106/137‑1 (MQ=255)
                                                                         ttcccACCGCCATGTGCAAAATTCCCTCCGCGGCGGCTGTCCGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                       <  1:94533/135‑1 (MQ=255)
                                                                         tacccAACGCCATGTAAAAAAACCCCTGCGCGGCGTCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                       <  2:226833/135‑1 (MQ=255)
                                                                          aTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGGAAAAACTTTTAACCCCCCCTTTTTTGTGGGGATCgct                                                     >  1:720512/1‑139 (MQ=25)
                                                                           cccAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGATTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                    <  1:833828/138‑1 (MQ=255)
                                                                           cccAACCCCATGCACAAAATCCCCTTCGCGGGCGCGATCGGCCAGCTTTTTTCCCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                    <  2:293012/138‑1 (MQ=255)
                                                                                  ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAATTTTAGCCCGCCATGTTTTGTGGGGATCgccgctcggt                                              >  1:541730/1‑134 (MQ=37)
                                                                                         aaaaaaTGCCTTCGCGGCGGCTGTCTGCCATCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                       <  1:521710/137‑1 (MQ=255)
                                                                                                  ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTTAAACCCACATGTTGTGGTGGAGACGCTGCTGGGAAAAAAGAGGGAAgggg                              >  2:514868/1‑138 (MQ=255)
                                                                                                  ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGACGGGAAAAAAATTTTAGACCCCCATGTTGTGGGGGGATCGGTGCTGGGGAAAAAGTGGGAAggg                               >  2:71135/1‑138 (MQ=255)
                                                                                                      cgcgGCGGCTTTCGGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                            <  2:627512/137‑1 (MQ=255)
                                                                                                        cggcggCTGTCGGCCAGCTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                        <  1:720429/139‑1 (MQ=255)
                                                                                                                 tttGCCTGCTATTTTCCCCCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                 <  1:622715/135‑1 (MQ=255)
                                                                                                                 tcggCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                 <  2:191368/134‑1 (MQ=255)
                                                                                                                 tCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                 <  1:79097/137‑1 (MQ=255)
                                                                                                                 cccgCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                 <  2:929053/134‑1 (MQ=255)
                                                                                                                         ttATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                       <  1:337536/43‑1 (MQ=255)
                                                                                                                         cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                       >  2:337536/1‑44 (MQ=255)
                                                                                                                           tttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc      <  1:661759/137‑1 (MQ=255)
                                                                                                                           aTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc      <  1:125453/138‑1 (MQ=255)
                                                                                                                           aTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc      <  2:27789/138‑1 (MQ=255)
                                                                                                                           aTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc      <  2:619942/138‑1 (MQ=255)
                                                                                                                               tcccgcccCACCAAACCCGCCGTTCTCTCCAGTAGTGTCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg  <  1:643276/136‑1 (MQ=255)
                                                                                                                                          |                                                                                                                               
CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCG  >  NZ_CP009273/353464‑353727

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