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

AGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NZ_CP009273/4594034‑4594298
                                                                                                                                     |                                                                                                                                    
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                                                                                                                                   <  1:8955/65‑1 (MQ=255)
                                                                        cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                   >  2:8955/1‑65 (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)
                                                                                                                                     |                                                                                                                                    
AGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGC‑GATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAT  >  NZ_CP009273/4594034‑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: ‑