Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,6010CA54.2% 39.4 / 98.9 72L267F (TTG→TTTcynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base C (24/9);  new base A (0/39);  total (24/48)
Fisher's exact test for biased strand distribution p-value = 4.85e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-07
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG  >  NZ_CP009273/353515‑353689
                                                                                       |                                                                                          
tttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                         >  2:201721/1‑90 (MQ=255)
     gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                   >  1:1050918/1‑90 (MQ=255)
     gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                   >  1:374234/1‑90 (MQ=255)
     cccAACACCCGCGCATTTACCCACCCCCATGTGCAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAAccgc                                                                                   <  2:1398016/89‑1 (MQ=18)
        aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                >  1:1510351/1‑90 (MQ=255)
        aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                >  2:412571/1‑90 (MQ=255)
          cacccGCGTTTTTATCCAACGCCATGTGAAAAAATGCTTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                              <  2:1536116/87‑1 (MQ=21)
          cacccGCGCTTTTCCCCACCCCCATGTGAAAAAATGCTTCCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                              <  2:1500902/87‑1 (MQ=16)
          cacccGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                              <  2:315297/87‑1 (MQ=21)
          cacccGCGCATTTACCCACCCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                              <  2:485970/87‑1 (MQ=21)
          cacccGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGGGGGGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTtc                                                                              <  2:428226/87‑1 (MQ=25)
          cacccGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                              <  2:736293/87‑1 (MQ=21)
          cacccGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTtc                                                                              <  2:242019/87‑1 (MQ=25)
          cacccGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                              <  2:495391/87‑1 (MQ=21)
          cacccGCCCTTTTACCCACCCCCATGTCCAAAAAGCCCTTCGCGGGGGTTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                              <  2:1440128/87‑1 (MQ=255)
          cacccGCCCTTTTACCCACCCCCATGTAAAAAAACGCCTCCGGGGGGGCTGTCGGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                              <  1:1537778/87‑1 (MQ=255)
          cacccGCACTTTTACCAAACCCCATGTAAAAAAAGCCCTCCGCGGCGGTTGTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                              <  1:1856702/87‑1 (MQ=255)
          caacCGCGCATTTACCCACCGCCATGTAAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                              <  2:1437928/90‑1 (MQ=21)
               gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                         >  2:1653567/1‑90 (MQ=255)
               gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                         >  1:72320/1‑90 (MQ=255)
                      tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                  >  1:667024/1‑90 (MQ=255)
                        cccACCGCCTTGTGAAAAAATGCTTCGGGGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                <  2:402309/89‑1 (MQ=255)
                        cccACCGCCATGTGAAAAATTCCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                <  2:563682/89‑1 (MQ=25)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                 <  1:886223/69‑1 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcccccg                                                                                 >  2:886223/1‑69 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                            >  1:604127/1‑90 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                            >  1:1593260/1‑90 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                            >  2:1206596/1‑90 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                            >  2:741729/1‑90 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                            >  2:803942/1‑90 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                            >  1:686628/1‑90 (MQ=255)
                            aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                            >  1:724591/1‑90 (MQ=255)
                              cccATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                          <  2:281471/89‑1 (MQ=37)
                                cATGTAAAAAACTCCCTTCGGGGGGGCTGTCTGCAATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                        <  2:851077/90‑1 (MQ=255)
                                 aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAaga                                                       >  2:1190617/1‑90 (MQ=255)
                                 aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAaga                                                       >  2:1695431/1‑90 (MQ=255)
                                     aaaaaaaTCCCTCGGGGGCGGTTTTCTGCCGTTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                  <  1:763022/88‑1 (MQ=255)
                                      aaaaaaTGCTTCCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                 <  2:1510351/89‑1 (MQ=255)
                                       aaaaaTGCTTCCGCGGGGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                                <  2:1658534/90‑1 (MQ=255)
                                       aaaaaTGCCTTTGGGGGGGTTTTCGCCCCAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                 <  1:585235/90‑1 (MQ=255)
                                         aaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGACGAGgaaaaaa                                              >  1:1747293/1‑86 (MQ=38)
                                                 ccggggcggtTGTCCCCCATTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                      <  2:1323485/86‑1 (MQ=255)
                                                 ccggggcggcTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                      <  2:511845/86‑1 (MQ=255)
                                                 ccgcgGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                      <  2:1691112/89‑1 (MQ=255)
                                                    gggcgccTGTCTCCCATTTTTTTTCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                   <  1:847630/89‑1 (MQ=255)
                                                      ggggcTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                 <  2:456396/88‑1 (MQ=255)
                                                      gcggcTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                 <  2:593424/90‑1 (MQ=255)
                                                         gttGTCCGCCAGTTTTTTTCCCCCCCCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                              <  2:1294855/88‑1 (MQ=255)
                                                         gcTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                              <  2:289988/90‑1 (MQ=255)
                                                                gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                              >  1:644791/1‑36 (MQ=255)
                                                                gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                              <  2:644791/36‑1 (MQ=255)
                                                                      tttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                  <  2:1278619/87‑1 (MQ=255)
                                                                        tttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  2:1210641/89‑1 (MQ=255)
                                                                        tttttcCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  1:657089/89‑1 (MQ=255)
                                                                        tttttcCCGCAACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  1:1459943/89‑1 (MQ=255)
                                                                        tttttcCCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  1:382352/89‑1 (MQ=255)
                                                                        tttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  2:40410/89‑1 (MQ=255)
                                                                        tttttcCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  1:657101/89‑1 (MQ=255)
                                                                         tttttCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  2:1030152/90‑1 (MQ=255)
                                                                         ttttcCCCCCCCAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                <  1:1027539/90‑1 (MQ=255)
                                                                          tttcCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg              <  2:1160422/90‑1 (MQ=255)
                                                                          tttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg              <  2:978433/90‑1 (MQ=255)
                                                                           ttcccgccgCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc             <  2:684476/87‑1 (MQ=255)
                                                                            tcccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct            <  2:526748/88‑1 (MQ=255)
                                                                               cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                >  1:942620/1‑36 (MQ=255)
                                                                               cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                <  2:942620/36‑1 (MQ=255)
                                                                               cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                >  1:1656728/1‑36 (MQ=255)
                                                                               cgccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                <  2:1656728/36‑1 (MQ=255)
                                                                                gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCCTCCTGTTTTGGGGGGATCGCTGCTgggt        >  2:336369/1‑88 (MQ=255)
                                                                                gccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAATTTTAAACCCTCCTGTTGTTTGGGGAACGCTGCTgggg        >  2:773900/1‑88 (MQ=255)
                                                                                gccgaAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa        <  1:156304/90‑1 (MQ=255)
                                                                                cccccaaaaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa        <  2:641814/86‑1 (MQ=255)
                                                                                 ccgcAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa       <  1:369692/90‑1 (MQ=255)
                                                                                    caaaaaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa    <  1:1314875/90‑1 (MQ=255)
                                                                                      acaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg  <  2:134471/90‑1 (MQ=255)
                                                                                      acaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg  <  2:714347/90‑1 (MQ=255)
                                                                                       |                                                                                          
TTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG  >  NZ_CP009273/353515‑353689

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: