Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1610CT40.9% 45.7 / 20.4 45S62F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (9/17);  new base T (18/0);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 4.94e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NZ_CP009273/4594026‑4594298
                                                                                                                                       |                                                                                                                                          
gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttt                                                                                                                                         >  2:620058/1‑139 (MQ=255)
        agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCa                                                                                                                                 <  1:84200/139‑1 (MQ=255)
         gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCAt                                                                                                                                <  2:293490/139‑1 (MQ=255)
                 ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                                                                        <  2:247490/139‑1 (MQ=255)
                  gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGTGAAGt                                                                                                                         >  1:319890/1‑137 (MQ=255)
                  gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATgaggagg                                                                                                                         >  2:442356/1‑136 (MQ=255)
                  gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTGCATGGGAAGTgg                                                                                                                       >  1:295524/1‑138 (MQ=255)
                       gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTTCCATGAGAAGTGAgcg                                                                                                                    >  1:393012/1‑137 (MQ=255)
                            tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTGAGAAGTGGGCGCagagag                                                                                                             >  2:308725/1‑139 (MQ=255)
                             ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTGCCCTGAGAAGTGAGCGCAGAgggg                                                                                                            >  2:25867/1‑139 (MQ=255)
                               aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAggg                                                                                                           <  2:393012/138‑1 (MQ=255)
                                    gAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGTGGGCGCCGAGAGGGATTAtt                                                                                                     >  1:609403/1‑138 (MQ=255)
                                        gTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                 <  1:590691/139‑1 (MQ=255)
                                                          taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttttttttct                                                                                                                                     >  2:163274/1‑85 (MQ=255)
                                                          taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                     <  1:163274/85‑1 (MQ=255)
                                                          taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGGGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGC‑GAt                                                                               >  1:377418/1‑138 (MQ=255)
                                                          taAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GAtt                                                                              <  1:88667/139‑1 (MQ=255)
                                                               tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATTTTGCCTGTCTGAACTGGC‑GATTAAGtt                                                                         >  1:176545/1‑139 (MQ=255)
                                                               tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAAGGAGCGCAGGGAGGGATTATTTTTCCTGTCTGAATTGGC‑GATTAAGtt                                                                         >  2:155602/1‑139 (MQ=255)
                                                                 agtACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTc                                                                       <  1:626328/137‑1 (MQ=255)
                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttcttc                                                                                                                                   >  1:376681/1‑76 (MQ=255)
                                                                   tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                   <  2:376681/78‑1 (MQ=255)
                                                                        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTCTTTTTTTTCCCTGGGAAGTGAGGGCAGGGAGGGGTTTTCTTGCCTGGCTGGACTGGC‑GGTTAAATTTCCGCGTaa                                                                >  1:67994/1‑138 (MQ=255)
                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCATGAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTTCGCGTAcc                                                               >  2:154183/1‑139 (MQ=255)
                                                                         taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTAcc                                                               <  1:189545/139‑1 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGAGGGCAGAGAGGGATTATTTTTCCTGACTGAACTGGA‑GATTAAGGTTTCGCGCACCCgggg                                                          >  1:619512/1‑137 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGAGCGCCGAGAGGGGTTATTTTGTCTGACTGAAATGGC‑GATTAAGTTTCCGCGGACCCgggt                                                          >  2:626328/1‑137 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAATGGC‑GATTAAGTTTCCGCGTACCCGGCg                                                          >  2:143931/1‑139 (MQ=255)
                                                                              gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGGGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGG‑GATTAAGTTTCCGCGTACCCgggg                                                          >  1:622970/1‑137 (MQ=255)
                                                                                cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGGCATGGGAAGGGAGGGGAGAGGGGGGTTATCTTGGCTGACTGAAATGGC‑GGTTAAGGTTCCGCGGACACGGCggt                                                        >  1:164704/1‑137 (MQ=255)
                                                                                cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGC‑GGTTAAGTTTCCGCGTACCCGGCGtt                                                        >  2:158377/1‑139 (MQ=255)
                                                                                cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCCGAGGGGGGTTTTCTTTCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGtt                                                        >  2:211990/1‑139 (MQ=255)
                                                                                cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGGGCAGGGAGGGGGTTTTTTGCCTGACTGAAATGGC‑GATTAAGTTTCCGCGTACCCGGCGtt                                                        >  2:353809/1‑139 (MQ=255)
                                                                                cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                   >  2:8955/1‑65 (MQ=255)
                                                                                cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                   <  1:8955/65‑1 (MQ=255)
                                                                                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTTGCATGAGATGTGAGTGCAGCAAGGGGTCGTTTTGCCTGACCGAACTCGCGGATTAGGTGT‑CGCGGACTCGGTGTTTAg                                                     >  1:253693/1‑139 (MQ=255)
                                                                                        aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCCGAGAGGGATTATCTTGCCTGAATGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                >  1:432999/1‑139 (MQ=255)
                                                                                         tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTg                                               <  1:620058/139‑1 (MQ=255)
                                                                                            cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCg                                             <  1:476444/138‑1 (MQ=255)
                                                                                             aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGTGAGCGCAGAGGGGGGTTATTTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTTTTTCggg                                           >  1:526469/1‑138 (MQ=255)
                                                                                             aGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTATTTTTCTTCCATGAGAAGGGAGGGGAGAGAGGGGTTATTTTTCCTGACTGAACTGGC‑GATTAAATTTCCGCGTACCCGGCGGTTAGTTATTTCGga                                           >  1:617143/1‑139 (MQ=255)
                                                                                                  gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATCGAAAGTGGGCACCGAGAGGGCTTATTTTGGCTTCCCGGACTTGC‑GGATTAGTTTCCGCGCGACCGGCGTGTATGTGTTGGGggtgatg                                      >  1:388891/1‑139 (MQ=255)
                                                                                                  gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGGGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTAcc                                                               >  1:241624/1‑114 (MQ=255)
                                                                                                  gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTAcc                                                               <  2:241624/114‑1 (MQ=255)
                                                                                                      tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGatat                                  <  1:94840/139‑1 (MQ=255)
                                                                                                             gATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCt                           <  2:617143/139‑1 (MQ=255)
                                                                                                                  cGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATTGGAAGTGAGCGCAGAGGGGGGTTATTTTTCCTTGTTGAAATGGC‑GGTTTAGTTTCCGCGTAACCGGCGTTTTGTTTTTGCGGGTGGTGATATTGAAGCTCGGCt                      >  2:513412/1‑139 (MQ=255)
                                                                                                                                aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg         <  2:295524/138‑1 (MQ=255)
                                                                                                                                aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTg         <  2:319890/138‑1 (MQ=255)
                                                                                                                                       ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  1:353809/138‑1 (MQ=255)
                                                                                                                                       ctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt  <  2:643870/138‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
GCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NZ_CP009273/4594026‑4594298

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