Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273353,5910TC49.1% 80.2 / 12.4 60K271E (AAA→GAA) cynRtranscriptional regulator CynR
Reads supporting (aligned to +/- strand):  ref base T (25/5);  new base C (0/29);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 5.49e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTG  >  NZ_CP009273/353502‑353680
                                                                                           |                                                                                          
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                            >  1:962643/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                            >  1:741159/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                            >  2:1106696/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                            >  2:1268791/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                            >  2:198271/1‑90 (MQ=255)
cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGGCTGCCAGCTAttttt                                                                                            >  2:183889/1‑90 (MQ=255)
     gtgtTTCATTCCTGCCAACAGCCGCGCATTTATCCAACGCCATGTGCAAAAAGGCCTTCGCGGCGGCTGTCTGCCGGCTATTTTTccgcc                                                                                       <  2:940731/90‑1 (MQ=17)
     gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                       <  1:780847/90‑1 (MQ=255)
     gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                       <  1:242091/90‑1 (MQ=255)
     gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                       >  2:323115/1‑90 (MQ=255)
       gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                     >  2:351031/1‑90 (MQ=255)
       gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                     >  2:927936/1‑90 (MQ=255)
       gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                     >  1:328676/1‑90 (MQ=255)
             tttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                               >  2:1321712/1‑90 (MQ=255)
                  gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                         >  1:674060/1‑90 (MQ=255)
                     aacaacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                      >  1:317947/1‑90 (MQ=255)
                       cacccGCGTTTTTACCCACCGCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                    <  1:1106696/87‑1 (MQ=18)
                       cacccGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                    <  1:1003054/87‑1 (MQ=21)
                       cacccGCGCTTTTACCCACCGCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                    <  1:338644/87‑1 (MQ=21)
                       cacccGCGCTTTTACCCACCCCCATGTCAAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGCAACAAAACCGCCGTtc                                                                    <  1:1321712/87‑1 (MQ=19)
                       cacccGCGCTTTTACCCAACGCCAGGTGAAAAAATGCTTCCGCGGGGGCTGTCTCCCACTTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                    <  1:1258468/87‑1 (MQ=255)
                       cacccGCGCTTTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTTTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTtc                                                                    <  1:981640/87‑1 (MQ=16)
                       cacccGCGCATTTATCCACCGCCATGTGCAAAAATGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                    <  1:107961/87‑1 (MQ=21)
                       cacccGCGCATTTACCCACCCCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                    <  1:609776/87‑1 (MQ=21)
                       cacccGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTtc                                                                    <  1:1279174/87‑1 (MQ=21)
                       cacccGCCCATTTACCCAACCCCTTGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTAATTTTCCGCCGCAACAAAACCGCCGTtc                                                                    <  1:609727/87‑1 (MQ=25)
                       caacCGCGCTTTTCCCCACCCCCATGTAAAAAAATGCCTTCGCGGCGGTTTTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                    <  1:676825/90‑1 (MQ=16)
                       agcccGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTtc                                                                    <  1:431170/87‑1 (MQ=19)
                          ccGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTtctct                                                                 <  1:392862/90‑1 (MQ=21)
                            gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                               >  2:1167636/1‑90 (MQ=255)
                            gcgcATTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGTTTTCGCCAAGCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTcc                                                               <  1:32817/90‑1 (MQ=21)
                                   tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagct                                                        >  2:431170/1‑88 (MQ=255)
                                   tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                        >  1:568259/1‑90 (MQ=255)
                                   tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                        >  1:668002/1‑90 (MQ=255)
                                   tATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGGCTGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTTTCCAATAGTg                                                        >  2:913298/1‑90 (MQ=25)
                                    aTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagctg                                                       >  1:1084770/1‑87 (MQ=255)
                                       cAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                    >  1:378902/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  1:1310314/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  2:153820/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  1:295071/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  1:1244953/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  2:1328338/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  2:326281/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  2:801036/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                  >  1:436404/1‑90 (MQ=255)
                                         aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggcg                                                  >  2:524451/1‑88 (MQ=255)
                                           gCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg                                                >  2:613056/1‑88 (MQ=255)
                                           cccATGTGCAAAAATCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                <  1:230331/89‑1 (MQ=37)
                                           cccATGTAAAAAAATCCCTCCGCGGCGGCTGTCGGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGGGGGGCa                                                <  1:958626/89‑1 (MQ=255)
                                                gaaaaaaaaTCCCCCCGGGGCGTTTCTCTGCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                           <  2:1086158/86‑1 (MQ=255)
                                                   aaaaaaTGCCTCCGGGGGGGCTGTCGCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                        <  1:1202120/89‑1 (MQ=255)
                                                   aaaaaaTGCCTCCGCGGGGGTTGCCGGCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                        <  2:826072/89‑1 (MQ=255)
                                                   aaaaaaTGCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                        <  1:351031/89‑1 (MQ=38)
                                                    aaaaaTTCCTTTCGGGGGGGCTGTCGGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                        <  1:483559/90‑1 (MQ=255)
                                                    aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                       <  1:198271/90‑1 (MQ=255)
                                                    aaaaaTGCCTCCGGGGCGGTTTTCGGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                       <  1:316184/90‑1 (MQ=255)
                                                    aaaaaTCCCTCCGGGGCGTTTGTCGGCCAGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                       <  2:204367/90‑1 (MQ=255)
                                                     aaaaTGCCCTCGGGGGGGGTTGCCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa                                       <  1:647616/90‑1 (MQ=255)
                                                          gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAGGtttt                                 >  1:916849/1‑89 (MQ=255)
                                                            cTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                               <  1:518326/90‑1 (MQ=255)
                                                                ccggggcggcTTTCGGCCAGTTTTTTTGGCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                             <  2:505587/86‑1 (MQ=255)
                                                                     gcggcTGTCTGCATGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                        <  1:247529/90‑1 (MQ=255)
                                                                     gcggcTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                        <  1:476700/90‑1 (MQ=255)
                                                                        gttGTCCGCCAGTTTTTTTCCCCCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                     <  1:224851/88‑1 (MQ=255)
                                                                        gcTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                     <  1:54955/90‑1 (MQ=255)
                                                                            tcccccAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                 <  1:581748/86‑1 (MQ=255)
                                                                                      tttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc       <  1:664329/89‑1 (MQ=255)
                                                                                      tttttcCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc       <  2:610531/89‑1 (MQ=255)
                                                                                      tttttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc       <  1:425647/89‑1 (MQ=255)
                                                                                       tttttCCCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc       <  1:1112938/90‑1 (MQ=255)
                                                                                        tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg     <  1:339108/90‑1 (MQ=255)
                                                                                           tCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTTAACCCCCCCTGTTTTGTGGGGAACGctggtg  >  1:278264/1‑90 (MQ=255)
                                                                                           |                                                                                          
CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTG  >  NZ_CP009273/353502‑353680

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